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REVIEW ARTICLE 


THE NATURE AND PROPERTIES OF TRANSITION METAL HYDRIDES 


G. G. LIBOWITZ 


Atomics International, A Division of North American Aviation, Inc., Canoga Park, California, USA 


The basic chemical and physical properties of transition 
metal hydrides are reviewed with particular emphasis 
on the hydride phases rather than dilute solutions of 
hydrogen in metals. The review includes discussions 
of the phases present in metal-hydrogen systems, their 
crystal structures, pressure-composition-temperature 
relations, thermodynamic properties, and electrical and 
magnetic properties. 

The nature of transition metal hydrides is considered. 
Evidence is given to show that they are definite 
chemical compounds (rather than solid solutions) 
having some degree of ionic character. The large 
deviations from stoichiometry exhibited by these 
compounds can be explained on the basis of lattice 
defects. 


Les propriétés chimiques et physiques fondamentales 
des hydrures des métaux de transition sont passées en 
revue en insistant particuliérement sur les phases 
hydrures plutot que sur les solutions diluées d’hydro- 
géne dans les métaux. Cette mise au point discute des 
phases présentes dans les systémes métal-hydrogéne, 
de leurs structures cristallines, des relations entre 
pression, composition et température, des propriétés 
thermodynamiques et des propriétés électriques et 
magnétiques. 


La nature des hydrures des métaux de transition est 
examinée. Il est démontré que ce sont des composés 
chimiques définis (pluté6t que des solutions solides) 
possédant un certain caractére ionique. Les grands 
écarts & la stoechiométrie que présentent ces composés 
peuvent étre expliqués en se basant sur les défauts 
réticulaires. 


Es werden die grundlegenden chemischen und physi- 
kalischen Eigenschaften von Hydriden der Uber- 
gangsmetalle besprochen, wobei grosserer Wert auf 
Hydridphasen als auf feste L6sungen von Wasserstoff 
in Metallen gelegt wird. Der gegebene Uberblick 
enthalt Erérterungen tiber die Phases in Metall- 
Wasserstoff-Systemen, uber ihre Struktur, tber 
Zusammenhinge zwischen Druck, Konzentration und 
Temperatur, ttber thermodynamische Eigenschaften 
und iiber elektrische und magnetische Eigenschaften. 

Der Aufbau der Hydride der Ubergangsmetalle wird 
betrachtet. Dabei werden Hinweise dafiir angegeben, 
dass es sich um definierte chemische Verbindungen 
(und nicht um feste Loésungen) mit dem Charakter 
ahnlich Ionenkristallen handelt. Die grossen Abwei- 
chungen dieser Verbindungen von der stdchiome- 
trischen Zusammensetzung kénnen durch Gitterfehler 
erklart werden. 
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1. Introduction 


Since hydrogen is theoretically the most 
effective thermal neutron moderator, it becomes 
highly desirable to find materials which contain 
a high concentration of hydrogen atoms per 
unit volume. This objective can be accomplished 
most efficiently through the utilization of 
metallic hydrides. Therefore, the hydrides of 
the transition metals have been the subject of a 
great deal of study since the advent of nuclear 
power. Table 1 shows the number of hydrogen 


TABLE 1 


Ny Values for hydrogen containing compounds 


Compound Nu X 10-22 (atoms/cc) 
Liq. He (— 253° C) 4.2 
Liq. CHa (— 164° C) 6.3 
H20 | 6.7 
Li 5.9 
ZrHe 7.2 
UHs 8.2 
Tie 9.5 


atoms per cubic centimeter, Ny, for some more 

familiar hydrides as compared with water, 

liquid hydrogen, and a saturated organic 

molecule. 

ee es of Hydride) 
M.W. of hydride 

molecule) x 6.02 x 1023. 


The Ny values shown for the hydrides and 
water are at room temperature. Although Ny 
for the transition metal hydrides are not 
significantly greater than for liquid hydrogen 
or water, the value of the hydrides as nuclear 
moderators lies in their stability at higher 
temperatures. Whereas water and organic 
molecules vaporize (or decompose) at higher 
temperatures, the Ny values for the transition 
metal hydrides do not change appreciably until 
their dissociation temperature (usually several 
hundred degrees) is reached. From the stand- 
point of a moderating material, a useful hydride 
should have a high hydrogen to metal ratio, high 
density, low molecular weight, and a low 
dissociation pressure. 

It is the purpose of this paper first to review 


(No. of H atoms per 


the basic chemical and physical properties of 
transition metal hydrides (and related hydrides 
of interest) which could lead to an elucidation 
of the fundamental nature of these materials. 
This discussion will cover the phases present in 
metal-hydrogen systems, their crystal struc- 
tures, pressure-composition-temperature rela- 
tions, thermodynamic properties of hydrides, 
and, where available, their electrical and mag- 
netic properties. Secondly, some ideas on the 
bonding and nature of these substances will be 
presented, and correlated with experimentally 
determined properties. 

In 1947, Smith 1) presented a rather complete 
review on occlusion of hydrogen by metals 
(predominantly the transition metals), and in 
1952, Hurd 2) published a book on the chemistry 
of hydrides in which some space was devoted 
to a review of hydrides of the transition metals. 
Therefore, in this paper, emphasis will be placed 
on work performed subsequent to these dates, 
and also on the hydride phases rather than 
dilute solutions of hydrogen in metals. 


2. Properties of transition metal hydrides 


2.1. HypDRIDES OF THE GROUP IV A METALS 


Titanium Hydride 

Phase Relations: A detailed study of the 
pressure-composition-temperature relations in 
the titanium-hydrogen system was first carried 
out by McQuillan 3). The phase diagram he 
obtained from this work is essentially the one 
shown in Fig. 1 with the dotted line representing 
McQuillan’s extrapolations of his measured data 
(full line). Presumably, as hydrogen is dissolved 
in titanium metal, the « to # transformation 
temperature is lowered giving rise to the two 
phase region («+) observed at the lower 
hydrogen compositions (<0.7 H/Ti) in the 
pressure-composition isotherms. Unfortunately, 
there have been no X-ray studies of this system 
at high temperatures and low hydrogen contents 
to establish if the second phase is actually the 
high temperature body centered cubic 6 modi- 
fication of titanium metal. At higher hydrogen 
contents, a y or hydride phase is formed which 
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Titanium-hydrogen phase diagram. 


e McQuillan 
® Haag and Shipko (Ti-D) 
Fig. 1. 


approaches the stoichiometric formula TiH»- 

Fig. 1 also includes the points obtained from 
the subsequent P-C-7' study of this system by 
_ Haag and Shipko 4). These authors also studied 
the titanium-deuterium and titanium-tritium 
systems and observed no isotopic effect with 
respect to either dissociation pressures or phase 
boundaries. Except for the boundary between 
the 6+¥y phases and the y phase, the agreement 
with McQuillan’s work seems to be quite good. 
Gibb, McSharry and Bragdon*) report an 
additional phase at close to stoichiometric TiHs 
from a study of pressure-composition isotherms. 
Although it is shown as 6 in fig. 1, the existence 
of this phase is questionable since the authors 
calculated their compositions on the basis of 
their starting material being stoichiometric 
TiHs, and no actual experimental points were 
obtained on the two phase plateau regions of 
the isotherm. 


Thermodynamic Properties: The heat 
of formation of the nonstoichiometric hydride 
present in the two-phase region at 400—500° C 
(TiH;,4) is —22.6 kcal/mole as calculated from 
the data of Haag and Shipko. 

The dissociation pressure of the hydride to 
6-Ti and hydrogen is given by the relation f 


+ Except where explicitly given this form of 
relation will be used and the values of A and B 
given for the phase under discussion. 
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log Pip (mm) = a ea 
where A=8630 and B=12.43. These values 
were calculated from the data of McQuillan. 


Crystal Structure: An X-ray structure 
determination of the hydride phase by Hagg ®) 
revealed that the titanium atoms were in a 
face-centered cubic structure with a lattice 
parameter ranging from 4.40 to 4.46 A with 
increasing hydrogen content. This structure was 
confirmed in an X-ray and neutron diffraction 
study by Sidhu, Heston and Zauberis ’). The 
neutron diffraction results on titanium deuteride 
showed that the hydrogen (or deuterium) atoms 
occupied the tetrahedral positions in the face- 
centered cubic titanium lattice to form a 
fluorite structure. The lattice parameter of 
TiDi.971 was 4.440 A. As the composition 
approached TiDe, the cubic unit cell deformed 
to a face-centered tetragonal cell. A further 
study of this deformation was carried out by 
Yakel 8) who demonstrated the occurrence of a 
second order transformation below 310° K. At 
315° K both TiDi,9g and TiHj.99 are cubic with 
lattice parameters of 4.440 A and 4.454 A, 
respectively. As the temperature is decreased 
below 310° K, there is a continuous expansion 
of the crystal in two directions (a and b axes) 
and a contraction in the third (c¢ axis) to form a 
face-centered tetragonal structure. The mini- 
mum c/a ratio reached at 79° K was 0.945. 

A high temperature high pressure X-ray study 
of titanium hydride (at hydrogen contents close 
to TiH2) by Goon and Malgiolio %) revealed a 
maximum in the lattice parameter of the f.c.c. 
structure with temperature at constant hydrogen 
pressure. The temperature at which the maxi- 
mum occurred varied with hydrogen pressure 
(e.g., at 100 psi, 280 psi, and 435 psi, the 
temperatures at which the maximum occurred 
were about 350°, 450° and 500° C, respectively). 
The initial increase in lattice parameter appears 
to be normal thermal expansion. Goon proposes 
that the subsequent contraction is due to the 
filling of octahedral sites in the f.c.c. lattice 
with a corresponding formation of vacant 
tetrahedral sites. Wolfarth 1°) has explained 
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maxima in magnetic susceptibility isotherms 
obtained by Trzebiatowski and Stalinski™) as 
due to antiferromagnetism. However, Wilkinson 
and Shull 12) at Oak Ridge National Laboratory 
found no evidence of antiferromagnetic ordering 
in a neutron diffraction study of titanium 
hydride. The maxima in magnetic susceptibility 
isotherms may be associated with the maxima 
in lattice parameter observed by Goon. There 
is an obvious necessity for further work to 
explain both the tetragonality at lower temper- 
atures and the maxima in lattice parameter and 
magnetic susceptibility at higher temperatures. 


Zirconium Hydride 


At the present time, the hydride of most 
interest from the nuclear reactor standpoint is 
zirconium hydride because of its relatively 
high value of Ny, high temperature stability, 
and the low capture cross section (0.18 barns) 
of zirconium. For these reasons, the zirconium 
hydrogen system has been studied more than 
any other metal-hydrogen system in recent 
years. 


Phase Relations: A composite phase 
diagram for this system is shown in fig. 2. Al] 
eRe) T T T T 1 1 T 
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——-—Vaughan and Bridge 
Fig. 2. Zirconium-hydrogen phase diagram. 


points on the diagram represent boundaries 
between phase fields as determined by different 
methods. The diagram is similar to the titanium- 
hydrogen phase diagram in form, with a lower- 
ing of the « to # transition temperature by the 
presence of hydrogen. The « and f phases are 
the low and high temperature modifications of 
zirconium metal, respectively, containing dis- 
solved hydrogen. The hydride phases, which are 
actually compounds of zirconium and hydrogen 
are designated as the 6 and ¢ phases. Solubility 
studies of hydrogen in «-zirconium by Gul- 
bransen and Andrew!’), and Mallett and 
Albrecht 14) led to the boundary between the 
« and «+6 phases shown in fig. 2. Ells and 
McQuillan 4) measured equilibrium pressures 
from 500° to 950° C for hydrogen contents up to 
ZrHo.s, thus obtaining the limits of the two 
phase («+ ) region. A high temperature X-ray 
phase study of this system by Vaughan and 
Bridge 16) gave very poor agreement for the 
boundary between the two phase («+ /) region 
and the one phase 6 region as shown by the 
dashed curve in fig. 2. The large discrepancy 
between the two investigations is probably due 
to the difference in experimental methods. The 
equilibrium pressure study is inherently more 
precise since this method determines the actual 
points on the phase boundary from breaks in 
the isotherms. On the other hand, the X-ray 
method depends on an estimation of the bounda- 
ry between two regions containing points 
representing the phases present. In addition, 
a phase present in low concentration cannot be 
detected by X-ray diffraction analysis. 

The points representing the boundary between 
the 6 and two phase 6+6 regions have been 
obtained from pressure-composition-tempera- 
ture measurements by Edwards, Levesque and 
Cubicciotti!’), and Hall, Martin and Rees 18), and 
also from a heat content study by Douglas 19) 
as well as from the work of Ells and McQuillan, 
and Vaughan and Bridge. The results of all five 
investigations give very good agreement with 
each other for this boundary. However, there 
was considerable disagreement between three 
different determinations of the boundary 
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between the 6+ 6 region and the hydride region 
as shown in fig. 2. It is doubtful that any of 
the three sets of data shown are entirely correct. 
With respect to the work of Vaughan and 
Bridge, the inherent inaccuracy of X-ray phase 
boundary determinations has already been 
discussed above. The zirconium sample used by 
Hall, Martin and Rees contained a_ large 
amount of oxygen impurity (0.4 wt %), which 
should have a considerable effect on the uptake 
of hydrogen. In addition, the negative slope of 
this boundary is opposite to what one would 
expect in such a phase diagram. The isotherms 
of Edwards et al. may be in error in the high 
hydrogen content region as indicated by the 
rapid decrease in slope of their log P versus 
1/7 isochores 2°) at H/Zr values above 1.6. It 
has been pointed out by Gilbert 24) that this is 
thermodynamically improbable. A current in- 
vestigation 2?) of this system indicates that the 
(6+ 0)/hydride boundary has a positive slope 
and occurs at lower hydrogen contents than 
the boundary given by Edwards et al. 

The X-ray data of Vaughan and Bridge 
indicate that the («+ 6)/6 boundary is vertical 
below 450° C. Further evidence for the position 
of this boundary is given by the heat content 
study of Douglas. However, since it is sometimes 
difficult to reach equilibrium at lower tempera- 
tures, there may be some question as to the 
accuracy of this boundary. 

The two phase 6+ .¢ region has been observed 
at room temperature by several investigators. 
An X-ray study by Gulbransen and Andrew 23) 
indicated that this region extended from 
H/Zr=1.55 to about 1.62. Trzebiatowski and 
Stalinski’s X-ray work 24) gave an approximate 
range for the 6+<« region of H/Zr=1.63 to 1.73. 
The range of this two phase region can be also 
determined from thermoelectric power measure- 
ments on zirconium hydride carried out by 
Bickel 25) at this laboratory. From these data, 
the two-phase region was found to exist between 
H/Zr=1.57 and 1.73. An average of these three 
investigations indicates that the 6+ ¢ region lies 
approximately between H/Zr=1.6 to 1.7 as 
shown in fig. 2. There is no evidence for the 


Or 


existence of this two phase region at higher 
temperatures. The pressure-composition iso- 
therms of Edwards et al.1’), and a more recent 
P-C-T study 22) indicate no plateaus or breaks 
in this composition range as would be required 
by a two phase region. Therefore, the 6+ ¢« two 
phase region must become narrower with 
increasing temperature until it disappears 
entirely at higher temperatures. Whitwham 
et al.26) believe that the ¢ phase is formed by a 
pseudo-martensitic transformation from the 
6 phase. 

There is a strong possibility that the two 
phase 6+e region in the Zr-H system may 
actually not exist. Edwards and Levesque 2’), 
in a study of the Zr-H-O system, observed that 
the two phase 6+. region, present below about 
15 at % oxygen, narrowed with decreasing 
oxygen content. They believed that this two 
phase region is still present in the binary Zr-H 
system. However, since their conclusion was 
based on extrapolated data, the possibility that 
the two phase 6+e region in the binary Zr-H 
system may be due to oxygen impurity should 
be considered. 


Crystal Structures: There has been some 
confusion as to the structures of the phases 
present at high hydrogen compositions in the 
zirconium hydrogen system. Higg 6) reported 
the five phases shown in table 2. The f and y 


TABLE 2 


Zirconium-hydrogen phases as reported by Hagg °) 


Composition : ; Lattice 
(H/Zr) Designation | Structure Paraanions 
0-0.05 Oo Cp.) |= 3-2 0306 

0.25 B £.C.C. a=4.66 

0.50 y h.c.p. | @=3.55; c=5.45 

1.0 0) ABYC: a=4.78 
(fluorite) 

2 € f.c.t. a=4,96 c=4.44 


phases reported by Hiagg have not been 
observed by subsequent investigators. They are 
probably non-equilibrium phases; although the 
f.c.c. 6 phase may really be 6 phase hydride, 
with the reported composition ZrHo.25 being in 
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error. The 8 phase shown in fig. 2 is actually the 
high temperature b.c.c. form of Zr, with the 
lattice parameter increasing 16) from 3.51 to 
3.68 A at ZrHo.s2. 

Gulbransen and Andrew 23) have confirmed 
the 6 and e phases of Hage. In addition, they 
observed a new phase designated as y’, which 
had f.c.t. structure with c/a=1.08 (a=4.61 A, 
c=4.975 A), and co-existed with the « and 6 
phases. From the work of Vaughan and Bridge, 
and the P-C-7 data at high composition, it 
appears that there is a single hydride phase at 
high temperatures; a face-centered tetragonal 
structure with continuously varying c/a ratio. 
The hydride at the (f+hydride)/hydride 
boundary (i.e., the hydride co-existing with 
B-Zr) has a c/a ratio greater than one. As the 
hydrogen content increases, the c/a ratio 
decreases becoming one (f.c.c. structure) at 
about H/Zr=1.5. Further addition of hydrogen 
causes the c/a ratio to become less than one. 
Gulbransen’s y’ phase, therefore, may be the 
high temperature hydride which co-exists with 
6-Zr at high temperatures and is metastable at 
room temperature. 

On the basis of three experimental points 
(at H/Zr=1.22, 1.50, 1.86), Vaughan and 
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Fig. 3. Variation in c/a ratio with hydrogen content 


in e-zirconium hydride. 


Bridge show a linear decrease in c/a ratio with 
hydrogen content over the whole hydride range 
at room temperature. Actually, the range over 
which the c/a ratio changes is much narrower. 
Their last two points completely bypass the 
two phase (6+) region shown in Fig. 2. Also, 
the point at H/Zr=1.22 is in the two phase 
(x+6) region. Trzebiatowski and Stalinski 2+) 
have shown that there is a decrease in c/a ratio 
with hydrogen content at room temperature in 
the single phase ¢ region as shown in fig. 3. This 
decrease was also observed by Sofina and 
co-workers 28), 

An X-ray and neutron diffraction study ?°) 
of ZrHz confirmed a f.c.t. structure for the 
zirconium atoms with the hydrogen atoms 
located in the tetrahedral positions to give 
a slightly tetragonal “fluorite” structure. The 
lattice parameters were a=4.978 A and c= 
= 4.449 A. Yakel8) and Goon 9) agreed with 
these values and also observed an increasing 
c/a ratio with temperature similar to that in 
titanium hydride. 


Thermodynamic Properties: From the 
dissociation pressure data in the two phase 
o+6 region of Gulbransen and Andrew 2°), the 
heat of formation of 6-zirconium hydride from 
hydrogen saturated «-Zr is calculated to be 
—48.3 kcal per mole of Hz absorbed in the 
temperature range 400°-550° C. According to 
fig. 2, the composition of the 6 hydride phase 
formed at temperatures below 450° C is ZrH1.38. 
At 400° C, the reaction can be written: 


Zr Ho .02(x) + 0.68 Hy > Zr H.38(6). 


Using the value of — 27.6 kcal/mole He for the 
heat of solution of hydrogen in «-Zr given by 
Douglas 19), and re-calculating per mole of 
hydride, the heat of formation, AH°aoo°c, of 6 
zirconium hydride is —33.1 kcal/mole hydride 
(ZrH1.3s) at 400° C. From heat content measure- 
ments, Douglas gives a value of — 42.7 kcal/mole 
He for the conversion of hydrogen saturated 
«-Zr to 6 zirconium hydride at 550° C. This is 
5.6 kcal less than the value derived from 
Gulbransen’s 23) data. 

Assuming the composition of ZrH1.5, and an 
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approximate value of 2.25 cal/degree C for the 
heat capacity of the hydride, Gulbransen and 
Andrew 28) give the following relation for the 
free energy of formation of the hydride: 


AF° = — 34 930— 2.25 T In T'+ 43.88 7. 


For the heat of formation of 6-hydride from 
hydrogen saturated f-zirconium, Edwards ?°) 
obtains a value of —49.6 kcal/mole Hz from 
dissociation pressure measurements in the 
range 600°-900° C. Douglas’ value is —50.7 
keal/mole Hz at 550° C. 

For the two phase $+ 6 pressure plateau, the 
data of Edwards for dissociation pressure give 
values of +10850 and 12.22 for A and B, 
respectively. Hall, Martin and Rees 18) measure- 
ments yield values °°) of +10 300 and 11.75. 
For the two phase «+6 region, the values of 
A and B are +10 540 and 11.80, respectively, 
as calculated from Gulbransen’s data. 


Hafnium Hydride 


The phases present in the hafnium-hydrogen 
system at room temperature were determined 
in an X-ray study by Sidhu and McGuire 21). 
The first hydride phase appeared at a very low 
concentration of hydrogen. It was a slightly 
deformed f.c.c. hydride which co-existed with 
x-hafnium metal up to a hydrogen content ratio 
of about H/Hf=1.6. Between H/Hf 1.6 and 1.8, 
the cubic hydride was no longer deformed and 
existed as a single phase. At about an H/Hf 
ratio of 1.8, a face centered tetragonal hydride 
was formed which co-existed with the cubic 
phase up to an H/Hf ratio of about 1.86. From 
an H/Hf ratio of 1.86 to 2, f.c.t. hydride exists as 
a single phase, the c/a ratio decreasing 32) with 
increasing hydrogen content as in e zirconium 
hydride. The lattice parameters of the f.c. 
tetragonal Hf He are 3) a9 = 4.919 and co = 4.363, 
and that of the f.c. cubic phase a9=4.708 A. 
A neutron diffraction study’) of the corre- 
sponding hafnium deuterides, revealed that the 
deuterium (or hydrogen) existed in the tetra- 
hedral positions of a fluorite-type lattice similar 
to the structure of zirconium hydride. 
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Fig. 4. Hafnium-hydrogen phase diagram. 

A pressure-composition-temperature study of 
the hafnium-hydrogen system has been carried 
out by Edwards and Veleckis 34). The phase 
diagram derived therefrom is shown in fig. 4 
along with the phases obtained from X-ray 
studies 3!) (represented by the dotted line). The 
two phase region («-Hf metal-+f.c.c. hydride) 
shows up clearly as pressure plateaus in the 
pressure-composition isotherms, while the two- 
phase hydride region (f.c.c.+f.c.t.) appears as 
a break in the increasing pressure high compo- 
sition portions of the isotherms. The phase 
boundaries obtained from the P-C-T studies 
agree fairly well with those from the X-ray 
study. It has been observed by Sidhu e¢é al.’) 
that the position of the («-Hf + f.c.c.)/f.c.c. 
phase boundary is quite sensitive to the con- 
ditions under which the samples were prepared. 
This may explain the slight discrepancy in the 
position of this boundary as determined by 
the two different investigations. 

From the P-C-7' studies the dissociation 
pressure of the cubic hydride can be expressed 
by the usual relation if 4 =7250 and B=8.65. 
From these data the heat of formation of the 
nonstoichiometric cubic hydride from hydrogen 
saturated «-Hf is calculated as — 33.2 kcal/mole 
hydrogen. 
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The effect of pile-irradiation on the hafnium 
hydrides was investigated by Sidhu, Campos 
and Zauberis 5), Samples of both the cubic and 
tetragonal hydrides were irradiated for almost 
a month in a heavy water reactor under a 
thermal neutron flux of 4.5 <10!2, an epicad- 
mium neutron flux of 3 x 104 and a gamma-ray 
flux of 1600 r/sec. The samples were studied by 
X-rays about a year after removal from the 
reactor. No change in structure of hydrogen 
composition was observed indicating that the 
hydrides undergo no permanent changes under 
reactor radiation. 


2.2. RARE EARTH HYDRIDES 


Except for ytterbium and europium, the 
rare earth metals appear to form two hydrides 
corresponding to the stoichiometric formulas 
MHz: and MHz. Pressure-composition-tempera- 
ture relations of the hydrides of the first four 
metals in the rare earth series were studied by 
Mulford and Holley 3*), and an X-ray and 
neutron diffraction investigation of the same 
four hydrides plus samarium hydride was made 
by Holley et al.°’). Hydrogen has a limited 
solubility in these rare earth metals (La, Ce, 
Pr and Nd). Definite dihydrides are formed 
which co-exist with the corresponding hydrogen 
saturated metal. The dihydrides actually have 
hydrogen contents slightly less than the stoi- 
chiometric value of two. They have a fluorite 
structure, with the lattice parameters shown in 
table 3. Further addition of hydrogen yields a 
continuous solid solution which approaches the 
stoichiometric formula HM3. The neutron dif- 
fraction measurements indicate that the ad- 
ditional hydrogen enters the octahedral sites 
of the fluorite structure, with contraction of 


TABLE 3 


Lattice constants of rare earth dihydrides 


Hydride | Lattice parameter (A) 


LaHe 5.667 + 0.001 
CeHe 5.581 + 0.001 


PrHe | 5.517 -+ 0.001 
NdH2 | 5.470 + 0.001 
SmH2 5.376 + 0.003 


the lattice. Korst and Warf 8) observed that 
the lattices of these hydrides began to contract 
at about MHj.s, indicating the occupation of 
octahedral sites before all the tetrahedral sites 
were filled. 

From the pressure-composition-temperature 
data of Mulford, the dissociation pressures of 
the hydrides and their heats of formation from 
the hydrogen saturated metal, 4p, can be 
calculated in the temperature range 500°—800° C, 
These are shown in table 4 along with the 
corresponding values obtained from P-C-7' 
measurements by Korst. The values given for 
A and B represent the terms in the dissociation 
pressure equation given above. Except for the 
case of cerium hydride, the agreement between 
the two different P-C-7' investigations is 
reasonably good. The value given by Korst is 
probably the more accurate one, since there is 
no apparent reason for the Ap of cerium 
hydride to differ so much from the other three 
hydrides. 

Although the stoichiometric composition MH3 
was not attained by either Mulford or Korst, 
Dialer and Rothe 39) claimed the preparation 
of a cerium hydride corresponding to the exact 
stoichiometric composition CeH3.000+0.002. Such 
accuracy is probably not warranted 38) in view 
of the experimental procedures used. More 


TABLE 4 


Dissociation pressures and heats of formation of rare earth hydrides 


Hydride A B Ay (keal/mole He) 
Mulford Korst Mulford Korst Mulford Korst 

LaHe 10 860 10 850 10.76 10.64 —49.7 —49.6 
CeHe 7 420 10 760 Metts! 10.63 —33.9 —49.2 
PrHe 10 450 10 870 10.23 10.53 —47.8 —49.7 
NdHe 9 800 11 030 9.37 10.48 —44.8 —650.5 


a 
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recently, Mikheeva and Kost) report a 
cerium hydride with a hydride content up to 
CeH3 16. 

A study*) of the gadolinium-hydrogen 
system revealed the formation of two hydrides 
approaching the stoichiometric formulas MHz 
and MHs3. The dihydride has a fluorite structure 
(similar to the hydrides discussed above) with 
a lattice parameter of 5.303 A. The trihydride, 
however, has a structure different from the 
tribydrides of the rare earths previously dis- 
cussed. Rather than a cubic structure derived 
from GdH2 by filling octahedral sites, the tri- 
hydride of gadolinium has a hexagonal structure 
with lattice parameters dp = 3.73 and co= 6.71 A 

From P-C-T studies, the heat of formation of 
GdHz from hydrogen saturated gadolinium 
metal is —46.9 kcal/mole in the temperature 
range 600-800° C, and the dissociation pressure 
of the dihydride is given if A=10 250 and 
B=9.72. The trihydride has a much lower 
thermal stability than the dihydride, and 
probably dissociates (Hz pressure reaches | atm) 
at about 400° C. 

The di- and tri-hydrides of erbium have also 
been prepared 42), they have fluorite and 
hexagonal structures, respectively. The dissoci- 
ation pressure of the dihydride, ErHe, is given 
if A=11780 and B=10.91 and the heat of 
formation is about —54 kcal/mole. 

The hydrides of ytterbium and europium 
differ considerably from the hydrides of the 
other rare earth metals. An X-ray diffraction 
study #*) showed that these metals formed 
hydrides approaching the stoichiometric formula 
MHz, having orthorhombic structures iso- 
morphous with the alkaline earth metal hydrides. 
The lattice parameters of the hydrides were 
a=5.871, b=3.561, c=6.763 A for ytterbium 
deuteride, and a=6.21, 6=3.77, c=7.16 -for 
europium deuteride. The reason for the different 
properties of these two hydrides lies in the 
electronic configuration of the metals. In 
ytterbium, the 4f shell is completely filled thus 
stabilizing these electrons and permitting only 
the two 6s electrons to engage in bonding. 
Similarly, in europium all seven 4f orbitals 


contain one electron (according to Hund’s rule 
of maximum multiplicity) resulting in a stabi- 
lized electron configuration. These two elements, 
therefore, form stable divalent compounds since 
their electron configuration is similar to the 
alkaline earth elements. 

Electrical resistivity measurements 44-46) of 
some rare earth hydrides indicate that they may 
be semiconductors since their resistivities fall 
in the range 0.1-10? ohm cm, and they exhibit 
a negative temperature coefficient of resistivity. 


2.3. HyDRIDES OF THE ACTINIDE SERIES 


ELEMENTS 


The actinide series hydrides do not resemble 
each other as do the hydrides of the two series 
discussed above. The three metals which are of 
most interest to nuclear reactor technology, U, 
Th and Pu, form hydrides of different hydrogen 
to metal ratios and structures. The hydrides 
will be discussed in order of atomic number. 


Actimum Hydride 


An actinium hydride having the stoichio- 
metric formula AcH2 was prepared by Farr 
et al.4”). X-ray structure studies indicated that 
the actinium atoms were arranged in a face 
centered cubic structure with a lattice para- 
meter of 5.670 A. Presumably, the hydrogen 
atoms are situated in the tetrahedral positions 
to form a fluorite structure. 


Thorium Hydride 


A pressure-composition-temperature study of 
the thorium-hydrogen system by Nottorf 48) 
revealed that thorium reacted readily with 
hydrogen to form two definite hydrides, one 
approaching a stoichiometric composition ThHe 
and the other approaching the H/Th ratio of 
3.75. The dissociation pressure of the lower 
hydride corresponds to A=7700 and B=9.54; 
that of the higher hydride to A=4220 and 
B=9.50. The heats of formation of hydrogen 
deficient ThHs from hydrogen saturated thorium 
metal was about — 35.2 kcal/mole He; and the 
heat of formation of hydrogen deficient ThHs3.75 
from non-stoichiometric ThHz was —19.3 
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kcal/mole Hz. A subsequent P-C-7' study of the 
dihydride by Mallett and Campbell 4%) gave 
A=7500 and B=9.35 for the dissociation 
pressure of the nonstoichiometric di-hydride, 
and — 34.3 kcal/mole He for its heat of formation. 
These values are in good agreement with those 
of Nottorf. It should be pointed out, however, 
that in both investigations the pressure of the 
two-phase regions (Th+ThHe) and (ThH2+ 
ThHs3.75) were not constant as required by the 
phase rule, but increased slightly with increasing 
hydrogen content of the system, thus leading 
to sloping plateau regions of the isotherms. Since 
this is due to either impurity in the samples or 
thermal gradients 5°) along the sample during 
measurements, the values given may be some- 
what in error. 

From X-ray and neutron diffraction studies ?9), 
the crystal structure of ThHs was found to be 
isomorphous with zirconium hydride having a 
face centered tetragonal structure with a= 
5.80 A and ¢)=5.03 A (or b.c.t. with ap = 4.10 A). 
The hydrogen atoms are situated in the tetra- 
hedral positions. The structure of the higher 
hydride was determined by Zachariasen 5!) 
from X-rays to be body centered cubic with a 
lattice parameter of 9.11-+ 0.02 A. The unit 
cell contains 16 thorium atoms and 60 hydrogen 
atoms which definitely establishes the chemical 
formula Th4Hj5 for the stoichiometric hydride. 


Protactinuum Hydride 


A protactinium hydride was prepared by 
Sellers e¢ al.52) by the reaction of protactinium 
metal with hydrogen at 250°. The hydride is 
cubic and isostructural with uranium hydride 
(see below) having a lattice parameter of 
6.648 A. 


Uranium Hydride 


The first comprehensive investigation of the 
uranium-hydrogen system was carried out by 
Spedding and _ co-workers 53), Uranium was 
found to react readily with hydrogen to form a 
single pyrophoric hydride with the stoichio- 
metric formula, UH 3. From pressure-compo- 
sition-temperature studies the dissociation pres- 
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Fig. 5. Pressure-composition isotherms for the ura- 
nium-hydrogen system at 357° C (Spedding e# al.®3)). 


sure corresponds to 4=4500 and B=9.28 and 
the heat of formation was calculated to be 
— 30.8 kcal/mole hydride. In studying the 
357° C pressure-composition isotherm of this 
system, Spedding observed that the equilibrium 
pressures on hydriding were higher than those 
obtained on dehydriding. This hysteresis effect 
has also been observed in other metal-hydrogen 
systems 1). As shown in fig. 5, the dehydriding 
isotherm also had a dip in the high composition 
end of the isotherm which could be reproduced 
in subsequent experimental runs. No expla- 
nation for these effects are given. 

A dissociation pressure study of UHs3 by 
Flotow and Abraham *4), in the temperature 
range 300°-420° C, yielded the following dis- 
sociation pressure equation: 


6270.1 
fi 


log 1oP (mm) =— — 6.0146 log 107’ + 28.9235, 


and a heat of formation equal to — 30.5 keal/ 
mole hydride. A high pressure dissociation 
study 55) (up to 60 atm of Hz) gave the usual 
simpler expression for the dissociation pressure 
in the range 450°-650° C with A=4410 and 
B=9.14 from which a heat of formation of 
— 30.3 kcal/mole hydride was obtained. This 
value is in good agreement with the values 
obtained by Spedding and Flotow and Abraham, 
and also with a calorimetrically dtermined 5) 
value of —30.35 kcal/mole hydride at 25°C. 
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Therefore, it can be seen that the heat of forma- 
tion of UH3 remains reasonably constant over 
the wide temperature range 0—650° C. 

From an X-ray diffraction study’), the 
uranium atoms in the UHs structure were found 
to form the cubic f-tungsten structure with 
8 U atoms per unit cell at (0, 0, 0), (4, 4, 4), 
Gs 2 0), (0, L B)» (25 0, t)s (a; B> 0), (0, i 3), 
and (4, 0, ?). The lattice parameter of uranium 
hydride is 6.645 A. The positions of the hydrogen 
atoms in the UHs lattice were determined by 
neutron diffraction 58). They are situated in 
the center of deformed tetrahedra equidistant 
from four U neighbors, the U-H distance being 
2.32 A. Each uranium atom is surrounded by 
twelve hydrogen atoms. 

When uranium hydride is prepared at low 
temperatures, a new phase (x) appears ®9), the 
lower the temperature the greater the pro- 
portion of «-phase in the /-phase (f-phase being 
the phase described above). The «-phase is 
stable at 100°C but converts to f-phase at 
250° C. This low temperature phase can also be 
prepared by an electrolytic method ®). The 
structure of this phase is body centered cubic 
with two uranium atoms per unit cell at 
(0, 0, 0) and (4, 4, 4) and six hydrogen atoms 
at + (}, 0, 4), G, £0) and + (0, 3, 3), and 
a lattice parameter of 4.155-+ 0.05 A. 

Magnetic susceptibility measurements by 
several investigators 61-63) showed that uranium 
hydride becomes ferromagnetic at about 173° K. 
Some rough electrical conductivity measure- 
ments *3) on uranium hydride indicated that its 
resistivity was of the same order of magnitude 
as uranium metal. 


Neptunium Hydride 


Neptunium metal is reported 4) to react 
rapidly with hydrogen at 50° C to form a hydride 
having the composition NpHs3.6-3.3. If this 
hydrogen to metal ratio is correct, it may be 
similar to the higher thorium hydride. 


Plutonium Hydrides 


Plutonium metal forms two hydrides having 


the stoichiometric formulas PuHs, and PuHs. 
Pressure-composition-temperature studies ®) 
yielded A=8165 and B=10.01 for the dis- 
sociation pressure of the dihydride, and a 
corresponding heat of formation from hydrogen 
saturated plutonium metal of —37.4 kcal/mole 
He. The dihydride has a fluorite structure 
similar to the rare earth hydrides with a lattice 
parameter of 5.359 A. 

Mulford and Sturdy ®) observed that PuHe 
will dissolve hydrogen in the octahedral inter- 
stices of its lattice with a corresponding decrease 
in lattice parameter. At about PuHe.7 (the 
actual H/Pu ratio varies with temperature) a 
hexagonal phase is formed which seems to co- 
exist with the cubic PuHz up to a hydrogen 
content of about PuHe.9 at which composition 
only the hexagonal phase is present. This 
behavior is analogous to that of gadolinium 
hydride except that the boundaries of the two 
phase (cubic hydride + hexagonal hydride) region 
are not as sharply defined. The lattice para- 
meters of the hexagonal tri-hydride are ad9= 3.78 
and co=6.76 A. An X-ray study of PuHe.7 by 
McDonald and Fardon §’) gave values for the 
lattice parameter of the hexagonal hydride 
which were identical with those of Mulford. 
Their value for the cubic hydride was 5.34 A 
which was the value given by Mulford at the 
composition of PuHg.5. 


Americium Hydride 


Americium metal reacts ®8) violently with 
hydrogen at 50°C to yield a hydride of the 
formula AmH2,7+09.3. A lower hydride, AmHe 
has also been reported ®9). 


2.4. HYDRIDES OF THE GROUP VA METALS 


In general, the Group VA metals form 
hydrides approaching the stoichiometric formula 
MH. The change from metal to hydride seems 
to be a continuous one with no sharp change in 
structure. T'wo phase regions, if they exist, seem 
to occur at much lower temperatures than they 
do for hydrides of the Groups III A and IV A 


metals. 
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Vanadium Hydride 


The hydride having the highest hydrogen to 
vanadium ratio was VHo..94 prepared by 
Trzeciak et al.7°) from iodide process vanadium. 
X-ray diffraction showed the hydride to be body 
centered tetragonal with a—3.02 A and 
¢o=3.36 A. Vanadium metal has a _b.c.c. 
structure 7) with aj=3.028 A. A neutron 
diffraction study 7) of VDo.75 gave a b.c.c. 
structure for the hydride with ap = 3.148 A. The 
positions of the deuterium atoms could not be 
determined above 207° K. Below this temper- 
ature, the deuterium atoms ordered themselves 
into a primitive cubic lattice with a= 6.30 oe 
twice the lattice parameter of the vanadium cell 
in the hydride. Therefore, it appears that the 
b.c.c. vanadium lattice expands with addition 
of hydrogen, eventually distorting to a tetra- 
gonal lattice at high hydrogen contents. 


Niobium Hydride 


Because of its high temperature properties 
and relatively low cross section, niobium has 
recently become of interest to the nuclear 
reactor field. Pressure-composition-temperature 
and X-ray studies on the niobium-hydrogen 
system have been carried out by Albrecht 
et al.”8). The P-C-T studies show that there is 
continuous solubility of hydrogen in niobium 
above 300° C. Below this temperature, however, 
a hydride phase is formed. Albrecht showed that 
this hydride or 6-phase is body centered cubic 
with ag= 3.42 A. The x-phase is also b.c.c. with 
a lattice parameter 7!) of 3.301 A. At room 
temperature the two phases, hydrogen saturated 
a-Nb (a9=3.308 A) and £-niobium hydride, 
co-exist between the approximate compositions 
0.09 to 0.7 H/Nb. Above H/Nb=0.7, only f- 
phase is present, its lattice parameter increas- 
ing 4) to 3.447 A at NbHo.s3. 

According to Brauer and Hermann”), the 
B-phase is actually a deformed b.c.c. structure 
and can better be considered as face centered 
orthorhombic with aj=4.84, bo=4.90, and 
co=3.45 A at NbHo.59. 

The heat of solution of hydrogen in niobium 
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increases from 16.1 kcal/mole at H/Nb=0.01 
to 23.3 keal/mole at H/Nb=0.70. 


Tantalum Hydride 


Horn and Zeigler 7) found hydrogen to be 
continuously soluble in tantalum metal (b.c.c.”), 
do = 3.303) up to H/Ta=0.49 with the lattice 
parameter increasing to 3.399 A. Brauer and 
Hermann ”) report a solubility of H/Ta=0.2 
with the lattice parameter increasing to 3.34 A. 
The f-hydride forms at this point. The solubility 
reported by Stalinski 7’) was H/Ta=0.37, the 
lattice parameter of the metal increasing to 
3.363 A. The reason for this large variation in 
reported solubilities can be partially explained 
from fig. 6 which is a partial phase diagram of 
the Ta-H system determined from heat capacity 
data, X-ray measurements, and electrical resist- 
ivity measurements by Waite, Wallace and 
Craig 78). The slope of the boundary between 
x-tantalum and the «+f two phase region is 
fairly shallow at about room temperature thus 
causing the solubility limit of hydrogen in 
tantalum to be sensitive to slight temperature 
variations. 

Tantalum metal is body centered cubic, the 
lattice parameter of the a-phase expanding 
with increasing hydrogen content according to 
the relation: 
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Fig. 6. Partial phase diagram of the tantalum- 
hydrogen system. 
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ao = 3.303 + 2.5 x 10-3 N, 


where N is the atomic percent of hydrogen. The 
B-phase is body centered tetragonal at about 
Manes with do=3.98 A and c=3.41 A. The 
hydrogen atoms are believed to occupy the 
octahedral interstices (4, 4, 0). With increasing 
hydrogen content, the (f-phase gradually 
becomes orthorhombic. At TaHo.s, the lattice 
parameters %) of the orthorhombic structure 
are do = 4.728, bo =4:778 and co=—3.428 A. 

Nuclear magnetic resonance studies ) of 
tantalum hydride indicated that the hydrogen 
atoms had.strong internal motion above 215° K 
but become tightly bound below this tempera- 
ture. 


2.5. HYDRIDES OF METALS ABOVE GROUP VA 


The transition metals above Group V do not, 
in general, form hydrides. Some of them such as 
Mn, Fe, Ni will dissolve small quantities of 
hydrogen. However, hydrides of Pd, Cr and Cu 
have been prepared. 


Palladium Hydride 


Historically more study has been devoted to 
the palladium-hydrogen system than any other 
metal-hydrogen system. Since an_ excellent 
review of this system is given by Smith 4), 
only a short discussion will be given here to 
include some of the more recent work. 

Hydrogen dissolves in face centered cubic 
palladium with an expansion 7) of the lattice 
from 3.891 A to 3.894 A at PdHo.95. At about 
this composition, a new phase is formed having 
the same f.c.c. structure as palladium but a 
larger lattice parameter, 4.027 A. This new 
phase can be considered a hydride phase but is 
usually designated as f. The two phases co-exist 
until an H/Pd ratio of about 0.7, above which 
only the hydride phase exists. The (-phase 
continues to expand with addition of hydrogen 
to about 4.07 A. Pressure-composition-temper- 
ature studies 89:81) indicate that the stoichio- 
metric composition would be PdH. The dis- 
sociation pressure of the hydride is given by the 
usual relation if d=1877 and B=7.48. From 


a calorimetric study of the heat of absorption 
and desorption of hydrogen in palladium, Nace 
and Aston 82) have calculated the heat of 
formation of palladium hydride (of hydrogen 
content PdHo5) from pure palladium as 
—9.440 keal/mole Hz. Nace and Aston 8%) have 
also made a low temperature thermodynamic 
study of this system. 

From a neutron diffraction study, Worsham 
et al.84) have established that the hydrogen 
atoms in palladium hydride are situated in the 
octahedral interstices of the f.c.c. palladium 
lattice to give a hydrogen deficient sodium 
chloride structure. 


Chromium Hydride 


By electrodeposition of chromium at high 
current densities, Snavely ®) has claimed the 
formation of two chromium hydrides; a hexa- 
gonal hydride with limiting stoichiometric 
formula CrH, and a face centered cubic hydride 
whose limiting stoichiometric formula is CrHs. 
The maximum H/Cr ratio observed was 1.7. 
Trzeciak et al.”°) prepared a chromium hydride 
of maximum H/Cr equal to 1.6. It decomposes 
above about — 10°C. 


Copper Hydride 


The direct combination of copper and 
hydrogen gas does not yield a hydride. However, 
a hydride, CuH, can be prepared by the reaction 
of copper sulfate with hypophosphorous acid 86). 
The structure is the hexagonal wurtzite structure 
with lattice parameters a = 2.89 and cp = 4.61 Ae 
A neutron diffraction study 8’) showed that the 
hydrogen atoms were at the (0, 0, 2), (4, 3, 3) 
positions and the copper atoms at the (0, 0, 0) 
and (4, %, 4) positions. 


2.6. ALLOY HYDRIDES 


Although there has been a fair amount of 
work on the solubility of hydrogen in alloys, 
until recently virtually no work has been done 
on the formation of hydrides from alloys. 
However, because of the advantages of develop- 
ing new hydrides with desirable properties for 
nuclear power applications, studies are being 
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performed on various alloy-hydrogen systems, 
at the present time. 

Trzeciak et al.7°) have investigated the hydro- 
gen retention of a large number of titanium and 
zirconium alloys. Under hydriding conditions 
which gave a H/Ti ratio of 1.75 for pure titanium, 
the addition of vanadium caused an increase of 
the hydrogen to metal ratio to 1.96 at 14 at % V. 
At 69 at % V, a maximum Ny equal to 10.1 
was reached. The addition of Nb, Mo, Mn and 
Zr also caused some increase in the hydrogen 
retentive capacity of titanium. The maximum 
H/M ratio of 2.16 was obtained with 18 at % 
Mo. Chromium caused a decrease in H/M ratio. 

For the case of zirconium, Nb, V, Mo, Sn 
and Cr all caused a decrease in the hydrogen to 
metal ratio. On the other hand, addition of 
2.3 at % Ni caused an increase of H/M from 
1.92 to 2.17. However, the presence of 33 at % 
Ni caused a decrease to an H/M of 1.04. 

A study of the effect of alloying elements on 
the thermal stability of zirconium hydrides 
indicated no large effect caused by low concen- 
trations of alloying elements. The formation of 
intermetallic compounds, however, seemed to 
result in a decrease in thermal stability. 

Gulbransen et al.88) studied the reaction of 
50 wt percent U-Zr alloy with hydrogen. In all 
cases, the only hydride phase formed seemed 
to be zirconium hydride. 

Libowitz, Hayes and Gibb 89) showed that 
the intermetallic compound ZrNi formed a 
definite hydride approaching the stoichiometric 
composition ZrNiH3. The maximum H/Zr ratio 
observed was 2.93. This hydride has the 
structure of the intermetallic compound, which 
is othorhombic %), rather than one resembling 
ziconium hydride. The dissociation pressure is 
considerably higher than that of zirconium, 
reaching one atmosphere at about 275° C. The 
pressure-composition isotherms exhibited a very 
large hysteresis effect. 

At the temperature studied (100-250° C), 
close to one hydrogen atom per zirconium atom 
dissolved in the alloy before formation of the 
trihydride. It has not been established whether 
a monohydride, ZrNiH, exists at lower temper- 


atures. In any case, the heat of the reaction: 


ZrNiHl + }(2—2)H» > ZrNiHs-z 


was found to be —18.4 kcal/mole He. 

A study of the reaction of the intermetallic 
compound Zr2Ni with hydrogen was carried out 
by Gibb and Krutenat 91), The sample contained 
some excess nickel. Assuming this was present 
as ZrNi as indicated by the phase diagram of the 
Zr-Ni system 92), the maximum composition of 
the hydride formed would be ZrgNiH4.35. This 
is higher than the value Zr2NiH3.12 obtained by 
Trzeciak et al.5’), Further investigations are 
being carried out on the Zr-Ni-H system. 


3. The Saline Hydrides 


Although the hydrides of the alkali and 
alkaline earth metals are not classed as tran- 
sition metal hydrides, a brief discussion of their 
properties will be given here as a basis for 
comparison with transition metal hydrides. 
These hydrides are also of interest to nuclear 
technology mainly for shielding applications. 

The hydrides of the Group IA and ITA 
metals are predominantly ionic and are, there- 
fore, often referred to as saline or salt-like 
hydrides. The hydrides of the alkali metals have 
the face centered cubic sodium chloride structure 
and resemble the alkali halides in their physical 
properties. The lattice parameters are shown 
in table 5. Lithium hydride is the only hydride 
which can be melted at normal hydrogen 
pressures, the others dissociating before reaching 
their melting point. The dissociation pressures 
are related to temperature by the values of 
A and B shown along with heats of formation 
in table 5. 

The alkali hydrides are also similar to the 
alkali halides in the formation of color centres. 
Pretzel 98) has observed several absorption bands 
in lithium hydride one of which is due to the 
F center, its maximum occurring at 5660 A 
thus giving rise to a bluish color in the crystal. 

The alkaline earth hydrides from calcium 
through barium form orthorhombic crystals 
whose lattice parameters are given in table 6. 
These hydrides also dissociate before melting. 
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| Their heats of formation and dissociation 
pressures are also given in table 6. 
Magnesium hydride, MgHe, was prepared 1°) 
_ by heating magnesium in hydrogen under high 
pressures (up to 70 atm). Its dissociation pres- 
sure is quite high (reaches one atm at about 


considering the energy necessary to dissociate 
the hydrogen molecule, and also the high 
exothermic heats of formation of the hydrides, 
it appears that rather strong metal-hydrogen 
bonds are formed in a hydride. Alternatively, 
the nature of the saline hydrides supports the 


TABLE 5 


Properties of alkali metal hydrides 


; a mtien Dissociation pressures | Heats of formation 
Hydride attice parameter %) (A) Al (kcal/mole) hydride 
Lit 4.083 8224 9.926 9) — 21.67 9) 
NaH 4.897 5958 10.47 94) —13.49 %) 
KH 5.708 6300 11.86 97) —13.82 %) 
RbH. 6.049 4930 9.51 97) —11.3 97) 
CsH. 6.389 4410 9.25 97) —10.1 97) 
TABLE 6 
Properties of the alkaline earth hydrides 
: Lattice parameters %) (A) Dissociation pressures |Heat of formation 84) 
Hydride 
ao bo Co A B (keal/mole) 
CaHe 6.851 5.948 3.607 10870 11.49 97) —45.1 
SrHe 7.358 6.377 3.882 10400 11.10 100) —42.3 
BaHe 7.845 6.801 4.175 6450 8.20 100) —40.9 


285° C) and can be expressed by the relation 10): 


3 
= a — 2.441 logioT’ 


+4,524x 10-37 — 2.277 x 10-6 T2 4 12.214. 


logio P (atm) = 


The heat of formation of MgH2 is — 18.6 kcal/ 
i mole at 25°C. Its crystal structure is body 
| centered tetragonal with aj=4.517 A and 
Co= 3.021 A. There are two magnesium atoms 
per unit cell at (0, 0, 0) and (4, 4, 4) and four 
hydrogen atoms at + (x, x, 0),+ (x+4, 4-2, 4) 
where x=0.306. 


4. The Nature of Transition Metal Hydrides 


In reviewing the phase relations of transition 
metal hydrides, it is seen that they can all 
deviate from their stoichiometric compositions 
to a considerable extent. It has generally been 
assumed that hydrogen enters a transition metal 
lattice interstitially as protons with the extra 
electrons being distributed in the d-band (or 
f-band) of the metal-like lattice. However, 


notion that the crystal can be regarded as made 
up of hydride anions and metal cations. 
The particular case of LiH being almost 
literally this structure. The nature of the metal- 
hydrogen binding and the reason for the 
deviations from stoichiometry are probably the 
two problems of most interest in transition 
metal hydrides. 


4.1. THEORETICAL STRUCTURAL CONSIDER- 


ATIONS 


By comparing experimentally determined 
pressure-composition-temperature relations with 
those derived theoretically on the basis of a 
proposed model, insight can be gained into the 
structure of hydrides. Using the assumption 
that the energy of absorption of hydrogen 
increases as the number of sites already occupied 
increases, Lacher 103) derived, from statistical 
mechanical considerations, the pressure-compo- 
sition-temperature relations for the absorption 
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of hydrogen by a metal. For the hydrogen 
pressure, P, as a function of the fraction of 
interstitial sites occupied, 0, and temperature, 
T, he obtained an expression which can be put 
into the form: 

47 
fl 


InP =InPy+2In;", + 420-26) (1) 


where Po is the pressure at 0 =} (i.e., the pressure 
of the two-phase region, metal+ hydride), and 
T'. is the critical temperature or temperature at 
which the metal and the hydride are miscible. 
The agreement between the isotherms calculated 
from this equation and those observed for the 
palladium-hydrogen system were quite good. 
Harasima et al.04) extended Lacher’s derivations 
by considering the effect of the expansion of the 
lattice on the energy of the interstitial protons 
during the absorption of hydrogen. This resulted 
in an additional 62 term in eq. (1) and asym- 
metry of the isotherms with respect to the 
critical point. This asymmetry is actually 
observed experimentally. 

Using a thermodynamic treatment and as- 
suming that the nonstoichiometry of hydrides 
was due to the solid solution of metal in the 
hydride phase, Messer 15) derived the following 
expression for pressure-composition isotherms 
in metal-hydrogen systems: 

InP =In Po+ =n + Se(1-2%) (2) 
s s 


S—n st 


where n is the hydrogen to metal ratio at 
pressure P, and s is the hydrogen to metal ratio 
in the stoichiometric hydride. Since 0=n/s, this 
equation agrees with that of Lacher except for 
the 1/s factor of the last two terms in Messer’s 
equation. 

A somewhat broader statistical mechanical 
approach to the problem of interstitial solid 
solutions was taken by Rees 1%), He assumed 
that all the interstitial sites are not crystallo- 
graphically equivalent as hydrogen atoms are 
added. For sites of the first kind, the absorption 
energy is /,; for the second kind Hs, ete. The 
number of sites of any particular kind available 
for occupation depends the upon number of the 
previous kind already occupied. For example, 
sites of the 12 kind are created by occupation 
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of sites of the (i— 1) kind. Using this assump- 
tion and a method similar to Lacher’s, Rees 


derived component isotherms for each kind of — 


site as follows: 


60; Do — 2H; — 46 Lu 
D's 9 
In F = eine Sg ee +In f(7’) 
bP= Ine 
6; —84 41) 
1 Bus 
+ oF | Po— 2E w41) —4 aie Buscn | 


+Inf (2) 


where Dp is the dissociation energy of molecular 
hydrogen, H; is the energy of a hydrogen atom 
in a site of the 2t4 kind, Hy, represents the 
attractive interaction energy of two interstitial 
H atoms in sites of the ith kind, and f(7’) is a 
slowly varying function of temperature. The 
total pressure-composition isotherm derived 
from these component isotherms gives rise to a 
separate plateau region for each kind of site 
present. Therefore, these equations are particu- 
larly applicable to metal hydrogen systems in 
which more than one hydride phase is formed, 
or in which the f-phase of the metal is formed 
in addition to a hydride phase. The assumption 
is made that the effect of changes in lattice 
symmetry and cell dimensions are negligible. 

Martin and Rees 107) have applied these 
equations to the P-C-7' data on the zirconium- 
hydrogen system 18). A value of 66.6 kcal/mole 
was calculated for the energy of hydrogen atoms 
in f-Zr sites and 70.4 kcal/mole for H atoms in 
the hydride phase. The interaction energies were 
2.52 and 4.57 kcal/mole, respectively. The heat 
of formation of ZrH»z was calculated to be 
—41.4 keal/mole. This is slightly lower than 
the experimental values of Douglas 19) and 
Edwards 1°). 

These equations have also been applied to the 
titanium-hydrogen system by Kant 198). The 
values of H; and H;; were 63.9 and 1.0 kcal/mole, 
respectively for the f-phase, and 66.9 and 
3.5 kcal/mole for the y-phase. A calculated value 
of —32.9 kcal/mole was obtained for the heat 
of formation of TiHs. 

Except in Messer’s thermodynamic treatment, 
these methods assume interstitial solid solution 
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of the hydrogen in a metal lattice. Because of 
their large deviations from stoichiometry, 
metallic appearance, and high electrical conduc- 
tivity, transition metal hydrides have been 
referred to as solid solutions or interstitial 
compounds. This would imply that the structure 
of the metal lattice in the hydride has the same 
structure as in the pure metal. That this is not 
the case for the Group III A and IV A metals 
can be seen from table 7 where only cerium and 
actinium do not change structure on hydriding. 
Therefore, transition metal hydrides appear to be 
definite chemical compounds (i.e., compounds 
having a characteristic structure different from 
that of the component elements). Even for the 
case of cerium and actinium, the large dis- 
continuous increase in lattice parameter indi- 
cates the formation of a definite compound 
rather than a solid solution. 

Nuclear magnetic resonance studies 7%) on 
several hydrides gave broad resonance line 
widths at room temperature for TiHe, ZrHe 
and ThHg, of the same order of magnitude as for 
the saline hydrides CaHz and NaH. The broad 
line width indicated that the hydrogen atoms 
were tightly bound in the hydride lattice. On 
the other hand, the line widths of tantalum, 
vanadium 19), and palladium 11°) hydrides indi- 
cate internal motion of the hydrogen atoms. 

If the group III A and IV A hydrides are 
definite chemical compounds, their non-stoichio- 

TABLE 7 


Structures of metals and corresponding hydrides 


Metal | Metal structure |Hydride| Hydride structure 
fie |b .e=p. TiH2 f.c.C. 
Zit, | b.e.p. ZrHe fe Cate 
Jahr || laegay HfHe HCO be 
La_ | h.e.p. LaHe TOKO, 
Ce | f.c.c. (ao =5.16 A) | CeHe f.c.c. (a9 =5.58 A) 
IPAs laneryey PrHe tiOH OE 
Nd | h.c.p. NdHe | f.c.c. 
Sm _ | rhombohedral SmHe | f.c.c. 
Caeeeiresp. Gdilig® || f.c-c: 
Ac | f.c.c. (a@0a=5.31 A) | AcHe | f.c.c. (ao =5.67 A) 
Ahe a f.c.c, ThHa fhe 
ThaHis b.c.c: 
Paw aibvest. PaHs3 b-W (cubic) 
UW orthorhombic UH3 B-W (cubic) 


metry can be attributed to lattice defects, 
either hydrogen vacancies or metal interstitials, 
since the hydrides are always hydrogen deficient. 
Assuming the presence of lattice . defects, 
Anderson 111) has derived pressure-composition- 
temperature relations for binary one-to-one 
compounds. Libowitz 2) modified Anderson’s 
treatment to include compounds of the form 
ABs, where s is not necessarily equal to unity. 
For the case of anion vacancies in the lattice, 
the following expression was obtained: 


ln Pin Po-e 2 In 


nN 2H pp (/2n 
s—n | kT (a iy: 1) -) 
where 2Hpp/z is the attractive interaction 
energy between two defects, and z is the number 
of nearest neighbor hydrogen sites around each 
hydrogen site. On the other hand, assuming 
metal interstitials, the following expression was 
obtained 118) ; 

2 N+ Nx —S 
InP =In Po + =In{"="*—*] 


s—n 


(4) 


Bi 2 Epp [Pe] 


skT an 

where « is the proportionality constant between 
the number of interstitial positions and the 
number of metal atoms on normal lattice sites. 
x depends on the structure of the hydride. 
Eq. (4) differs from the equation previously 
given 2) for metal interstitials which was 
derived from incorrect assumptions. Eqs. (3) 
and (4) can also be derived from strictly thermo- 
dynamic principles by treating vacancies and 
interstices as chemical entities. 

If the composition in eq. (3) were expressed 
as the fraction of vacant sites rather than 
hydrogen to metal ratio, this equation would 
be essentially identical to the ones given by 
Anderson 11!) and Lacher. The agreement with 
Lacher’s equation is not surprising since Lacher 
was assuming interstitial hydrogen in a metal 
lattice which is equivalent to hydrogen vacancies 
in a hydric lattice. However, eq. (4) differs from 
the previous equations by the 1/s factor in the 
last two terms, even when the composition is 
expressed in fraction of interstitial sites occu- 
pied. The equation does becomes identical to 
Anderson’s for one-to-one binary compounds 
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(i.e., when s=1). The 1/s factor arises from the 
fact that the formation of an interstitial metal 
atom, which is in excess over the stoichiometric 
ratio, is associated with the loss of s hydrogen 
atoms from the lattice. Since Messer 1%) was 
implicitly assuming the formation of interstitial 
metal (by considering the activity of metal 
rather than the activity of vacancies) his equa- 
tion also contained the 1/s factor. Messer’s 
equation differs from eq. (4), however, because 
its derivation did not include interstitial 
positions as a chemical entity. 

It should be pointed out that eqs. (3) and (4) 
were derived assuming only one type of defect 
is predominant, and that the interaction energy 
term, Hpp, is independent of the concentration 
of defects. By comparing with experimental 
P-C-T data, these equations can be used to 
determine whether hydrogen vacancies or metal 
interstitials are the cause of the non-stoichio- 
metry of hydrides, as well as to calculate defect 
interaction and formation energies. 

The application 112) of these equations to the 
uranium-hydrogen system indicated that hydro- 
gen vacancies were present in the uranium 
hydride lattice. The interaction energy between 
two vacancies was calculated as 66. kcal/mole. 
During the dissociation of a hydride, the 
following processes can be considered to occur 
(assuming the formation of hydrogen vacancies) : 


Ep+ di (lattice site) + []+ A (surface) 
H (surface) > 4H2 (gas) +3Do 
—> 3(Nvz/N) (2+ 3(Nvz/N) (2H pp/z) 


(Ny2/N) 


where Hp is the energy necessary to form a 
hydrogen vacancy, [], Do is the dissociation 
energy of hydrogen, (Nyz/N) is the number of 
vacancy pairs per vacancy, Ny is the number of 
hydrogen vacancies in the lattice, and N is the 
number of hydrogen sites. When a lattice is 
dissociated Ny=N. Therefore, the total energy 
required for dissociation of a hydride is 
Ep—43Do—Epp. The dissociation equation for 
a hydride, MHs, is usually written: 


MH, -> M+ 4sHo+AH. 


Writing this equation per gram atom of 
hydrogen: 

1MH,~2M+4H,+@ 
where AH is the heat of formation per mole of 
hydride. Therefore, 

at = —(Ep—4Do— Epp). 
The energy necessary to form a vacancy can 
now be calculated from the heat of formation 
of the hydride, and the interaction energy 
between vacancies, (2Hpp/z). For uranium 
hydride a value of 68.3 kcal/mole was calculated 
for Ep. 

In general, the dissociation of a hydride can 
be described in the following manner. As 
hydrogen is removed from a hydride near its 
stoichiometric composition, defects are formed 
(either hydrogen vacancies or interstitial metal). 
These defects attract each other with an energy 
2Epp/z and agglomerate. As more hydrogen is 
removed, a point is reached at which the lattice 
cannot accommodate any more defects. Further 
removal of hydrogen causes the hydride lattice 
to break down or dissociate into two phases; 
the hydride saturated with defects, and the 
metal saturated with dissolved hydrogen (or a 
lower hydride). The nonstoichiometric homo- 
geneity range is determined by the attractive 
interaction energy between defects, 2 pp/z. The 
greater the value of Hpp, the smaller the concen- 
tration of vacancies necessary to break down 
the lattice. 


4.2, THE IONIC CHARACTER OF TRANSITION 
METAL HYDRIDES 


The nature of the chemical bond in transition 
metal hydrides is not known. However, there 
is evidence that the bond in the Group III A 
and IV A hydrides has some ionic character. 
Attempts to calculate a hydrogen atom radius 
from interatomic distances in these compounds 
on the basis of metallic or covalent bonding 
have led to variable values. Because of the 
similarity between the heats of formation of 
rare earth hydrides and_ saline hydrides, 
Dialer "4) suggested that the metal-hydrogen 
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_ bonds in the rare earth hydrides of composition 


MH are ionic. Libowitz and Gibb 45) extended 
the assumption of ionic bonding to all the 
hydrides of the Group IJI A and IV A metals 
including the lanthanide and actinide series and 
obtained a constant value of 1.29-+ 0.05 A for 
the hydride anion radius using Zachariasen’s 
cation radii and coordination number correc- 
tions 116) as shown in table 8. The only hydride 
which does not give good agreement with this 
picture is HfHz which has a calculated H- 
radius of 1.20 A. Regardless of the valence 
exhibited by the metal in the hydride, the radii 
used are for cations with the rare gas electronic 
structure, except for the lanthanide series where 
the filling of the 4f shell was considered to have 
negligible influence on the xenon electronic 
structure, and plutonium where the more 
stable tetravalent state was used. Therefore, in 
every case, the charge of the cation exceeds its 
valence in the hydride. The extra electrons 
which do not contribute to the formation of 
hydride anions probably take part in residual 
metallic bonding thus giving rise to the high 
electrical conductivity of these compounds. 
This ionic character would also explain the 
brittleness of these compounds. This type of 
bonding also seems to be characteristic of the 
transition metal chalcogenides. Studies 1’) on 
some rare earth monochalcogenides indicated 


TABLE 8 


Hydride anion radii in metallic hydrides 


: M-H C.N. : 
Hydride | distance (A) | correction H-radius, (A) 
TiHe 1.93 0.08 125 
ZrHe 2.09 0.08 1.24 
LaHe 2.45 0.08 e333 
CeHez 2.42 0.08 32 
PrHe 2-39 0.08 esl 
NdHe 2.37 0.08 1.30 
SmHe hve 0.08 1.28 
GdHe 2.30 0.08 We27 
AcHa 2.46 0.08 ea 
ThHe 2.41 0.08 1.34 
PaHs 2.32 0.19 2's 
UHs3 Mays 0.19 We2i 
PuHe eae 0.08 1.34 


Avs 1:29 48005,04 


that trivalent cations and divalent anions were 
present in the lattices, with the additional 
electron per atom contributing to the metallic 
nature of these compounds. Magnetic and 
crystallographic studies of thorium 8) and 
uranium 119) sulfides showed that tetravalent 
cations were present even in the subsulfides of 
these metals, with metallic bonding accounting 
for the additional electrons. 

The value of 1.29 A for the ionic radius of 
the hydride anion is in fair agreement with the 
hydride anion radii calculated in the saline 
hydrides if the electronegativity of the cation 
is considered. The hydride anion radius de- 
creases from about 1.53 A to 1.37 A in the series 
CsH to LiH, i.e., with increasing electro- 
negativity of the cation. Since the electro- 
negativities of the transition metals are all 
higher than lithium, a value for the radius of 
the hydride anion somewhat less than 1.37 A 
would be expected. 

Magnetic susceptibility studies of transition 
metal hydrides have usually been cited as 
evidence for the formation of protons on 
hydriding. Fitzwilliam et al.12°) have observed 
that the paramagnetic susceptibility of zir- 
conium decreases on the addition of hydrogen, a 
minimum occurring at H/Zr equal to two. 
Similar observations have been made on the 
palladium-hydrogen!), niobium-hydrogen 121), 
tantalum-hydrogen 7”) and lanthanum- hydro- 
gen 122) systems. In each case, the magnetic 
susceptibility decreases with increasing hydro- 
gen content, in palladium a constant value being 
reached at H/Pd=0.7. (It should be mentioned 
here that susceptibility measurements on the 
titanium-hydrogen system show an increase in 
paramagnetic susceptibility with composition!), 
The decrease in magnetic susceptibility with 
hydrogen content for most hydrides has been 
explained by the filling of holes in the d-band 
of the metal with electrons from the hydrogen 
atoms to form protons. However, the pairing 
of electron spins in the formation of hydride 
anions will also give decreasing magnetic 
susceptibility with increasing hydrogen content. 

The observed migration of protons in the 
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palladium-hydrogen 123), tantalum -hydrogen 
and titanium-hydrogen 124) systems has often 
been cited as evidence for the existence of 
protonic hydrogen. However, all these measure- 
ments were carried out at hydrogen contents 
low enough so that only the metallic phase 
containing hydrogen in solution was present. 
To date, there have been no results available on 
hydrogen migration in an electric field in the 
single phase hydride region. 

In some preliminary calculations by Gibb 15), 
in which the ability of a metal to form a hydride 
is deduced from lattice energies, good agree- 
ment with experimental fact was obtained for 
lattice energies calculated on the basis of an 
ionic model for transition metal hydrides. 


5. Conclusion 


A review of the properties of transition metal 
hydrides shows that, in general, the hydrides 
within any one group of the periodic table have 
similar properties. The Group IV A metals all 
form hydrides with a maximum hydrogen to 
metal ration of two. These hydrides all have a 
cubic fluorite structure at lower hydrogen 
contents and become increasingly tetragonal as 
the hydrogen to metal ratio approaches two. 
This effect is most pronounced in zirconium 
hydride and least in titanium hydride. 

The rare earth metals form cubic trihydrides 
in addition to dihydrides having the cubic 
fluorite structure (except ytterbium and euro- 
pium). The trihydride takes on a hexagonal 
structure for the case of gadolinium and 
hydrides of higher atomic number. On the other 
hand, the actinide series hydrides do not seem 
to have much in common with each other. The 
plutonium hydrides resemble the hydrides of the 
rare earth metals of higher atomic number, but 
the hydrides of uranium and thorium are quite 
different, although the lower hydride of thorium 
ThHe does have a fluorite structure. There has 
been insufficient work on the other hydrides in 
this series to draw any conclusion as to simi- 
larities in properties. 

The hydrides of the Group V A metals form 
hydrides approaching the stoichiometric ratio 


MH, and have structures identical or similar to 
the parent metal. The only other transition 
metals which appear to form hydrides are 
palladium, chromium, and copper and some 
alloys of hydride-forming with non-hydride- 
forming metals. 

In conclusion, it should be emphasized that 
the hydrides of the transition metals, in particu- 
lar those of Group III A and IVA metals 
(including lanthanides and actinides), are defi- 
nite chemical compounds rather than interstitial 
solid solutions. The large deviations from stoichio- 
metry in these compounds can be attributed to 
lattice defects as in most transition metal oxides 
or sulfides. The metallic hydrides appear to have 
some degree of ionic character along with their 
metallic properties. However, there is obviously 
need for further study in order to fully under- 
stand the nature and bonding of these inter- 
esting compounds. 
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The transformation temperatures of high-purity 
uranium were determined by repeated thermal 
analyses where the rates of heating and cooling were 
independently controlled. The solid state trans- 
formation temperatures and the logarithm of the 
heating or cooling rates are related linearly. The 
extrapolated functions intersect at a point where the 
disturbing effects of hysteresis, superheating and 
undercooling disappear, i.e. at the equilibrium 
temperature. The transformation temperatures are 
667.7 +1.3°C for «o =f and 774.8 + 1.6°C for 
6 =y. The mean temperature for melting and freezing 
is 1132.3 + 0.8° C. 


Les températures de transformation de uranium de 
haute pureté ont été déterminées par des analyses 
thermiques répétées en contrdlant indépendamment 
les vitesses de chauffage et de refroidissement. Les 
températures de transformation dans l’état solide 
varient linéairement en fonction du logarithme de la 
vitesse de chauffage ou de refroidissement. Les 
fonctions extrapolées se coupent en un point ou les 


1. Introduction 


Measurements of the solid state transfor- 
mation temperatures of uranium were made as 
early as 1942 and of its melting point in 1930. 
The early data were reviewed by Katz and 
Rabinowitch1). Few of these measurements 
have any significance today since they were 
made on rather impure metal under conditions 
where the formation of large quantities of oxide 
and nitride was inevitable. Nevertheless, Chip- 
man 2) reported in 1946, in summarizing the 
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effets perturbateurs de Vhystérésis de surchauffe ou 
de sous-refroidissement disparaissent, c’est a dire a 
la température méme de l’équilibre. Les températures 
de transformation sont 667.7 + 1.3° C pour « = 6 et 
774.8 + 1.6°C pour Bf =y. La température moyenne 
pour la fusion et la solidification est de 1132.3 + 0.8° C. 


Die Umwandlungstemperaturen von reinstem Uran 
wurden mit Hilfe wiederholter, thermischer Analysen 
unter unabhangiger Kontrolle der Erhitzungs- und 
Abkithlungsgeschwindigkeit bestimmt. Die Umwand- 
lungstemperaturen des festen Zustandes und der 
Logarithmus der Erhitzungs- und Abkthlungs- 
geschwindigkeit sind linear von einander abhangig. 
Die extrapolierten Kurven schneiden sich in einem 
Punkt, wo die st6renden Hinfliisse von Hysterese, 
Uberhitzung und Unterkiihlung verschwinden; dieser 
Punkt ist die Gleichgewichtstemperatur. Die Um- 
wandlungstemperaturen sind 667.7 + 1.3° C flra =f 
und 774.8 + 1.6° C fur B = y. Das Mittel der Schmelz- 
und Erstarrungstemperaturen is 1132.3 + 0.8° C. 


work of the Metallurgical Laboratory, an « = 6 
transformation temperature of 665° C, a B = y 
transformation temperature of 775°C and a 
melting point of 1130° C; data which approach 
present day measurements. 

In making heat content measurements, 
Moore and Kelley ?) found two sharp transitions 
at 662°+4+ 3°C and 772°+ 3°C with no 
evidence of hysteresis. Resistivity measure- 
ments by Dahl and Van Dusen 4), however, 
revealed a large hysteresis for the solid state 
transformations. Their data were: 


« —> B=667° C 
p > «=645° C 


B= y=772°C 
y > B=764° C. 


The data by Dahl and Van Dusen were 
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evaluated by Ginnings and Corrucini *) for the 
purpose of calculating heats of transition. They 
used 668° C and 774°C, probably the highest 
values obtained by Dahl and Van Dusen. 
Similar results were obtained by some of the 
investigators of binary phase diagrams §~13), 

While these results appeared to be fairly 
consistent, they were, nevertheless, looked upon 
with suspicion partly because of the impurity 
content of the metal and partly because of the 
large hysteresis effects, even at low rates of heat- 
ing and cooling. The very large effect of the 
cooling rate on the transformation temperatures 
was shown by Duwez 14) who, at a rate of 
8000° C/sec, obtained supercooling of the y > 6 
and $-—>« transformations by about 200° C 
and 300° C, respectively. 

Only one reliable datum existed for the y = L 
phase change. Dahl and Cleaves ) determined 
the freezing point (Ly only) by taking a 
sequence of cooling curves at a rate of 1° C/min. 
Successive determinations resulted in an in- 
crease of the freezing point from 1125°C to 
1133° C. Simultaneously, a reduction of the 
carbon content was observed in the melt; the 
highest freezing point was associated with a 
carbon content of 6 ppm. Most of the other 
observers reported data below 1133° C without 


specifying the carbon content or giving other 
analytical data. With the advent of high 
purity uranium 16, 17,18) and the development 
of precise analytical methods for impurity 
determinations it became possible, and indeed 
desirable, to determine the transformation 
temperatures. 

Differential thermal analyses on high purity 
uranium were made by Baumrucker 19). He 
found 7m=1130.9 + 0.4°C at rates ranging 
from 0.60° to 3.88° C/min for melting and 1.13° 
to 2.80° C/min for freezing. For « = f he found 
666.8 + 0.4° C for heating and 657.7 + 0.5° C 
for cooling at variable rates ranging from 0.8 
to 2° C/min, with a hysteresis of 9.1° C. For 
6 — y he found 772.2 + 0.4° C for heating and 
767.8° + 1.5° C for cooling at variable rates of 
0.72 to 4.5° C/min, with a hysteresis of 4.4° C. 
Dilatometric measurements of the solid state 
transformation temperatures at the constant 
rate of 1°C/min by Lloyd 1%), showed much 
less hysteresis than those by Baumrucker. 
Lloyd found 662.6 + 0.9° C for « + 6, 660.5 + 
0.6° C for B + «, 764.0 + 1.0° C for 6B + y and 
768.0 + 1.2 for y—+. Lloyd’s dilatometric 
data are lower than those by Baumrucker. 
Although these data were obtained on high 
purity uranium, the disturbing undercooling 


TABLE 1 


Composition of high purity uranium 


Remarks 


Spec. Chemical analysis in ppm Spectrochemical analysis * in ppm 
ee AH C N GO| | Al Cr Cul Beethte ain es 
B-721 | — | 26-35 <10 — creel gt Oe alin Meal 10 
— | 28-37 | <10-19 1 Ui el 2 yell 20 
B-724 | — | 26-35 <10 — hs i Perl 10 
— | 22-25 <10 1 7 2 2 <1 10-15 
B-726 | — | 34-40 <10 6.6 Opies 2 Ue ell 20 
— | 40-50 <10 2 Tee Ue Geil 10 
B-728T | .21 5 <10 22-25 5 1 1 We Mele al 15 
B-734 | — 7-10 12 25-26 Dig el ae 1 <l 10 


Before thermal analysis 
After thermal analysis 
Before thermal analysis 
After thermal analysis 
Before thermal analysis 
After thermal analysis 
Used in solid state only 


Used in solid state only 


eS 


* 


All other elements below limits of spectrochemical detection. 
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and superheating effects persisted and called 
for a closer control of the rate factor. It became 
necessary to clear up these questions and 
establish precise equilibrium temperatures as a 
preliminary to a study of the constitution of 
the low carbon U-C alloys by B. Blumenthal, 
where the effect of carbon on the transformation 
temperatures was expected to be small. 


2. Materials 


The material used in this investigation was 
remelted high purity electrolytic uranium; the 
pertinent analytical data are given in table 1. 


3. Experimental Apparatus 


The experimental apparatus was a modified 
version of the high vacuum Globar resistance 
furnace described previously 1’). Its vacuum 
system was left unchanged. The heating system 
was rebuilt to permit repeated thermal analyses 
at a predetermined, yet variable rate unaffected 
by load oscillations caused by the common on- 
off controls. The furnace was heated by two 
separate sets of Globars; one provided the 
manually adjusted base load (about 2/3 of the 
total), the other one the control load. The latter 
was supplied by a Variac, which was driven by 
an air operated diaphragm motor (Conoflow 
Corp.) in conjunction with a Minneapolis- 
Honeywell triple range pneumatic controller 
(2 to 8, 6 to 12, and 10 to 16 MV). The program 
control was imposed by a cam cut to suit the 
desired cycling pattern. Since the control 
thermocouple was touching the outside of the 
furnace tube, it did not register the true 
temperature inside the furnace. The temperature 
difference due to the changing conditions of the 
experiments was about 25°C. To compensate 
for this difference the cam was cut for a 
temperature that much higher. With this 
compensation it was possible to impart any 
desired temperature pattern to the specimen. 
Since the furnace had a low heat capacity, it 
followed the predetermined cycle with little lag. 

For melting point determinations, a specimen 
of about 300 grams was contained in a urania 


crucible; a Pt/Pt-10 % Rh thermocouple pro- 
tected by a urania tube was located in its 
center. Care was taken to make sure that the 
thermocouple was secured in its central position 
throughout the melting and freezing cycles. 

The rather rigid suspension is shown in fig. I. 
For the solid state transformations a similar 
quantity of bare metal was suspended in 
the furnace as shown in fig. 2. The cold 
junctions of both the measuring and the control 
thermocouples were held at 0° C in well insulated 
Thermos bottles. Suitable precautions were 
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Fig. 1. Suspension of crucible for thermal analysis. 
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Fig. 2. Suspension of specimen for solid state 
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Fig. 3. Effect of cooling rate on the y—>f and 

Bf +o transformation temperatures of high-purity 

uranium (at left data by Blumenthal, at right data 
by Duwez). 


taken to shield crucible and melt from the cold 
sections of the furnace and to establish a uniform 
temperature gradient around the melt or the 
solid specimen. Despite the sensitivity of the 
urania refractories to thermal and mechanical 
stresses, the arrangement proved to be highly 
successful provided the system was not cooled 
to room temperature during an experiment. 
Many cycles at various rates in the range of 
+ 40°C around the melting point could be 
made over a period of days. A pressure of 
5 to 7x 10-7? mm Hg was maintained throughout 
each series of cycling experiments. 

The specimen temperature was recorded by 
a specially built Leeds and Northrup Speedomax 
recorder with a span of 3 MV and ranges 
beginning at 0, 2.5, 5, 7.5, and 10 MV. Its 
sensitivity was 0.003 MV and its error within 
0.009 MV. The instrument was calibrated for 
recording of the correct EMF and for mechanical 
differences in recording at rising or falling 
temperatures, before and after each experiment 
by means of a Rubicon potentiometer. In the 
later runs the paper position was checked by 
continuously tracing a straight line with a sepa- 
rate, fixed pen. Under these circumstances it was 
possible to read the EMF within + 0.005 MV. 

The measuring thermocouple was calibrated 
before and after each group of experiments over 
the range of 500° to 1250° C against a standard 
couple which had been calibrated by the 
National Bureau of Standards. The corrections 
applied to the actual temperature measure- 


ments were obtained from smoothed out curves 
of these two calibrations. Care was taken to 
obtain a strain free thermocouple by annealing 
it in air at 1500° C. Thus four corrections were 
applied to each individual EMF measurement: 
one for errors in instrument recording, a second 
for errors due to paper position, a third for the 
calibration of the measuring thermocouple 
against the standard couple and a fourth for 
the calibration of the standard couple by NBS. 
Ultimately, the EMF was converted to degrees 
centigrade using the tables of NBS Circular 
No. 5610-45 

Occasionally a heating or cooling curve, which 
generally confirmed the findings of the instru- 
ment, was made by direct potentiometric 
measurement. 


4. Data 


4.1. MELTING AND FREEZING 


The data are given in table 2. The results 
were obtained at a constant rate of 0.74° C/min. 
Good arrests of about 22 minutes duration were 
obtained during the 7 to 11 cycles which were 
made on each specimen. Although the graphs 
produced by the automatic recorder gave no 
indication that superheating or undercooling 
took place, a direct potentiometric measurement 
on ingot B-727 showed very clearly that both 
effects did occur. The freezing temperatures 
were higher than the melting temperatures. The 
reason for this isnot known. It may be explained, 
however, by the presence of a temperature 
gradient between the melt and the thermocouple 
which is shielded by a rather poorly conducting 
refractory (urania). This temperature gradient 
persists through the long arrests and causes the 
measured temperature to be lower during 
heating and higher during cooling. It is assumed 
that the gradients were the same on heating and 
cooling. The mean of the melting and freezing 
temperature is then the equilibrium solid = 
liquid transformation point. This point is at 
1132.3 + 0.8° C. 

The difference of 1.4° C between the present 
results and those by Baumrucker cannot be 
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TABLE 2 


Melting and freezing temperatures of high purity 
uranium at a rate of 0.74° C/min 


Specimen Melting point | Freezing point 

No. (°C) (°C) 
1133.0 DSS ol 

1132.3 1134.5 

1132.3 1134.1 

B-721 1131.9 1133.8 
1132.0 1133.9 

1132.5 1133.9 

1133.6 

1131.8 1131.8 

1131.1 1131.3 

1130.9 1131-9 

B-724 1130.8 1131.9 
1131.4 1132.1 

| 1131.7 1131.8 

1131.6 1132.0 

1130.1 1133.3 

1130.3 1133.3 

1130.6 1133.8 

1130.8 1133.3 

1130.9 1133.4 

B-726 1131.0 1133.4 
1130.7 1133.8 

1130.6 1132.9 

1130.9 1132.8 

1133.6 

1133.0 

Mean of 3 ingots 1131.4 + 0.4 1133.1 + 0.4 
Mean of melting 
Pratmecing 1132.3 + 0.8 


explained easily. Since the melting temperatures 
agree within the limits of precision, the differ- 
ence is due to the difference in the observed 
freezing temperatures. It is probable that 
Baumrucker missed the recalescence during 
solidification and that his data are too low. 
Since the higher temperatures agree with the 
earlier measurements on less pure material by 
Dahl and Cleaves and generally higher purity 
materials give higher equilibrium temperatures 
than lower purity materials, the higher values 
are believed to be more nearly true than the 
lower ones. 


4.2. 

The very large hysteresis between heating 
and cooling, even at low rates, and the large 
effect of the rate on the magnitude of this 
hysteresis, made a new approach to the determi- 
nation of the solid state transformation points 
mandatory. A few direct potentiometric meas- 
urements showed clearly that not only under- 
cooling, but also superheating are real. Further- 
more, it is not at all sure that upon recalescence 
the latent heat evolved brings the specimen 
temperature back to the equilibrium temper- 
ature. Lag in response of the temperature 
sensing element is still another factor to be 
considered. The problem was, therefore, ap- 
proached by measuring the transformation 
temperatures on heating and cooling at various 
constant rates. The spread between heating and 
cooling of the « = £ transformation temper- 
ature decreased as the rate decreased. A plot of 
the transformation temperature versus the 
logarithm of the rate showed that the present 
data form a satisfactory continuation of 
Duwez’ 14) data for high purity uranium ob- 
tained at much higher cooling rates (fig. 3). The 
left end of the curve may be regarded as a 
straight line. By the method of least squares, 
functions of the type: 


THE ALPHA-BETA TRANSFORMATION 


Tr=a+b1 logior 
Tc = Ag — be logior 


were obtained, where 7'q and 7'¢ are the trans- 
formation temperatures in ° C on heating and 
cooling respectively, logior is the logarithm of 
the rate, and a, de, b; and be are constants. 
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Fig. 4. Effect of heating and cooling rate on the 
« =f transformation temperature of high-purity 
uranium (B-724 and B-728T). 
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The two lines intersect at the temperature 
T =Tg=T¢ which is the temperature where the 
disturbing factors of hysteresis, superheating 
and undercooling disappear and which may be 
regarded as the equilibrium temperature of the 
transformation under investigation. By this 
method the following equations and trans- 
formation temperatures were obtained: 

For « = 8 (Ingot B-724) 


T 4 = 671.79 + 0.064 logior 
To = 662.22 —0.151 logyor. 


The two lines intersect at 7’=669.0° C and at 
r=2.4x10-45° C/min. This very slow rate is 
the result of the metal having been contained 
in a urania crucible, a procedure that was 
subsequently abandoned. Afterwards only bare 
specimens were used as shown in fig. 2. 

For « = 6 (Ingot B-728T) 


T= 673.964 4.19 logior 
To =659.61—3.76 logior. 


The two lines intersect at 7’=666.4° C and at 
r= 0.016° C/min. The mean of the two values is 
667.7 + 1.3°C (fig. 4). The details of the 
measurements are given in table 3. 

Baumrucker’s mean value 666.8 + 0.4° C for 
heating agrees with this value within the limits 
of precision. 


4.3. THE BETA-GAMMA TRANSFORMATIONS 


The same extrapolation method was used to 
determine the §f =y transformation temper- 
ature (table 4). The two specimens (B-728T 
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Fig. 5. Effect of heating and cooling rate on the 
6B =y transformation temperature of high-purity 
uranium (B-734). 


and B-734) underwent superheating in a manner 
comparable to that observed for the « > 6 
transformation. The transformation tempera- 
tures on heating versus the logarithm of the 
rate may again be expressed by the relationship 
T= +b; logior. The pertinent equations are: 


B-728T 
B-734 


eee T y= 778.34 + 2.34 logior 
aac Ty=779.11+2.19 logior (fig. 5). 


TABiE 3 


Mean « =f transformation temperatures of high 
purity uranium 


a> fp p>« 
ye Geo neers No. of No. of 
No. (° C/min) Tome single wees single 
( C) meas. ( ) meas. 
2.0 672.3 3 661.3 3 
2 AO) 671.7 9 662.1 9 
Bete 0.667>-| 6714 | @ | 66at-\) 4 
0.50 671.5 2 663.0 Ps 
4.0 677.3 3 657.7 3 
2.0 675.7 8 658.5 8 
1.0 672.6 4 659.4 4 
B-728T 0.50 671.6 3 660.0 2 
0.333 GTIRT 4 661.4 4 
0.167 672.1 2 663.1 2 
TABLE 4 


Mean $6 =y transformation temperatures of high 
purity uranium 


. oy y>B 
Spogmnen pate No. of No. of 
No. (°C/min)| Temp. | |. Temp, | 
x single Z single 
ge) meas. | ' ©) meas. 
B-728T 2.47 780.2 3 774.7 3 
0.825 777.2 8 772.3 8 
0.412 didgied| 3 777.6 3 
0.206 775.7 3 772.4 3 
0.103 Ugoe 3 775.0 3 
B-734 2.47 780.3 4 774.9 4 
0.825 778.0 6 776.2 6 
0.412 779.0 3 775.8 3 
0.206 777.6 4 774.7 4 
0.103 776.9 3 774.8 3 
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No rate effect of the kind that was observed 
for the 6 > « transformation took place on 
cooling although recalescence was noted in 
several, but not all, cycles. Thus, the cooling 
data are best represented by the mean of the 
observed values which graphically form a 
straight line parallel to the logio 7 axis. 

The intersection of the lines for heating and 
cooling and, therefore, the equilibrium trans- 
formation temperature is identical with the 
mean of the cooling curve data. They are: 


Po 7ST Me TIED a 2:90) 
Pom Toe 7153 4 0.7° C (fig. 6) 


The 6 = y transformation temperature of high 
purity uranium is 774.8 + 1.6°C. 

This temperature is 7°C higher than the 
mean of Baumrucker’s cooling data, most of 
which were obtained at rates of 1.6° C/min and 
higher. It is possible that the tantalum thermo- 
couple protection tube which was used by 
Baumrucker conducted enough heat away from 
the specimen to produce the lower value. 

Similarly it is probable that heat transfer and 
heat conduction phenomena greatly influenced 
the measurements of Lloyd and account for the 
large difference even at comparable rates. On 
the other hand it is interesting to note that in 
the 6 — « case when undercooling and recales- 
cence are small and the insulating properties of 
the sleeve in the dilatometer had a favorable 
effect, the agreement between the dilatometric 
measurements of Lloyd and the present results 
is quite good. 


4.4. THE BETA PHASE ANOMALY 


The Beta phase anomaly encountered in the 
measurements of Duwez and Lloyd merits 
mention. A heat evolution was noted in the 
beta region on cooling. For the high purity 
specimens this occurred at temperatures ranging 
from 686° to 728° C, it also occurred in several 
of the runs made with impure materials. The 
phenomenon was characterized by a sudden 
increase in temperature resulting in a momen- 
tary slight expansion of the specimen, but with 


no change in the slope of the cooling curve. 
Other investigators have reported a break in the 
magnetic susceptibility curve of uranium at 
about 698° C. 

In spite of a diligent search, no such anomaly 
was detected in the present thermal analyses 
A similar search by Baumrucker also had 
negative results. 

Since Blumenthal and Baumrucker used cast 
material of large grain size and Lloyd and Duwez 
fabricated material of much smaller grain size, 
it is possible that the effect is associated with 
the release of some strain energy introduced 
during the y — f transformation. It is felt that 
the effect merits further investigation. 


5. Conclusion 


Because of the sensitivity of the solid state 
transformation to the rate at which measure- 
ments are made and because of the large effect 
which heat transfer and heat conduction 
phenomena have on the temperature measure- 
ments of all uranium transformations, only 
those measurements should be considered valid 
which take these factors into account. It is, 
therefore, felt that the present measurements 
and calculated transformation points come 
closest to true equilibrium conditions. The 
equilibrium transformation temperatures of 
high purity uranium are: 


x= ps 6677 Se V3 € 
pay: TI48+4£1.6°C 
yes Te: 1132.3 0.82 C. 
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Experiments have been done to determine the kinetics 
of the oxidation of graphite at 900—-1000° C in argon 
containing small amounts of oxygen, water vapour, 
* carbon dioxide, hydrogen and carbon monoxide. 

The reactivity of the gases decreased in the order 
oxygen, water vapour, carbon dioxide. The reaction 
rates increased with temperature and with the partial 
pressure of the oxidising gas, whilst hydrogen and 
carbon monoxide inhibited attack by water vapour 
and carbon dioxide respectively. 

The relationships between the reaction rates and 
the partial pressures of the gases are in good agreement 
/ with theory, previous work and the supposition that 
1 hydrogen atoms are adsorbed on graphite in contact 

- with hydrogen and water vapour. 


Des expériences ont été réalisées pour déterminer la 
cinétique de l’oxidation du graphite a 900-1000° C 
dans de l’argon contenant de petites quantités d’oxy- 
géne, de vapeur d’eau, de gaz carbonique, d’hydrogéne 
et d’oxyde de carbone. 

La réactivité des gaz diminuait dans l’ordre oxy- 
gene, vapeur d’eau, gaz carbonique. Les vitesses de 
réaction croissaient avec la température et avec la 
pression partielle du gaz oxydant tandis que l’hydro- 


1. Introduction 


The present design of the high temperature 
gas cooled reactor experiment to be built at 
Winfrith, England, has a graphite moderator 
and a fuel element of uranium and thorium 
carbides clad in a graphite can. The maximum 
temperature of the can is approximately 1000° C 
and the coolant is helium. Although the coolant 
is inert, compatibility problems exist as oxygen, 
carbon dioxide and water vapour, present as 
impurities, may react with the graphite to give 
carbon monoxide and hydrogen. The relevant 
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géne et loxyde de carbone inhibaient l’attaque par 
la vapeur d’eau et le gaz carbonique respectivement. 

Les relations entre les vitesses de réaction et les 
pressions partielles des gaz sont en bon accord avec 
la théorie, les travaux antérieurs sur le sujet et aussi 
avec ’hypothése d’une absorption des atomes d’hydro- 
géne sur le graphite en contact avec lhydrogéne et 
la vapeur d’eau. 


Die Kinetik der Oxydation von Graphit in Argon, dem 
kleine Mengen Sauerstoff, Wasserdampf, Kohlen- 
dioxyd, Wasserstoff und Kohlenmonoxyd zugesetzt 
waren, wurde bei 900—-1000°C untersucht. Dabei 
wurde der Angriff bein Ubergang von Sauerstoff auf 
Wasserdampf und dann auf Kohlendioxyd schwacher. 

Die Reaktionsgeschwindigkeit stieg mit wachsender 
Temperatur und mit zunehmendem Partialdruck des 
oxydierenden Gases an. Wasserstoff und Kohlen- 
monoxyd hemmten den Angriff von Wasserdampf 
bzw. Kohlendioxyd. 

Die Beziehung zwischen der Reaktionsgeschwindig- 
keit und dem Partialdruck der Gase steht in gutem 
Hinklang mit theoretischen Vorstellungen, mit friihe- 
ren Ergebnissen und mit der Annahme, dass bei 
Gegenwart von Wasserstoff und Wasserdampf Wasser- 
stoffatome am Graphit absorbiert sind. 


literature has been summarised and out-of-pile 
experiments done to determine the kinetics of 
the attack of graphite, in an inert gas containing 
small amounts of reactive impurities. 
Previous work on sorption and_ surface 
reactions has demonstrated the ability of carbon 
to retain oxygen atoms at certain active sites. 
Typical active sites are believed to be impurity 
atoms and the edges of graphitic planes. On 
an initially clean carbon surface oxygen atoms 
are adsorbed irreversibly from oxygen whilst 
carbon dioxide and water vapour are reduced 
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reversibly to carbon monoxide and hydrogen. 
These reactions may be expressed as 


20;+02 er 2Co 
and 
CO Het CO 
Cry aS Co+H,, (1) 


where Cy represents a free site and Co a site 
possessing an oxygen atom. k; and kz are the 
relevant rate constants. Previous work indicates 
that even at temperatures as low as 600° C 
exchange reactions are fast. 

Transfer of carbon to the gas phase may take 
place by the release of carbon monoxide from 
occupied sites and may be expressed as 


kg 
Op = SS CO] ae 


where » is an integer which normally has a 
value 1 but can be between 0 and 21). 

The rate of loss of carbon (dm/dt) can be 
expressed as 

am — ks (Co) A, (2) 

where (Co) is the number of occupied sites per 
unit area and A is the available surface area. 

Under steady state conditions the rate of 
formation of occupied sites will be equal to 
their disappearance. 


Bee a ky Po (C+) 2s ko Pr(Co) —— k3(Co) = 0, 


where po=the partial pressure of the oxidising 
gas 1.e. O2, COs, H2O. pr=the partial pressure 
of the reducing gas i.e. CO, He. 

If Cy is the total number of reaction sites 
per unit area 


(Ct) = (Cr) + (Co) 
and 


j= \ ees ki po Ct 
(Co) = kipo+kapr+ks° ) 
Substitution of equation (3) in (2) yields 


dm = kz ky Po CiA 4 
dt ki po +ke pr-+ks (4) 


It follows that if the products are removed 
the steady rate of carbon loss will be propor- 
tional to the partial pressure of the oxidising 
gas if kypo is small compared with k2p;y+ ks, 
and that the reciprocal of the reaction rate 
should be proportional to the partial pressure 
of carbon monoxide and hydrogen. 

Although a vast amount of data has been 
published on the oxidation of coals, cokes, etc. 
little is applicable to the oxidation of a pure 
graphite at 1000° C in an inert gas containing 
small amounts of oxidising impurities. The 
most relevant work concerns the oxidation of 
pure carbon filaments in low partial pressures of 
oxygen, water vapour and carbon dioxide 2: 3; 4). 
At most temperatures the reactivity of the gases 
decreases in the order oxygen, water vapour, 
carbon dioxide. The reactivity of graphite 
increases with temperature up to 1100-1400° C, 
and then steadily decreases to 1700° C due to 
the removal of active sites by an annealing 
process. 

There is some confusion as to the effect of 
pressure and the order of the reactions but most 
workers are agreed that oxygen, water vapour 
and carbon dioxide behave similarly; the rates 
increase with pressure and the reactions are 
first order at 1200°C and zero order at 
800-900° C. At the lower temperatures all the 
available active sites are therefore occupied and 
the dependence of the reaction rate upon the 
chemical reactivity of the gas must be explained 
by differences in the activity and hence availa- 
bility of sites. With a first order relationship 
the rates are directly proportional to pressure 
whilst with a zero order the gas penetrates the 
pores and the rates may depend upon pressure 
to a fractional powder. 

The most reliable work appears to be that 
of Duval?) and Boulangier 3) who interpreted 
their results from values for the ratio of number’ 
of gas molecules reacting to the number striking 
the carbon. The reaction rate (r) is related to 
this ratio (Ff) by the equation 


r= 3.513 x 102 -2_ x F molecules/sec em2 
VMT 
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where p= partial pressure of the gas in mm Hg. 
M=molecular weight of the gas 
T =absolute temperature of the gas 
r=number of gas molecules reacting/ 
sec cm2. 


F was invariant with pressure at 1000° C for 
carbon dioxide and oxygen over the range 
10-3-10-2 mm Hg and therefore the rate was 
directly proporticnal to pressure and _ the 
reaction first order. The dependence of rate on 
water vapour pressure was complicated as F 
decreased with increasing pressure over the 
range 10-4-10-2 mm Hg; the effect was ex- 
plained by postulating the adsorption of 
hydrogen atoms on active sites. 

Strickland-Constable 4) with 10-°-10-! mm 
Hg of water vapour and carbon dioxide 
obtained a fractional order reaction at 900- 
1100° C changing to first order at 1200° C. 

Blyholder and others >) showed that although 
the reaction with oxygen at pressures of 
10-8-10-! mm Hg and temperatures of 800- 
900° C was zero order the reaction rate varied 
with pressure to the one half power due to 
appreciable diffusion of gas into the pores of 
the graphite. 

Sabri Ergun 1) has demonstrated that carbon 
monoxide inhibits the attack of graphite by 
carbon dioxide and eq. (4) was found to explain 
the attack by CO-COz gas mixtures at one 
atmosphere pressure over the temperature 
range 700-1400° C. 


2. Experimental Procedure 


Two sets of experiments have been done, in 
the first, rates of attack by 5-10 vpm oxygen, 
0-1000 vpm water vapour and 0-500 vpm 
carbon dioxide in one atmosphere of argon 
were determined for three types of graphite at 
900-1000° C; in the second the effects of 0-3000 
vpm hydrogen on the oxidation by 300 vpm 
water vapour and of 0-5000 vpm carbon 
monoxide on the oxidation by 500 vpm carbon 
dioxide were determined for Pile Grade A 
graphite at 1000° C. 

A conventional silica balance was used to 


measure rates of weight loss for the samples 
of graphite in a flowing stream of argon con- 
taining known amounts of the reactive gases. 
The graphite specimens were parallelepipeds 
16x5x5 mm; they weighed 700 mg and were 
suspended from the spring on a thin platinum 
wire. Great care was taken not to contaminate 
the graphite: the machining was done at a 
factory specialising in pile grade graphite and 
experimentalists always wore gloves when 
handling the specimens. Three different graphites 
were used in the experiments, namely Pile 
Grade A, EY9 (made by Morgan Crucible Co.) 
and EY9 which had been impregnated with a 
furfuryl alcohol resin and heat treated to 
approximately 1000°C to make it almost 
impervious to gas. They had densities between 
1.7 and 1.8 g/cc. The gas was at atmospheric 
pressure and passed through the balance at 
100-300 std. cc/min. The total burn-off during 
the course of the experiments varied between 
20 and 150 mg, with the most reactive samples 
having the greatest burn-offs. 

Argon sufficiently pure for all the experiments 
was obtained by passing the gas through 
manganous oxide at 150° C to lower the oxygen 
content to <1 vpm followed by anhydrous 
magnesium perchlorate to lower the water 
vapour content to <15 vpm. The oxides of 
carbon were added to the gas by passing it 
successively through heated beds of magnesium 
carbonate and graphite. Hydrogen was added 
by a conventional flowmeter technique, and 
water vapour by passing the argon-hydrogen 
mixture through a heated bed of cupric oxide 
or by passing a proportion of the main gas flow 
over a saturated solution of calcium chloride. 
When heated beds were used their temperature 
controlled the partial pressure of the additive. 
Experiments were only done with the 5-10 vpm 
of oxygen present as an impurity in cylinder 
argon and injection techniques were therefore 
not needed for this gas. 

Oxygen was analysed by a Hersch meter, the 
oxides of carbon by an Infrared Gas Analyser 
and water vapour by an Electrolytic Water 
Analyser. 
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TABLE 1 


Attack of graphite by gaseous impurities 


nnn nn EEE 


: Temp Tpurity cones -VRia Rate of weight Flowrate 
paces 2G Os COs H.O rare) iH; loss mg/cm? h (std. cc/min) 
Pile Grade A 1000 5-10 <a Nee OD 0.06 + 0.002 150 
EY9 1000 5-10 | eos 0.025 + 0.002 75 
EY9 1000 5-10 OT VD 0.039 + 0.002 300 
EY9 1000 << <0 <a 0.003 + 0.002 
EY9 1000 <1 <li) 420 0.20 + 0.002 
EY9 1000 eae < 10 720 0.43 + 0.002 
EY9 impregnated 1000 =e <M eal) 0.003 + 0.002 
EY9 impregnated 1000 ant < 10 76 0.027 + 0.002 
EY9 impregnated | 1000 ull Ze IKY 167 0.11 + 0.002 
EY9 impregnated 1000 <l <a li () 315 0.26 + 0.002 
EY9 impregnated 1000 << ll < 10 545 0.51 + 0.002 
Pile Grade A 1000 <a} LOM al 5 0.003 + 0.002 
Pile Grade A 1000 —<el =a) 75 0.023 + 0.002 
Pile Grade A 1000 <1 < 10 238 0.117 + 0.002 
Pile Grade A 1000 <l Ean) 508 0.234 + 0.002 
Pile Grade A 900 < Il IO cals 0.003 + 0.002 
Pile Grade A 900 <l1 <<) 103 0.011 + 0.002 
Pile Grade A 900 <i < 10 250 0.034 + 0.002 
Pile Grade A 900 <I << II) 363 0.052 + 0.002 
Pile Grade A 900 ell < 10 654 0.096 + 0.002 
Pile Grade A oOe! LTE yr) ess 0.11 + 0.002 
Pile Grade A 1000 <1 MO ce 1s 0.003 + 0.002 
Pile Grade A 1000 < Il 250 ely 0.014 + 0.002 
Pile Grade A 1000 <i 500 | < 15 | 0.024 + 0.002 
EY9 1000 << <PLON a5 < 0.001 
EY9 1000 <u 260 | < 15 0.047 + 0.002 
EY9 1000 8 | 5380 | < 15 0.136 + 0.002 
Pile Grade A 900 <I << OR ee<ean ted < 0.002 
Pile Grade A 900 <<ell 270 | < 15 0.004 + 0.002 
Pile Grade A 900 <= i 500 | < 15 0.009 + 0.002 
Pile Grade A 1000 <1 500 | < 15 250 — 0.018 + 0.002 
Pile Grade A 1000 <el 500 | < 15 500 — 0.014 + 0.002 
Pile Grade A 1000 <i 500 | <15 | 5000 — 0.004 + 0.002 
Pile Grade A 1000 —<gl < 10 300 —_ 300 0.067 + 0.002 
Pile Grade A 1000 <a) < 10 300 as 750 0.044 + 0.002 
Pile Grade A 1000 << ll < 10 300 — 1500 0.021 + 0.002 
Pile Grade A 1000 aul < 10 300 — 2100 0.013 + 0.002 
Pile Grade A 1000 <l < 10 300 — 3000 0.007 + 0.002 
3. Results reaction rates obtained with the various gas 
A visual examination of the specimens compositions are summarised in table 1. The 


indicated that most of the attack took place 
at the geometric surface, although the pores 
of the Pile Grade A specimens were slightly 
enlarged. 

The weight loss-time graphs were linear for 
a given set of experimental conditions and the 


relations between reaction rate and impurity 
concentration for water vapour and carbon 
dioxide together with the effects of hydrogen 
and carbon monoxide are shown in figs. 1-4. 

Although the attack by carbon dioxide and 
water vapour was independent of flowrate, the 
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Fig. 1. Attack of graphite in argon containing water vapour. 
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Fig. 2. Attack of graphite in argon containing carbon dioxide. 
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rates for gas containing 5-10 vpm oxygen 
varied with flowrate and corresponded to all or 
a large proportion of the impurity reacting, and 
were therefore controlled by the oxygen supply. 


4, Discussion 

All the results are consistent with theory and 
previous work in so far as the reactivity of the 
impurities decreased in the order oxygen, water 
vapour, carbon dioxide, and the rates increased 
with temperature and the partial pressure of 
the oxidising gas, whilst hydrogen and carbon 
monoxide inhibited the attack. 

The reactivity of the different graphites 
decreased in the order: impregnated EY9, EY9 
and Pile Grade A. From considerations of the 
methods of manufacture, the difference between 
EY9 and Pile Grade A is most likely due to 
impurities in the EY9, although it could be 
explained by differences in surface roughness due 
to differences in the internal surface area, the 
respective values of the internal area being 3-5 
and 0.5 m2/g (according to B.E.T. measure- 
ments). The high rates for impregnated EY9 
are probably due to the reactive nature of the 
char obtained from furfuryl alcohol. 

The graphs of rate versus water vapour 
content at 1000°C are linear for the greater 
part of the concentration range but change in 
slope or become curved in the range 0-70 vpm. 
The change in slope or curvature could be due 
to inaccurate gas analysis although this is 
considered unlikely; it is best explained by a 
film of hydrogen atoms, similar to that postu- 
lated by Boulangier 3), inhibiting the attack. 
The flattening off of the graph for Pile Grade A 
at 900° C with high water vapour contents is 
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indicative of a zero order reaction and pene- 
tration into the pores. 

With carbon dioxide at 900-1000° C the rate 
appears to be directly proportional to partial 
pressure although the results cannot be con- 
sidered sufficiently accurate to prove that a 
fractional order relationship does not hold. 

A comparison between the rates of attack in 
the present and previous work may be made if 
it is assumed that the reactions are first order 
and little affected by the inert gas. The partial 
pressures in the experiments of Duval, Bou- 
langier, etc. correspond to an impurity concen- 
tration of 5-10 vpm in a carrier gas at one 
atmosphere pressure and therefore if the rates 
for water vapour and carbon dioxide are 
linearly extrapolated through the origin a 
comparison may be made for the higher 
impurity concentrations used in the present 
work. An examination of table 2 reveals no 
outstanding discrepancy between the various 
sets of data. 

The theoretical reasoning leading to eq. (4) 
predicts that the reciprocal of the reaction rate 
should be proportional to the partial pressure 
of the reducing gas. The accuracy of the results 
for the CO2-CO additions is not sufficiently 
good to prove the relationship but an indication 
of the validity of eq.(4) may be obtained for 
the H2,O-—He additions. The reciprocal of the 
rate for 300 vpm water vapour is plotted against 
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the hydrogen concentration in fig. 5; the graph 
is not linear but curved upwards indicating 
that hydrogen has a greater effect than theory 
predicts. It can be shown that with the excep- 
tion of the rate for 300 vpm hydrogen the 
results satisfy an equation of the form: 


1 
rate 


= Ap 3+0 
where A and C are constants. 
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Fig. 5. Effect of hydrogen on the reciprocal of the 
reaction rate for the attack of Pile Grade A by 300 
vpm water vapour at 1000° C. 


TABLE 2 


Comparison of rates at 1000° C 


Rate of attack Rate 
; corresponding 
es oad weight loss impurity to 100 vpm 
mg/cm? h content vpm mg/cm? h 

Boulangier ) COz SielOne 10 Sie l0m 

Strickland—Constable *) CO2 LAX Os 10 RY Moll SK AOSe 

Pile Grade A COz — — 4 x 10-3 

Boulangier 2) H2O 2b SCA O= 10 4x 10-3 

Strickland—Constable #) H20 Wise x One 10 ~~ 8 X 10° 

Blyholder °) H20 Dexa Ome 100 ie 10m! 

Pile Grade A H20 — — 2 x 10-2 
Duval 2) | O32 0.77 10 
Pile Grade A Oz > 0.06 5-10 
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The greater dependence on the partial 
pressure of hydrogen is consistent with the 
adsorption of hydrogen atoms on active sites 
and the replacement of reaction (1) by 


Cy+ H,0 ———— Co+3 He, 


where Cy is a reaction site possessing a hydrogen 
atom. It is therefore in accord with the results 
for water vapour alone and the work of 
Boulangier 8). 

It follows that eq. (4) may still be valid for 
COz-CO additions as the reaction rate appeared 
to be directly proportional to concentration 
over the whole range 0-500 vpm for carbon 
dioxide alone (fig. 2) and Bouiangier found the 
efficiency of collision (F') independent of the 
partial pressure of carbon dioxide. 

These out of pile results provide a sound basis 
on which to conduct in pile experiments and 
indicate the relationships between gas compo- 
sition and attack of graphite in a high temper- 
ature gas cooled reactor. They should be 
directly applicable to a zero energy assembly 
but can only be used tentatively for a power 
reactor, as irradiation may activate the im- 
purities. 


5. Conclusions 


The reactivity of the gaseous impurities which 
may be present in a high temperature gas cooled 
reactor decreases in the order oxygen, water 
vapour, carbon dioxide. The reaction rates 
increase with temperature up to 1000° C and 
with the partial pressure of the oxidising gas, 
whilst hydrogen and carbon monoxide inhibit 
the attack. 

The relationships between the reaction rates 
and the partial pressures of water vapour, 
carbon dioxide, hydrogen and carbon monoxide 
are in good agreement with theory, previous work 
and the supposition that hydrogen atoms are 
adsorbed on graphite in contact with hydrogen 
and water vapour. 
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The characteristic temperature for PuOz has been 
determined by experimentally measuring the reduction 
in intensity of X-ray diffraction spectra from a 
specimen as the specimen temperature is increased. 
In order to evaluate the technique employed, the 
characteristic temperature of aluminium was also 
determined. 

For aluminium, a characteristic temperature of 
393° K was obtained, which agrees very well with the 
value of 390° K obtained from specific heat measure- 
ments. The characteristic temperature of PuO2 was 
found to be 415° K. 


La température caractéristique de PuOzg a été déter- 
minée en mesurant expérimentalement l’abaissement 
dintensité du spectre de diffraction aux rayons X en 
fonction de Vlélévation de température pour un 
échantillon donné. 

Parallélement la température caractéristique de 
Valuminium a aussi été déterminée afin de pouvoir 


1. Introduction 


When discussing the thermal properties of 
crystal lattices it is convenient to use the 
Debye-Waller theory for correlating various 
thermal data. According to this theory the 
thermal properties of a crystal lattice are 
associated with a single parameter, i.e. a 
characteristic temperature, 9, which, when 
substituted into the appropriate formulas, 
allows the specific heat of different substances 
to be correlated with a common curve. 

This paper is concerned with the determi- 
nation of the characteristic temperature for 
PuOz through measuring the reduction in 
intensity of X-ray diffraction spectra from a 


apprécier les techniques employées. Par cette méthode, 
on obtient une température caractéristique de 393° K 
pour aluminium ce qui concorde bien avec la valeur 
de 390° K obtenue par la mesure de la chaleur spéci- 
fique. On a trouvé une température caractéristique de 
415° K pour PuOs. 


Es wurde die charakteristische Temperatur von PuOzg 
bestimmt, indem bei steigender Probentemperatur die 
Abnahme der Intensitat von Rontgenbeugungsbildern 
gemessen wurde. Um die angewendeten Unter- 
suchungsmethoden zu tberpriifen, wurde zusatzlich 
noch die charakteristische Temperatur von Aluminium 
ermittelt. 

Die fiir Aluminium festgestellte charakteristische 
Temperatur von 393° K stimmt sehr gut mit dem 
Wert von 390° K iiberein, welcher bei Messungen der 
spezifischen Warme gewonnen wurde. Die charakter- 
istische Temperatur von PuOsg ergab sich zu 415° K. 


specimen as the specimen temperature is in- 
creased. In order to evaluate the techniques 
employed, the characteristic temperature of 
aluminium was also determined. 


2. Theory 


In the following discussion let the temperature 
T, be greater than the temperature 72. Then it 
follows for B, a quantity related to the average 
lattice vibration, that B, is greater than Bg. 
The intensity, J, of an X-ray reflection is 
related to B in the following manner 


I ~ exp (—2B sin? 0/2?) (1) 
where @ is the Bragg angle of diffraction and 


+ Work done under the auspices of the U.S. Atomic Energy Commission. 
tft <A portion of this paper was presented at the 17th Pittsburgh Diffraction Conference, November, 1959. 
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A is the wavelength of the X-rays employed. 
By taking, for a given reflection, the ratio of 
the intensities at two different temperatures, 
eq. (2) is obtained: 


Io/I, =a exp (2(B,— Be) sin? 6/A?). (2) 


According to the Debye-Waller theory the 
quantity B is related to the characteristic 
temperature, 0, by 


B= 6h?T (d(x) + tax) /mk 02 (3) 
where 


h=Planck’s constant, 6.6252 x 10-2" erg sec, 
m=average mass of a lattice point, g, 
k=Boltzmann’s constant, 
1.3804 x 10-16 erg deg-1, 
T =absolute temperature, °K, 
2=O(T, and 
(x)= Debye function. 


Assuming that © does not change in the 
temperature interval 7’: to 7’; then from eq. (3) 
if follows that 


_ 6h2[T1(d4a1 +44) —T2(p4%2+4 42] 
Q2 = ae SEE ESE Re ca (4) 


Thus, the value of (B,—B:), obtained from 
experimentally determined intensity ratios 
according to eq. (2), when inserted in eq. (4) 
allows the characteristic temperature of a 
substance to be calculated. 

In the Debye-Waller theory it is assumed that 
the volume of the crystal remains constant as 
the temperature changes. In actual practice the 
pressure remains constant and the volume 
changes as a result of thermal expansion. 
Consequently, if 0 is plotted as a function of 
the average temperature of the temperature 
interval 72 to 71, a curve is obtained which 
exhibits negative curvature as the average 
temperature increases. In order to obtain a 
single parameter for the material, the true 
characteristic temperature, a correction must be 
applied to eliminate this high temperature 
curvature. 

This correction for high temperature curva- 
ture is made in the following manner. Again, 


let 71 be greater than 72, but also restrict 71 
and 7’, to values in the neighbourhood of, or 
greater than @. In this region the value of the 
function (¢(x)+42) of eq. (3) does not differ 
appreciably from 1.00 and, therefore, to the 
first approximation this function may be 
neglected. Then from eq. (3) 


By = 6)2T)/mkO? (5) 
and 
= 6h2T'o/mk O22. (6) 


Paskin !), using thermodynamic reasoning, has 
derived the expression 


O;/@2= (V2/ Vi)" (7) 


where V is the volume of the crystal at the 
appropriate temperature and y is the Griineisen 
constant, which is equal to V6/Cyk where V is 
volume; , the coefficient of volume expansion; 
Cy, specific heat at constant volume; and &, 
the compressibility. By substituting eq. (7) in 
eq. (5), eq. (8) is obtained. 


By = 


6h2 Ty ee 
mk O22 (7, 


(8) 


By subtracting eq. (6) from eq. (8) and re- 
arranging, eq. (9) is obtained. 


Px ea ict 


A comparison of eq. (4) (Debye) with eq. (9) 
(Paskin) shows that these two equations for 
computing a characteristic temperature of a 
material are similar. The Paskin equation, 
however, takes into consideration the thermal 
expansion of the crystal and thus provides for 
a more practical experimental determination of 
the characteristic temperature. 


3. Experimental Procedure 


The experimental samples consisted of pressed 
and sintered disks approximately 0.5 in. in dia. 
by 0.030 in. thick. The samples were placed in 
a tantalum holder which in turn was inserted 
into a high temperature, high vacuum furnace. 
mounted on a General Electric XRD-3 X-ray 
diffractometer. A more complete description of 


Ee 
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this furnace and vacuum system is available 
elsewhere 2). Cu Kx X-rays were used and the 
intensity data were obtained by means of a 
standard General Electric SPG No. 2 Detector 
utilizing an argon-filled proportional counter 
tube. Pulse height analysis was used to dis- 
criminate against unwanted radiation. 

Because of the normal day-to-day random 
variations in X-ray tube output, counter tube 
efficiency, electronic stability, etc., the intensity 
data gathered at various temperatures and 
inserted in eq. (2) may not always result in a 
value of 1.00 for the quantity a. Since a has 
no effect on the desired experimental quantity 
B,— Bz, its value may be adjusted as a compu- 
tational convenience. Consequently, the raw 
intensity data obtained at various temperatures 
have been normalized to the data obtained at 
the highest temperature reached by a given 
specimen. The normalized data are given in the 
Tables in the following section. 

The intensity ratios obtained at various 
temperatures were least-squares fitted to eq. (2) 
by utilizg an iterative procedure 3) on an 
IBM 704 computer. This procedure allows the 
quantity (B,— Bz) of eq. (2) and its standard 
deviation to be calculated simultaneously. 


4. Experimental Results 
4.1. ALUMINIUM 

The experimental data for aluminium are 
given in table 1. 


TABLE 1 


Relative integrated intensities for aluminium 
(Normalized to 318° C) 


hkl 18°C 136° C 212°C 318°C 
111 2967 2896 2815 2705 
200 1502 1451 1401 1334 
220 866 804 758 701 
311 908 808 747 667 
222 241 216 195 WAL 
400 72 62 55 46 
331 337 282 242 197 
420 379 304 262 211 
422 416 328 272 208 


Average limit of error on intensity measurements 


+ 2.0 %. 


— 


The values of (6; — Bg) given in table 2 were 
obtained by substituting intensity ratios deter- 
mined from the data of table 1 into eq. (2). 


TABLE 2 


Experimental values of (6:—Be2) for aluminium 


Temperagure interval (BeBe lo-tecdane 
(°C) 
18-136 0.3618 + 0.0169 
136-212 0.2463 + 0.0089 
212-318 0.3603 + 0.0103 


The characteristic temperature, 0, deter- 
mined from the values of (B,— Bz) of table 2, 
when substituted into the Debye (eq. (4)) or 
Paskin (eq. (9)) equations is given in table 3. 


TABLE 3 


Characteristic temperature, ©, for aluminium as 
determined from the Debye and Paskin equations 


cue erature Debye Temperature Paskin 
inter! (°K) (°C) (°K) 
(°C) 
18-136 367 + 18 18 393 + 20 
136-212 358 + 10 136 394 + 10 
212-318 351 + 11 212 393 + 13 


For use in Paskin’s equation, the volumes were 
determined by cubing the lattice constants. The 
lattice constant was determined as a function of 
temperature by noting the Bragg angle of the 
422 reflection while intensity measurements 
were being taken. The Griineisen constant for 
aluminium is 2.17, taken from the work of 
Slater 4). The value of 393° K thus obtained for 
the characteristic temperature of aluminium 
agrees very well with the value of 390° K 
obtained from specific heat measurements °). 


4.2. PLUTONIUM DIOXIDE 


The experimental intensity data for PuOz are 
given in table 4. 

The values of (6,— Bz) are given in table 5 
and were obtained by substituting into eq. (2) 
the intensity ratios determined from the data 
of table 4. 

The characteristic temperature, 0, determined 
from the values of (B,— Bz) of table 5 when 
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substituted into the Debye (eq. (4)) or Paskin 
(eq. (9)) equations, is given in table 6. 


TABLE 4 


Relative integrated intensities for PuOs 
(normalized to 617° C) 


hkl 16°C | 196°C | 342°C | 474°C | 617°C 
200 1048 1039 1030 1024 1013 
311 1102 1072 1050 1024 1000 
400 214 206 200 194 187 
420 — 330 — 306 292 
333 525 491 464 442 410 
531 — 498 — 437 408 
620 440 400 371 345 315 
622 470 425 389 358 323 


Average limit or error on intensity measurements 


Se heey Ap 
TABLE 5 
Experimental values of (B1—B2) for PuOz 


Temperature interval 
° (Bi—Bz2) x 10-16 em? 

(°C) 

16-196 0.1351 + 0.0039 
196-342 0.1148 + 0.0046 
343-474 0.1090 + 0.0025 
474-617 0.1302 + 0.0113 


For use in Paskin’s equation the volumes 
were determined by cubing the lattice constant. 
The lattice constant was determined as a 


TABLE 6 


Characteristic temperature, 9, for PuOz as determined 
from the Debye and Paskin equations 


7 naa Debye Temperature Paskin 
interval (°K) (°C) (° K) 

(°C) 

16-196 403 + 8 16 414+ 9 
196-342 397 +--11 196 416+ 11 
342-474 388 + 6 342 4146+ 7 
474-617 376 + 14 474 414+ 13 


function of temperature by noting the Bragg 
angle of the 622 reflection while intensity 
measurements were being taken. The Griineisen 
constant for PuOg is empirically 2.6, this value 
being required in order to remove the high 
temperature curvature of the Debye character- 
istic temperature. The characteristic temperature 
of PuOz was thus found to be 415° K. 
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Das Modell der korrelativen Stossprozesse eignet sich 
fur die Deutung des Strahlungswachstums von «-Uran 
und der damit in Verbindung stehenden Erschei- 
nungen. In der Umgebung eines durch Spaltung 
hervorgerufenen Storbereiches (Fission Spike) treten in 
den 3 kristallographischen Hauptrichtungen folgende 
Prozesse bevorzugt auf: 


[100] — Fokusonen; 
[010] — Dynamische Crowdions; 
[001] — Hohe Dichte an Frenkelpaaren. 


Der bei héheren Temperaturen beobachtete Uber- 
gang vom anisotropen Strahlungswachstum zur iso- 
tropen Strahlungsschwellung ergibt sich aus der 
Defokussierung der korrelativen Stossprozesse. 

Die Korngréssenabhangigkeit und die Méglichkeiten 
der Verminderung des Strahlungswachstums werden 
diskutiert. Auf die bisher entwickelten Modelle (Spike 
und Diffusion) wird hingewiesen. 


The correlated collision model is apt for the inter- 
pretation of the irradiation-induced growth of «- 
uranium, and phenomena associated therewith. In the 
neighbourhood of a fission spike the following processes 
occur preferentially in the three principal crystallo- 
graphic directions: 


[100] — focusons; 
[010] — dynamic crowdions; 
[001] — high concentration of Frenkel pairs. 


1. Einleitung 


Uran, der wichtigste Brennstoff der heutigen 
Reaktoren, erfahrt unter Bestrahlung viel- 
seitige Anderungen seiner mechanischen und 
physikalischen Higenschaften. Die auffallendste 
Erscheinung ist seine Forminstabilitat, welche 
sich in Langen-, Oberflachen- und Dichte- 
anderungen deutlich ausdriickt. Die unmittel- 
bare Bedeutung dieser Forminstabilitat fiir die 


43 


The transition from anisotropic growth to isotropic 
swelling under irradiation, as observed at elevated 
temperatures, is a consequence of the loss of focusing 
of the correlated collisions. 

The grain size dependence of irradiation-induced 
growth, and the possibility of reducing it, are discussed. 
The existing models of growth (spike and diffusion 
models) are considered. 


Le modéle des processus de choc corrélatifs se préte 
bien & l’explication de la croissance sous irradiation 
de uranium « et des phénoménes qui s’y rattachent. 
Au voisinage d’un domaine perturbé préférentielle- 
ment par une fission (pointe de fission), les processus 
suivants prennent naissance préférentiellement suivant 
les 3 directions cristallographiques: 


[100] — Focusons; 
[010] — Crowdions dynamiques; 
[001] — Densité élevée de paires de Frenkel. 


Le passage observé, aux plus hautes températures 
de la croissance sous irradiation, au gonflement 
isotrope sous irradiation provient de la défocalisation 
des processus de choe corrélatifs. 

La relation entre croissance et dimensions des 
erains ainsi que les possibilités de diminuer la crois- 
sance sous irradiation sont discutées, ainsi que les 
modéles développés jusquwici (pointe de fission, 
diffusion). 


Reaktorentwicklung und -technik waren in den 
letzten Jahren Anlass fiir umfassende For- 
schungsarbeiten auf diesem Gebiete in allen 
beteiligten Lindern, besonders in den USA, in 
England und in Russland. Anlasslich der 
beiden Genfer Konferenzen zur friedlichen Aus- 
nutzung der Atomenergie (1955 und 1958) 
sind die Ergebnisse zusammenfassend vorge- 
tragen worden. 
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2. Forminstabilitat des a«-Urans 


Bei Bestrahlung von Uranproben lassen sich 
3 Erscheinungen der Forminstabilitat unter- 
scheiden : 


Beds 


Turkalo!) und Paine und Kittel #) zeigten, 
dass Kinkristalle von x-Uran eine Ausdehnung 
in [010]-Richtung und eine Schrumpfung in 
[100]-Richtung erfahren. In [001]-Richtung 
tritt kaum eine Langenanderung ein. In poly- 
kristallinem Uran entspricht das Strahlungs- 
wachstum der vorhandenen Textur. Fiir kleine 
Langenanderungen ist das Strahlungswachstum 
proportional der ‘“Verbrennung” (burn-up). Als 
sein Mass wird eine Wachstumsgrésse G; 


definiert 
[=] [010] 


Spaltprozesse 
Anzahl der Uranatome 


STRAHLUNGSWACHSTUM (GROWTH) 


Gy = 


Sie hangt sehr von der Textur und der Tempera- 
tur bei der Bestrahlung ab. Die Ergebnisse von 
Paine und Kittel 2) fiir «-Uraneinkristalle sind 
in Tabelle 1 wiedergegeben. 


TABELLE | 
(nach Paine und Kittel 2)) 


Strahlungswachs- Thermischer 
Richtung | tumskoeffizient | Ausdehnungskoeffizient 
Gi a (25-125° C) 
[100] — 420 + 20 se allnd S< OSE 
[010] + 420 + 20 Sullyay <= 
f001j 0 + 20 + 23,2 x 10-6 
2.2. STRAHLUNGSSCHWELLUNG (SWELLING) 


Bei héheren Bestrahlungstemperaturen (ober- 
halb 400° C) tritt an die Stelle des anisotropen 
Strahlungswachstums eine isotrope Ausdehnung, 
die nicht nur auf die «-Phase des Urans 
beschrankt bleibt. Sie ist weitgehend unab- 
hangig von der Vorgeschichte der Proben. Unter 
dem Mikroskop erkennt man eine Porenbildung. 
Der Gitterparameter bleibt jedoch annahernd 
unverandert. 


2.3. STRAHLUNGSRUNZELN (WRINKLING) 


An Oberflichen grobkristalliner, urspriinglich 
glatter Proben zeigen sich nach Bestrahlung 
Runzeln. Ihre Bildung hangt deutlich erkennbar 
mit der Kornstruktur der Proben zusammen. 
Die Runzeln wird man daher als eine Folge- 
erscheinung des Strahlungswachstums ansehen 
kénnen. Einige Kérner rufen durch ihre bevor- 
zugte Lage an der Oberflache starkere Wachs- 
tumserscheinungen hervor als andere, die in 
ihrem Wachstum gehemmt sind. 


3. Korrelative Stossprozesse 
(Correlated Collisions) 

Die experimentellen und theoretischen Ar- 
beiten iiber die Strahlenschadigung (Radiation 
Damage) in Festkérpern fiihrten zu der Vor- 
stellung, dass ein Atom aus seinem Gitterplatz 
herausgestossen wird, wenn bei dem Stoss 
eine Mindestenergie (Threshold Displacement 
Energy; Ha 25 eV) iibertragen wird. Erhalt 
das primar herausgeschlagene Atom (Primary 
Knock-on) bei dem Stossprozess eine hohe 
Energie, so kénnen kaskadenartig sekundare, 
tertiare usw. Atome aus dem Gitterverband 
gelost werden. Im Energiebereich von etwa 
104 eV liegt nach den Arbeiten von Brink- 
man *) fir viele Materialien die Stossfolge 
in der Gréssenordnung der atomaren Ab- 
stande. Die resultierenden Bereiche hoher 
Unordnung (Displacement Spike) sind mit dem 
Schmelzzustand verglichen worden. Silsbee 4) 
wies zuerst darauf hin, dass durch die Anordnung 
der Atome im Gitter die Stossprozesse in korrela- 
tiver Weise in bestimmten Richtungen erfolgen. 
Nach seinen Vorstellungen wird deutlich, dass 
nicht allein die mittlere freie Weglinge der 
primar herausgestossenen Atome den Schadi- 
gungsbereich charakterisiert, sondern dass kor- 
relative Stossprozesse die Storung aus dem 
Zentralbereich weit in das Gitter hineintragen. 
Thompson 5) zeigte, dass bei Protonen-Bestrah- 
lung von Goldfolien eine Emission von Gold- 
atomen in Richtung dicht-gepackter Reihen 
erfolgt. Dieses Experiment muss als eine Verifi- 
zierung des Modells von Silsbee angesehen 
werden. Durch neuere Arbeiten von Leibfried 6) 
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und Vineyard ’?) werden diese Ideen erweitert 
und verfeinert. 

Nach Leibfried sind bei den weitreichenden 
Effekten der korrelativen Stossprozesse zwei 
Moglichkeiten fiir das Gitter zu unterscheiden: 


1) Fokussierende Stossprozesse (Focusing 
Collisions ) 


Sie werden dadurch charakterisiert, dass 
lediglich Energie iibertragen wird, aber kein 
Transport von Zwischengitteratomen erfolgt. 
Fiir diese Art von Stossprozessen schlagt Leib- 
fried den Namen ‘‘Fokusonen”’ (Focusons) vor. 


2) Crowdion-Stossprozesse (Crowdion 
Collisions ) 


Ein Crowdion entsteht dadurch, dass eine 
Atomreihe auf kurzer Strecke ein Atom mehr 
aufnimmt, als der idealen Gitterbesetzung ent- 
spricht. Als ““dynamisches Crowdion”’ wandert 
diese Struktur lings ihrer Gittergeraden. Als 
“statisches Crowdion’ ist diese Struktur 
wahrscheinlich nicht bestandig; vielmehr geht 
das damit gekoppelte anisotrope Verschie- 
bungsfeld in das Kubisch-symmetrische eines 
Zwischengitteratoms tiber. Das Crowdion stellt 
demnach einen kurzlebigen Gitterfehler dar. 

Bei der Behandlung der korrelativen Stoss- 
prozesse wird das steile Potential zwischen zwei 
Atomen im Gitter durch ein hartes Kugel- 
Potential ersetzt, wobei der Radius der als starr 
zu betrachtenden Kugeln (4) mit abnehmender 
kinetischer Energie # des bewegten Atoms 
anwachst. 

Fokussierende Stdsse sind nach Abb. 1 + dann 
zu erwarten, wenn Rk > 4Dr ist. Dp ist der 
Abstand zweier Atome in Richtung dicht- 
gepackter Atomreihen. Der Stosspunkt Pr liegt 
naher an der Gleichgewichtslage des stossenden 
Atoms. Dadurch wird bei jedem Stoss der 
Winkel zur dicht-gepackten Atomreihe kleiner, 


+ Es soll darauf hingewiesen werden, dass in der 
Darstellung die Atome mit doppeltem harten Kugel- 
radius gezeichnet wurden. Dadurch werden die 
Stossverhaltnisse anschaulich, besonders die Unter- 
scheidung von Fokuson und Crowdion. 


Abb. 1. a-b-Ebene des «-Urans mit Darstellung von 


Crowdion Stoss und Fokuson. 


und es tritt eine Fokussierung in diese Richtung 
ein. Da R eine Funktion von £ ist, sollte sich 
nach dieser Vorstellung eine maximale Energie 
Ey fiir das Einsetzen der Fokussierungs-Stdésse 
angeben lassen. In diesem Falle ware R=4D, 
so dass sich die Kugeln in der Darstellung der 
Abb. 1 gerade berithren. Durch Fokusonen kann 
eine erhebliche Energieiibertragung iiber 100 
Atomabstande eintreten. Im Kristallgitter wird 
jedoch die Energie durch Anregung thermischer 
Schwingungen dissipiert, ohne dass Gitterfehl- 
stellen gebildet werden. Werden jedoch gewisse 
schon vorhandene Gitterfehlstellen (Versetzung, 
Korngrenze, Zwillungsgrenze usw.) von einem 
Fokuson getroffen, so kann das letzte Atom, 
das keinen direkten Stosspartner dieser Reihe 
mehr vorfindet, aus dem Gitterverband auf 


einen Zwischengitterplatz herausgestossen 
werden, sofern nur auf dieses Atom eine aus- 
reichende Schwellenenergie (Threshold Dis- 


placement Energy) tibertragen wird. 
Crowdion-Stossprozesse sind dann zu _ er- 
warten, wenn der Stosspunkt Pc naher an der 
Gleichgewichtslage des angestossenen Atoms 
dieser Gitterreihe liegt. 
In Abb. 1 sind auch die Verhaltnisse fiir den 
Fall des Crowdions dargestellt. Im Gegensatz 
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zum Fokuson tritt in diesem Fall ein Platz- 
wechsel der Atome ein. Vernachlassigt man die 
Abhangigkeit des Radius R von der Energie H 
des stossenden Atoms, so wird der Winkel bei 
jedem Stoss grésser, bis das crowdionbildende 
Zwischengitteratom aus der Kristallreihe hinaus- 
gestossen wird oder mangels Energie in dieser 
Atomreihe als Zwischengitteratom “hangen”’ 
bleibt. 


4. Korrelative Stossprozesse im a-Uran 


Die Anwendung des Silsbee—Leibfried-Modells 
auf das «-Urangitter kann zu einer Deutung des 
Strahlungswachstums und der damit in Ver- 
bindung stehenden Erscheinungen dienen. 

Zunachst sei naher auf die bis 668°C 
bestaindige orthorhombische Struktur des «a- 
Urans eingegangen. In der Einheitszelle treten 
dabei folgende Atomabstinde auf, wie aus Abb. 2 
entnommen werden kann: 


2,762 A zweimal (A-B) 
2,852 A zweimal (C-D) 
3,261 A viermal (B-C) 
3,322 A viermal ) 


sae hae he Ae 


Man erkennt ferner aus der Abbildung, dass in 
der c-Richtung, [001] keine dicht-gepackte 
Reihe von Atomen existiert, wenngleich kleine 
Atomabstainde vorkommen. Dagegen liegen in 
der a—b-Ebene dicht-gepackte Atomreihen mit 
Atomabstainden von 2,852 A und 3,261 A vor. 
Die Struktur dieser Gitterebene ist in Abb. 1 


= 2,8524—- 
Abb. 2. 


Die Struktur des «-Urans (nach Wilkinson 
und Murphy 8)). 


GONSER 


bei der Abhandlung der korrelativen Stoss- 
prozesse oben schon verwandt worden. Die 
Struktur des «-Urans kann als verzerrtes 
hexagonales Gitter aufgefasst werden. 

Bei der Spaltung (Fission) eines Uranatoms 
werden etwa 10 Spaltungsprodukte (Spalt- 
fragmente, Neutronen, $-Strahlung, y-Strahlung 
und Neutrinos) gebildet unter Freisetzung von 
mehr als 108 Elektronen-Volt. Der grésste Teil 
dieser Energie entfallt auf die Spaltfragmente, 
die als Ursache fiir den Strahlungsschaden 
angesehen werden miissen. Der Weg dieser 
Fragmente ist durch Bereiche hoher Storung 
ausgezeichnet, die allgemein Spikes genannt 
werden (Thermal Spike, Displacement Spike). 
Die Theorie der Spikes ist besonders von 
Brinkman 9) ausgearbeitet worden. Danach 
koénnen in einem Displacement-Spike des «a- 
Urangitters maximal etwa 60 000 eV umgesetzt 
werden. Nach einem zugrunde gelegten Bild 
werden aus einem Zentralbereich (verdiinnte 
Zone) Atome mit relativ hoher Energie in eine 
verdichtete Zone hineingestossen. Die Zu- 
sammenstésse mit den Atomen der Matrix sind 
dabei zunachst voneinander unabhangig (R ist 
klein). St6sse energiearmerer Atome (~ 102 eV) 
in Richtung dicht-gepackter Atomreihen werden 
iiber korrelative Stossprozesse als dynamische 
Crowdions oder Fokusonen iiber viele Atom- 
abstainde Energie Gitter hinein  iiber- 
tragen. 

Die verschiedenen Gitter-Richtungen im «- 
Uran sollen nun einzeln diskutiert werden. 
Besonders interessieren die [100]-Richtung, in 
der eine Schrumpfung unter Bestrahlung statt- 
findet und die [010]-Richtung, in der Strah- 
lungswachstum eintritt. 

[100]-Richtung: In dieser Richtung liegt 
eine Atomreihe mit Abstainden von nur 2,852 A 
vor. Durch die ausserst dichte Besetzung dieser 
Reihe ist diese Richtung fiir Fokusonen bevor- 
zugt. Im idealen Gitter kénnen diese keine 
Fehlstellen hinterlassen. Treffen sie jedoch auf 
eine Korngrenze, so kann die Energie aus- 
reichen, um das letzte Atom der Reihe in das 
benachbarte Korn hineinzustossen, wo es zu- 
nachst einen Zwischengitterplatz einnimmt. 


ins 
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Einige dieser Zwischengitteratome werden bei 
erhéhter Temperatur durch Zwischengitter- 
diffusion der Korngrenze entfliehen und sich 
an einer anderen Korngrenze oder Versetzung 
ablagern. Die betrachtete Atomreihe hat also 
bei diesem Prozess ein Atom verloren. 

[010]-Richtung: Diese Richtung ist durch 
das Strahlungswachstum ausgezeichnet. Aus 
Abb. 2 erkennt man, dass zwei dicht-gepackte 
Atomreihen [110] existieren, die mit der [010]- 
Richtung einen Winkel von etwa 26° ein- 
schliessen. Bei Atomabstanden in diesen Reihen 
von 8,261 A ist die Besetzung also aufgelockerter. 
Daher sind in diesen Richtungen Crowdions 
bevorzugt. Die Druckwelle der Verdichtungs- 
zone um einen Spike lasst Atome tiber dyna- 
mische Crowdions weit in das Gitter hinein- 
stossen. Die benachbarten Reihen wirken der 
Defokussierung der Crowdion-Stossprozesse ent- 
gegen. Der damit verbundene Atomtransport 
endet im idealen Gitter mit emem Zwischen- 
gitteratom. Trifft ein energiereiches Crowdion 
auf eine Korngrenze, so ist es in der Lage, 
das letzte Atom als Zwischengitteratom in 
das Nachbarkorn zu stossen. Die meisten 
Crowdions bleiben jedoch bevorzugt an der 
Korngrenze hangen und lagern das crowdion- 
bildende Zwischengitteratom an der Korngrenze 
ab, was zu einem Wachstum in dieser Richtung 
fiihrt. 

Fiir die beiden Atomreihen der a—b-Ebene des 
Urans kann die Schwellenenergie fiir Fokusonen, 
Ex, wie folgt abgeschatzt werden: Das Wechsel- 
wirkungspotential von Uranatomen soll durch 
die folgende Form 1°) gegeben sein 


V(R) = Ae-Rla (1) 


mit den Konstanten A=8x104 eV und a= 
2,1x10-9 em. Wird der Radius entsprechend 
der Beriithrungsbedingung (R=4D) eingesetzt, 
so erhalt man mit der Beziehung V(R)=4# °) 


Ey[100]=180 eV 
Ey{110]= 67 eV. 


Danach k6nnen herausgestossene Atome mit 
einer Energie von weniger als 180 eV _ bei 


ginstigem Auftreffen in  [100]-Richtungen 
Fokusonen entwickeln. Entsprechendes gilt fiir 
die Atome mit 67 eV in [110]-Richtungen. Die 
Verschiedenheit der Schwellenenergien Hy fiir 
die Atomreihen ist das Merkmal dieses Wachs- 
tumsmodells. Energiereiche Atome im Bereich 
von 180-67 eV werden in [100]-Richtungen 
vorwiegend Fokusonen bilden, jedoch nur 
Crowdions in [110]-Richtungen. 

[001 ]-Richtung: Im «-Uran sind keine dicht- 
gepackten Atomreihen in dieser Richtung 
vorhanden. Die energiereichen Stossteilchen 
(Knock-ons) der Verdichtungszone der Spikes 
konnen ihre Energie nicht mittels weitreichender 
korrelativer Stossprozesse abgeben. Daher 
werden sie nur in engumschriebenen Bereichen 
Atome aus ihren Gitterplaitzen herausstossen. 
Dies fiihrt schliesslich zu einer héheren Dichte 
von Zwischengitteratomen in der [001]-Rich- 
tung. Diese Verdichtungszone kann spater unter 
Umstanden in die entstandene “‘Pore’’ (cluster) 
zuriickdiffundieren. 


5. Das Modell der korrelativen Stossprozesse 
zur Erklarung verschiedener Erscheinungs- 
arten der Forminstabilitat des «-Urans 


Die Forminstabilitat des «-Urans zeigt sich 
phanomenologisch auf verschiedene Weise. Ihre 
wichtigsten makroskopischen LErscheinungen 
lassen sich auf der Grundlage korrelativer 
Stossprozesse deuten. 


5.1. RIcHTUNGSABHANGIGKEIT DES 
STRAHLUNGSWACHSTUMS 


Nach dem im letzten Abschnitt dargestellten 
qualitativen Modell treten in den verschiedenen 
Richtungen folgende Prozesse bevorzugt auf: 


a) [100] Fokusonen — Energietransport — An 
der Korngrenze werden Atome in das 
Nachbarkorn gestossen. 

b) [010] Crowdions — Transport von Zwischen- 
gitteratomen — Atome werden an der Korn- 
grenze abgelagert. 

c) [001] Hohere Dichte an Zwischengitter- 
atomen und dadurch Abriegelung fiir weit- 
reichende Stossprozesse. 
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5.2. ABSCHATZUNG DES STRAHLUNGSWACHSTUMS 


Nach den Uberlegungen von Leibfried ist 
die Zahl der gebildeten Crowdions mit einer 
Reichweite von n Atomabstinden im kubisch- 
flachenzentrierten Gitter gegeben durch 

dal : 

C= 32F tae ie (2) 
F ist die Zahl der aus dem Gitterverband 
herausgestossenen Atome, die fiir einen Spalt- 
prozess von Brooks 1!) mit 325 000 angegeben 
wird. Nimmt man an, dass fiir das Urangitter 
mit 12 nahen Nachbarn (Verwandschaft zum 
hexagonalen System!) die Zahl der Crowdions in 
ahnlicher Weise abgeschatzt werden kann und 
dass davon etwa 1/3 in die makroskopische 
Wachstumsrichtung stossen, so ergibt sich mit 
einer mittleren Reichweite von 30 Atom- 
anstainden fiir einen Spaltprozess 


Cio10) 1000. 


Kin Strahlungswachstumskoeffizient, G;, von 
einigen Hundert, wie in der Praxis beobachtet, 
ware auf diese Weise verstandlich. 


5.3. STRAHLUNGSWACHSTUM ALS FUNKTION 
DER KORNGROSSE 


Nach SBeobachtungen verschiedener For- 
scher 2-12.18) jist das Strahlungswachstum in 
polykristallinen Kristallen mit ausgepragter 
Textur wesentlich starker ausgebildet als in 
grobkristallinen Materialien oder Kinkristallen. 
Diese Experimente lassen deutlich erkennen, 
dass die Korngrenzen eine bedeutende Rolle 
beim Wachstum spielen. Nach dem Modell der 
korrelativen Stossprozesse ist dieser Effekt auch 
zu erwarten. Die Ablagerung des crowdion- 
bildenden Zwischengitteratoms an Korngrenzen, 
Subkorngrenzen oder Zwillingsgrenzen ware als 
primare Stufe des Wachstums anzusehen. 
Danach wiirden alle gebildeten Crowdions aus 
der Randschicht mit der Breite von n Atomen 
eines Korns zum Wachstum beitragen. Das 
Strahlungswachstum sollte damit etwa um- 
gekehrt proportional zum Kornradius sein. 
Gentigend genaue experimentelle Ergebnisse 
liegen leider noch nicht vor, um diese Vermutung 
zu bestatigen. 


Strahlungswachs nef 


pee a4 


%Dichteanderung. 


Strahlungswachstum Gj 


%.Verbrennung 
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Temperatur C 


Abb. 3. Temperaturabhéngigkeit von Strahlungs- 
wachstum und Strahlungsschwellung (nach Seigle und 
Castleman !”)). 


5.4. TEMPERATURABHANGIGKEIT DES STRAH- 
LUNGSWACHSTUMS 


In Abb. 3 ist die Temperaturabhangigkeit des 
Strahlungswachstums dargestellt. Bei etwa 
300° C wird ein Maximum beobachtet, wahrend 
bei tiefen Temperaturen (< 0° C) und hohen 
Temperaturen (> 500°C) kaum noch Strah- 
lungswachstum stattfindet. 

Diese Erscheinungen kénnen wie folgt ge- 
deutet werden: Nach den_ beschriebenen 
Vorstellungen verlieren die Korner in der 
[100]-Richtung Atome, wahrend neue in der 
[010]-Richtung angebaut werden. An den Korn- 
grenzen entsteht ein Spannungsfeld. Der an- 
steigende Ast entspricht dann der thermisch 
aktivierten Relaxation der Korngrenzenbe- 
wegung. 

Auffallend ist der Kurvenbereich oberhalb 
von 400°C, wo das anisotrope Strahlungs- 
wachstum von einer anderen Erscheinung — 
Strahlungsschwellung —abgelést wird. Bei an- 
steigenden Temperaturen tritt durch die ther- 
mische Bewegung der Atome eine Defokus- 
sierung der korrelativen Stossprozesse in den 
dicht-gepackten Atomreihen ein. Die Reich- 
weite der Fokusonen und Crowdions (n) nimmt 
ab. Dadurch werden Strahlungseffekte mehr auf 
die nahere Umgebung der Spikes verlagert. 
Die anisotrope Forminstabilitaét nimmt einen 
isotropen Charakter an. 


5.5. 
Bei héheren Temperaturen kann die Energie 


STRAHLUNGSSCHWELLUNG (SWELLING) 
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nicht mehr iiber weitreichende Stossprozesse in 
bevorzugten Richtungen in das Gitter hinein- 
getragen werden. Dadurch wird die Energie- 
dichte in den Storbereichen grésser. Die Stor- 
bereiche nehmen kugelformige Gestalt an. Gefor- 
dert wird dieser Vorgang dadurch, dass das 
Gitter bei héheren Temperaturen der dicht- 
gepackten hexagonalen Struktur ahnlicher wird 
(siehe Ausdehnungskoeffizienten verschiedener 
Gitterrichtungen in Tab. 1). Nach diesen Vor- 
stellungen sollte das anisotrope Strahlungs- 
wachstum bei Temperaturen, wo die weit- 
reichenden korrelativen Stossprozesse  aus- 
bleiben, durch einen isotropen Wachstums- 
mechanismus abgeldst werden, der als nicht 
mehr strukturabhangige Strahlungsschwellung 
beobachtet wird. Die elektronenmikroskopisch 
beobachteten Poren kénnten ihren Ursprung in 
solechen Stérbereichen haben. Das Anwachsen 
dieser Poren ist nach Greenwood, Foreman und 
Rimmer !4) durch einen thermisch aktivierten 
Diffusions-Mechanismus zu deuten. Besonders 
der Abbruch des Strahlungswachstums bei 
hoéheren Temperaturen und der gleichzeitige 
Ubergang zur Strahlungsschwellung kann als 
Stiitze dieses Modells angesehen werden. 


5.6. STRAHLUNGSWACHSTUM VON URAN- 


LEGIERUNGEN 


Das Strahlungswachstum ist vom technischen 
Standpunkt aus eine sehr unerwiinschte Higen- 
schaft der Brennelemente. Man kann versuchen, 
das Wachstum durch Legierungsbildung zu 
vermindern. Gemass dem Modell der korrela- 
tiven Stossprozesse kénnen zwei ‘‘Rezepte”’ 
angegeben werden, um das Strahlungswachstum 
zu reduzieren : 


a) Zusatze, die zu einer Herabsetzung der 
Anisotropie fiihren, bewirken, dass die auf- 
tretenden korrelativen Stossprozesse mehr iso- 
tropen Charakter annehmen und die Aus- 
zeichnung einer Wachstumsrichtung weniger in 
Erscheinung tritt. So bemerkten Konobeevsky 
und Mitarbeiter 15), dass durch Zugabe von 
Molybdin und Plutonium die Unterschiede 
zwischenatomarer Abstiande kleiner werden, das 
Gitter also dem idealen hexagonalen System 


zustrebt. Kittel und Paine 1*) beobachten, dass 
das Strahlungswachstum von Uran mit diesen 
Legierungskomponenten tatsaichlich merklich 
reduziert wird. 


b) Durch Zusiatze leichter Atome tritt bei 
den korrelativen Stossprozessen in einer dicht- 
gepackten Reihe ein erheblicher Energieverlust 
ein, wie sich unmittelbar aus der Stossformel fiir 
die maximal iibertragbare Energie ergibt 


4M,Me 


oe? (Mi +M2)? (3) 


M, und Mz sind die Massen der beiden Stoss- 
partner. Die Reichweite der korrelativen Stoss- 
prozesse wird kleiner und die Anisotropie damit 
herabgedriickt. 


6. Vergleich der korrelativen Stossprozesse mit 
anderen vorgeschlagenen Mechanismen zur 
Deutung der Forminstabilitat 


Kine Reihe von Mechanismen zur Deutung 
der Forminstabilitat des o«-Urans sind vorge- 
schlagen worden 17-18), Sie konnen auf zwei 
Grundvorstellungen reduziert werden: Spikes 
und Diffusion. 


6.1. SPIKES 


Diese Hypothese geht von der Vorstellung 
aus, dass die umgebende Matrix eines Spikes 
(Thermal Spike, Fission Spike) wahrend der 
kurzen Aufheizung und Abkiihlung plastische 
irreversible Verformung erfahrt. Fiir die Aniso- 
tropie konnen nun entweder die Verschiedenheit 
der Ausdehnungskoeffizienten oder die Unter- 
schiede in den plastischen Eigenschaften in den 
3 Richtungen herangezogen werden. Die Span- 
nungsrelaxation der Spike-Mechanismen kann 
durch Gleitung (Shp) oder durch Zwillungs- 
bildung (Twinning) erfolgen. Gegen einen Spike- 
Mechanismus sprechen folgende Argumente: 
a) Die raumliche Ausdehnung und die Lebens- 


dauer eines Spikes sind zu_ kurz, 
plastische Verformung hervorzurufen. 


um 


b) Die Korngréssenabhangigkeit ist schwerlich 
mit einem Spike-Mechanismus zu erklaren. 
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6.2. Dirrusion 


Seigle und Opinsky 19) schlugen zur Erklarung 
des Strahlungswachstums vor, dass die unter 
Strahlung gebildeten Zwischengitteratome und 
Leerstellen (Frenkelpaare) anisotrop diffun- 
dieren. Nach dieser Theorie sollten Zwischen- 
gitteratome besonders in [010]-Richtung und 
Leerstellen in [100]-Richtung wandern, wahrend 
die Diffusionsgeschwindigkeit in [001 ]-Richtung 
ungefahr fiir beide Fehlstellen gleich sein sollte. 
An den Korngrenzen bleibt die betreffende 
Fehlstelle hangen, was zu einer entsprechenden 
Ausdehnung oder Verkiirzung des Kristalls 
fiihrt. 

Die Anisotropie der Wanderung von 
Zwischengitteratomen und Leerstellen haben 
die Diffusionstheorie und das hier vorge- 
schlagene Modell zur Erklarung des Strahlungs- 
wachstums des o«-Urans in gewissem Sinne 
gemeinsam, jedoch liegt bei der Diffusion ein 
statistischer und thermisch aktivierter Prozess 
vor, waihrend die Energie fiir die korrelativen 
Stossprozesse aus der Wechselwirkung von 
Strahlung und Materie und deren Folge- 
erscheinungen entnommen wird. 

Experimentell sollte eine Verifizierung der 
korrelativen Stossprozesse dadurch méglich sein, 
dass die Emission von Uran-Atomen einer unter 
Spaltung stehenden LEinkristallprobe in den 
verschiedenen Richtungen untersucht wird. In 
Richtungen mit ausgepragter Langenanderung 
zeigen sich unter Umstianden ‘Muster’, die 
denen von Thompson>) entsprechen. Bei 
hoheren Temperaturen sollte beim Ubergang 


von Strahlungswachstum — Strahlungsschwel- 
lung (# 400°C) die anisotrope Emission zur 
einer isotropen werden. 
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Ageing experiments at 500°, 600° and 700°C on 
thorium-rich thorium-aluminium and thorium-urani- 
um alloys are described, together with the results of 
similar tests at 600° C on five more complex thorium 
alloys; hot hardness tests at temperatures up to 
800° C have been made on three samples of thorium 
and on several Th-Al, Th-U and Th-Zr alloys in the 
as-cast, solution treated and aged conditions. 

Aluminium additions improved the hot hardness of 
thorium and the alloys showed considerable hardening 
due to precipitation of TheAl. The rate of ageing 
increased directly with aluminium content and with 
temperature while the maximum hardness for a given 
composition increased inversely with temperature. 
Little precipitation hardening was observed in thorium- 
uranium alloys containing up to 3 at % uranium and 
their hot hardness behaviour showed little significant 
improvement over unalloyed thorium. The results on 
thorium-zirconium alloys suggest that zirconium 
additions within the «-thorium solid solution range 
lead to solution hardening. In an alloy containing 
22 at % zirconium, higher hot hardness values, up to 
600° C, were obtained after solution annealing in the 
B (body-centred cubic) phase-field than after solution 
treatment in the (« + #) field. 

Additions of aluminium and uranium in excess of 
the solubility limit restrict the grain growth of 
thorium at 1150° and 1200° C respectively. 


Des expériences de vieillissement a 500°, 600° et 700° C 
sur des alliages Th-Al riches en thorium et des alliages 
Th-U sont décrites, en méme temps que sont donnés 
les résultats d’essais analogues & 600° C sur 5 alliages 
de thorium plus complexes. Des essais de dureté a 
chaud & des températures jusqu’a 700°C ont été 
effectués sur trois échantillons de thorium et sur 
plusieurs alliages Th-Al, Th-U, et Th-Zr dans l'état 
coulé, homogénéisé en solution solide et vieilli. 

Des additions d’aluminium améliorérent la dureté 
& chaud du thorium et les alliages montraient un 


durcissement considérable di a la précipitation de 
TheAl. La vitesse de vieillissement croissait propor- 
tionnellement a la teneur en aluminium et avec la 
température tandis que la dureté maximum pour une 
composition donnée croissait de fagon inversement 
proportionelle a la température. Un faible durcissement 
par précipitation observé dans des alliages Th-U 
contenant jusqu’a 3 at % U et leur comportement 
a la dureté a chaud ne montrait pas une nette améli- 
oration vis a vis du thorium non allié. Les résultats 
sur les alliages Th-Zr suggérent que les additions de 
zirconium dans le domaine de la solution solide du 
thorium-« conduisent & un durcissement structural. 
Dans un alliage contenant 22 at % Zr, des valeurs de 
dureté a chaud plus élevées (jusqu’a 600° C) furent 
obtenus aprés recuit de mise en solution dans le 
domaine f (cubique centré) qu’aprés traitement de mise 
en solution dans le domaine (« + 8). 

Des additions d’aluminium et d’uranium au-dela de 
la limite de solubilité réduisent le grossissement de 
grain du thorium & 1150° et 1200° C respectivement. 


Es wurden Aushartungsversuche an thoriumreichen 
Thorium-Aluminium- und Thorium-Uran-Legierungen 
bei 500, 600 und 700° sowie an 5 komplizierter auf- 
gebauten Legierungen bei 600° vorgenommen. Warm- 
harte-Prifungen wurden bei Temperaturen bis 800° 
an drei Thorium-Proben und an_ verschiedenen 
Thorium-Aluminium, Thorium-Uran und Thorium- 
Zirkon-Legierungen durchgefiihrt. Die Proben be- 
fanden sich dabei in gegossenem, gegliihtem oder 
ausgehartetem Zustand. 

Aluminium-Zusatze verbesserten die Warmhirte 
von Thorium. Diese Legierungen zeigten eine betracht- 
liche Hartezunahme infolge Ausscheidung von TheAl. 
Die Aushartungsgeschwindigkeit stieg direkt mit dem 
Aluminiumgehalt und mit der Temperatur an, 
wahrend die maximal erzielbaren Hartewerte mit 
steigender Temperatur abnahmen. Bei Thorium- 
Uran-Legierungen mit Gehalten bis zu 3 At % Uran 
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wurde nur eine geringfiigige Héirtesteigerung durch 
Ausscheidung beobachtet. Auch die Warmharte 
dieser Legierungen lag nur unwesentlich tiber der von 
unlegiertem Thorium. Die an Thorium-Zirkon-Legie- 
rungen gewonnenen Ergebnisse lassen vermuten, dass 
Zirkonzusatze innerhalb des Bereichs von «-Thorium 
zu einer Hartesteigerung infolge Mischkristallbildung 
fiihren. Bei einer Legierung mit 22 At % Zirkon 


1. Introduction 


The solubilities of aluminium and of uranium 
in thorium have previously been shown to be 
temperature dependent !: 2), and ageing experi- 
ments on these binary alloys are now described 
together with the results of similar tests on four 
ternary and one quaternary thorium-base alloy. 
The hot hardness testing of metals and alloys 
is a useful, but not unambiguous, method of 
assessing their strengths at elevated tempera- 
tures and their deformation characteristics 
during hot working; in addition, the results of 
such tests on an alloy may indicate a range of 
temperature in which good creep resistance 
might be expected. An evaluation of some of 
these properties has been made by hot hardness 
tests on three samples of thorium and several 
thorium-base alloys at temperatures up to 
800° C. 


WATER COOLED PLATE 
= 


wurden nach einer Gliihung im f-Bereich (kubisch 
raumzentriert) hdhere Warmhartewerte (bis zu 600°) 
gefunden, als nach einer Gliihbehandlung im « + f- 
Bereich. 

Uberschritten die Aluminium- und Uran-Gehalte 
die entsprechenden Léslichkeitsgrenzen, so wurde das 
Kornwachstum von Thorium bei 1150 und 1200° 
abgeschwacht. 


2. Materials and Specimen Preparation 


Two grades of thorium were employed: (a) 
production bar material having a total thorium 
content of 99.5-99.6 wt % and containing 
400-700 ppm carbon, 3000-4000 ppm oxygen 
and total metallic impurities of less than 300 
ppm; (b) electrolytic thorium with a carbon 
content of 100-150 ppm and oxygen of 2000-— 
3000 ppm, metallic impurities being again less 
than 300 ppm. High purity aluminium, mag- 
nesium-reduced uranium billet (with a carbon 
content of approximately 300 ppm), iodide 
titanium and iodide zirconium were used as the 
alloying materials. 

All the alloys were prepared by arc-melting 
in a zirconium-gettered argon atmosphere; 
specimens for ageing tests were cast into a 9mm 
dia. mould while those for hot hardness testing 
were cast into cylinders 12 mm dia x 12 mm 
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long or “buttons” 25 mm in diameter. For 
solution treatment they were wrapped in 
tantalum foil, sealed in evacuated silica capsules, 
annealed at 900°-1200° C (depending on the 
material) and quenched by breaking the capsule 
under water. 


3. Ageing Tests 


For ageing treatments of up to 67 h the alloys 
were placed in a vertical, continuously evacuated 
silica tube and were heated by an external 
furnace which was interchangeable with a 
bucket of cold water for cooling. For longer 
annealing times, specimens were sealed in 
evacuated silica capsules, placed in a furnace 
and cooled by removal and immersion in a 
bucket of cold water without breaking the silica. 

The hardness of the aged materials was 
determined using a standard Vickers pyramid 
hardness tester with a 10 kg load; all results 
are the average of at least three values. 


4. Hot Hardness Test 


The hot hardness apparatus, shown in fig. 1, 
consisted essentially of a water cooled copper 
chamber evacuated by a rotary backing pump 
and an oil diffusion pump. The sample to be 
tested was sited on a platform in the bottom 
of the chamber and surrounded by a wire wound 
resistance furnace and tantalum reflector. 
Movement of the specimen was made by a 
stainless steel tube positioned in a small hole 
in the edge of the specimen and worked through 
a Wilson seal in the side of the chamber. A 
thermocouple was placed in the tube for direct 
reading of the specimen temperature and 
hardness impressions were made with a 136° 
pyramidal diamond loaded by a rod and scale 
pan outside the tank, the rod operated through 
a Wilson seal in a top plate of the chamber 
which was water cooled. The impressions were 
measured with an eyepiece micrometer after 
removing the specimen from the apparatus. 

A load of 5 kg was used from room temper- 
ature to 800° C. Two indentations were made 
at room temperature and from 100° C—800° C 
at intervals of 100°C. The spacing between 


each indentation was 2-3 times the diagonal 
of the impression and the load was applied for 
30 seconds; results on any one specimen at any 
temperature were reproducible to + 5 VPN. 
The heating time to the test temperature was 
approximately ten minutes and the sample was 
held at temperature for an equal time before 
the indentation was made. 


5. Metallographic Examination 


Standard methods of metallographic prepa- 
ration were employed and grain sizes were 
revealed by ionic bombardment in an argon 
atmosphere 3), All specimens were examined 
when testing was completed, and in addition 
the binary thorium-aluminium and thorium- 
uranium alloys prepared from production 
thorium were examined at intervals during the 
ageing tests. 


6. Ageing Results 


Hardness values are plotted against log time 
in figs. 2-7 and the details are described below. 
After solution treatment, alloys prepared from 
production bar thorium were about 5 VPN 
harder than corresponding specimens made 
from electrolytic material. 


6.1. THORIUM-ALUMINIUM: (figs. 2-4) 


All the thorium-aluminium alloys showed 
marked precipitation hardening behaviour. The 
hardness in the _ solution-treated condition 
increased with aluminium content and there 
was little difference between the ageing results 
obtained on 1 at 9% aluminium alloys prepared 
from the two types of thorium. The rate of 
ageing increased directly with aluminium con- 
tent and with temperature in the range 
500—700° C, while the maximum hardness for a 
given composition increased inversely with 
temperature. At 300° C age-hardening had not 
begun after 1000 h. Over-ageing occurred 
rapidly at 700° C for each composition but at 
500° C no significant drop in hardness of the 
1 or 2 at % aluminium alloys had occurred 
1000 h. The maximum hardness obtained was 
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Fig. 2. Ageing curves for thorium/1 at % aluminium, prepared from electrolytic 
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Fig. 4. Ageing curves for thorium/5 at % aluminium, prepared from production thorium. 
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Fig. 5. Ageing curves for thorium/2 at % uranium, prepared from electrolytic thorium. 
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Fig. 6. Ageing curves for thorium/3 at % uranium, prepared from production thorium. 
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199 VPN with 5 at % aluminium after 264 h 
at 500° C. 

Grain sizes of approximately 3 mm and 
0.38 mm were observed with 1.0 and 2.0 at % 
aluminium respectively. Since unalloyed thorium 
also has a grain size of approximately 3 mm 
after similar thermal treatment and the solu- 
bility of aluminium in thorium at 1300°C is 
0.85 at %, it appears that 1.0 at % aluminium 
is sufficiently close to the solubility limit to 
have little effect on grain growth. 


6.2. 

These alloys showed only a slight increase of 
hardness on ageing, the hardness in the solution 
treated condition increasing with the amount 
of alloying element. There was no significant 
difference between the results obtained on 
2 at % uranium alloys prepared from the two 
types of thorium. The rate of ageing varied 
with temperature but did not alter with compo- 
sition; overageing had begun after +h at 700° C 
but at 500°C annealing times approaching 
1000 h were required before peak hardness was 
exceeded. The hardness increase did not alter 
with composition or ageing temperature and 
maximum hardness values were 112 and 123 
VPN for 2 at % and 3 at % uranium respec- 
tively. Metallographic examination showed that 
uranium in excess of the solubility limit reduces 
grain growth at 1200°C. 


THORIUM-URANIUM: (figs. 5 and 6) 


6.3. 


The hardness of all the complex alloys 
increased on ageing at 600° C. They each aged 
at the same rate for the first 8 hours and soften- 
ing occurred slowly with further annealing of 
all except the thorium-aluminium-uranium 
alloy; this gave the highest peak hardness of 
163 VPN and had a hardness of 156 VPN after 
256 hours. 


COMPLEX ALLOYS: (fig. 7) 


7. Hot Hardness Results 


The variation of hardness of thorium and 
thorium alloys with temperature is shown in 
figs. 8-11 and the more important features are 
discussed below. 


7.1. THorium (fig. 8) 

The initial higher hardness of the production 
thorium in the “as-received”? condition is 
probably associated with the light reduction 
by swaging given to this material after sintering 
and the lower hardness of electrolytic thorium 
below 500° C, to its lower carbon content. All 
three materials show a tendency to a hardness 
plateau in the temperature range 200—400° C, 
the effect being most pronounced in the 
“as-received”? production bar. Figs. 9 and 10 
give an indication of the distribution of the 
oxide inclusions in the arc-melted samples of 
production and electrolytic thorium, and these 
are also typical of alloys prepared from these 
materials; the corresponding grain sizes of the 
two types of unalloyed thorium, after testing, 
were 30 uw and 80 wp respectively. 


7.2. ALLOYS PREPARED FROM PRODUCTION 


THORIUM (fig. 11) 


These results, obtained by R. C. Burnett and 
S. F. Pugh, have been included here since they 
have not been reported previously. There is little 
difference between the hot-hardness values for 
the thorium-zirconium, thorium-zirconium-ura- 
nium and thorium-zirconium-titanium alloys 
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Fig. 8. Hot hardness curves for thorium. 
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Fig. 9. Arc-melted production thorium. Mechanical Fig. 10. Arc-melted electrolytic thorium. Mechanical 
polish. ThOzg dendrites dark. x 150. polish. ThOz dendrites dark. x 150. 


tested in the arc-melted condition; all showed 


A PRODUCTION Th AS ARC-MELTED a tendency to a hardness plateau in the range 
° ‘ Th/12 Yo Ze Rg gee 200°-500° C, (with 16 at °% zirconium the 
x ” ThH/i2% Zr 185% U AS ARC-MELTED e : 

; Z : ca plateau extends to 600° C) and all were signifi- 
Oo " Th/\1% Zr 8 %o Ti AS ARC-MELTED i. 
° : th/2% w Kae ree cantly harder than unalloyed thorium over the 


1 4 PRODUCTION Th/l6% zr AS ARC-MELTED whole range of temperature. The reason for the 
; drop in hardness of the thorium-uranium alloy 
between 200° C and 300° C is not clear; much 
of the uranium was present as a discontinuous 
grain boundary network and in view of the long 
annealing times required to redistribute uranium 
in thorium below 800°C?) no such changes 
would be expected to occur during hardness 
testing. 
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7.3. ALLOYS PREPARED FROM ELECTROLYTIC 
THORIUM 
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These were all examined in the following 
conditions: (a) as-cast (b) as solution treated, 
and (c) as solution treated and aged for seven 
days at 600° C. 

(i) Thorium-zirconum: (figs. 12-14). Each 
of the three alloys was harder than unalloyed 
thorium in all conditions and below 600° C 
was softer after solution-annealing than as- 
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Fig. 11. Hot hardness curves for alloys made with pared from electrolytic (fig. 13) and from 
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cant differences in spite of the different carbon 
contents of the materials. After solution treating 
at 900° C, the 2.5 and 12 at % zirconium alloys 
showed a tendency to a hardness plateau in the 
temperature range 100°-500° C; with 22 at % 
zirconium this range was increased to 100°— 
600° C. 

Solution treatment of the 22 at % zirconium 
alloy in the body centred cubic (8) phase field 
at 1100° C produced a marked increase in hard- 
ness at temperatures up to 600°C but on 
subsequent ageing at 600°C the hardness 
dropped appreciably although it remained 24 
times higher than that of unalloyed thorium. This 
composition is the minimum zirconium addition 
to enable # solution treatment at 1100°C; at 
900° C all compositions fall in the («+ f)field. 

(ii) Thoriwm-aluminium : (figs. 15-17). Below 
700° C, in the as-cast condition, there was a 
considerable increase in the hot hardness with 
aluminium content. This might be related to the 
amount of Th-TheAl eutectic present in the 
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aluminium. 


grain boundaries since the solubility of alu- 
minium in thorium is less than 1 at % at 
1300° C1). After solution treating there was no 
significant hardness differences between the 
three compositions, possibly because the TheAl 
in the eutectic had started to “ball-up” and 
each alloy contained aluminium in excess of the 
solubility limit. After ageing for seven days at 
600° C the 1 and 5 at % alloys had higher hot 
hardness values than in the solution treated 
condition at all temperatures tested. Both alloys 
showed evidence of a strain-ageing peak at 
300°-400° C and by analogy with steel 4) 
might have particularly good creep properties 
at 250°-350° and be significantly better than 
unalloyed thorium up to 500° C. Metallographic 
examination showed that considerable grain 
growth had occurred during solution treatment 
of each alloy and grain sizes of the 1, 5 and 8 
at % aluminium alloys after testing were 
approximately 2 mm, 0.5 mm and 0.25 mm 


respectively. With 1 at % aluminium there 
were a few particles of Th2Al which had not 
been in solution at 1150° C and a fine precipitate 
of Th2Al was observed in the grain boundaries 
and also within the grains. Less fine precipitate 
was present in the 5 and 8 at °% alloys possibly 
because of the initial presence of many TheAl 
particles which could act as nuclei for subse- 
quent precipitation. 

(iii) Thorium-uranium : (fig. 18). The alloy 
containing 2 at °% uranium was slightly harder 
than that containing 3 at % uranium in the 
as-cast condition. The reason for this is not 
clear but might be due to a different amount of 
superheat and hence cooling rate in the two 
melts. The hot-hardness results for both 
materials were similar after solution treating 
and after ageing; in each case the aged curve 
was lower. The curves for the 3 at % uranium 
alloy are given here as typical. 
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uranium, 


Grain sizes determined after the completion 
of testing were approximately 1 mm and 0.1 mm 
respectively for the 2 at % and 3 at % alloys; 
precipitation occurred in the grain boundaries 
and within the grains of both alloys. With 
3 at % uranium, the excess uranium present 
during solution treatment was associated with 
the grain boundaries, and with both 2 and 
3 at °4 uranium the boundaries showed some 
association with ThOs particles. 


8. Discussion 


Alloys of thorium and uranium exhibit a 
considerably smaller hardness increase on ageing 
than those of thorium and aluminium although 
the change in solubility with temperature is 
greater for the former. In thorium-uranium 
alloys the precipitate is B-uranium at 700° C 
and «-uranium below 660° C; neither modifi- 
cation can form coherently on the «-thorium 
matrix and a considerable degree of misfit must 
occur. The precipitating phase in thorium- 
aluminium alloys is TheAl. Since «-thorium, 
like aluminium, is face-centred-cubic at the 
ageing temperature and ThoAl is isostructural 
with CuAlg, there is a possibility that precipi- 
tation in thorium-aluminium alloys is multi- 
stage, as is the case with aluminium-rich alloys 
in the aluminium-copper system; if it were so, 
then one stage might be the formation of a 
coherent precipitate with consequent increase 
in hardness. It was not possible to determine 
categorically if there was an initial small 
hardness peak in the present series of alloys; 
a slight increase was observed at 300°, 500° and 
600° C before the main hardening began but 
this was only slightly greater than the experi- 
mental error. In thorium-aluminium alloys, 
grain boundary precipitation was first detected 
metallographically after ageing for 24 h at 
600° C, i.e. when the hardness had just started 
to increase. 

Except for the thorium-aluminium-uranium 
alloys, all the ternaries showed evidence of 
slight grain-boundary liquation after solution 
treatment at 1200° C and therefore, assuming 
maximum solubility at the eutectic tempera- 
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ture, the amount of alloying elements in solid 
solution was less than this maximum. Since the 
constitution of none of the complex alloys has 
been studied over the relevant range of compo- 
sition, the eutectic temperatures are not known 
nor is it possible to predict the nature of the 
precipitating phases. 

No attempt has been made to determine the 
effect of the rate of cooling from 1200°C on 
the subsequent ageing behaviour of any of these 
alloys. In aluminium-rich aluminium-copper 
alloys Prakash and Entwistle >) have shown that 
slower rates of cooling from the solution anneal 
lead to failure to retain vacancies and hence a 
lower hardening rate. The effect of grain size 
on the rate of ageing has not been investigated ; 
the complex alloy which overaged least rapidly 
had the coarsest grain size, but it is not known 
whether the slower rate of overageing was due 
to the coarse grain size or to the nature of the 
material. 

Since alloys which show appreciable ageing 
at the upper temperatures of hot hardness 
testing would give a hot hardness result which 
was a function of thermal history, the conditions 
of solution treatment, of ageing and of hot 
hardness testing were standardised in order 
that the hardness results might give a useful 
comparison between different materials. An 
ageing treatment of seven days at 600° C prior 
to hot hardness testing was deliberately chosen 
so as to overage the alloys, in an attempt to use 
the hot hardness tests to assess the possible 
long-term creep behaviour of these materials at 
lower temperatures. 

The hardness results after solution treatment 
at 900°C suggest that zirconium additions 
within the «-thorium solid solution range lead to 
solution hardening, while the similar hardness 
of alloys after ageing at 600° C indicates that 
no capacity for age-hardening remains after this 
treatment. (The solid solubility of zirconium in 
thorium at 900° C and 600° C is 8.2 at % and 
<5 at % respectively °).) 

A few previous hot-hardness results are 


available on thorium and thorium alloys. 
Goldhoff, Ogden and Jaffee’) examined a 
series of binary alloys, (all of which had been 
annealed 2 h at 850° C), at room temperature 
and at 300°C. Their values for unalloyed 
thorium at the two temperatures were inter- 
mediate between those obtained on production 
and electrolytic thorium in the present investi- 
gation; figures for a 1 at % aluminium alloy 
show rather greater temperature dependence 
than the present results on solution treatment 
material and no relation to the figures at 300° C 
on aged material, while their values for 1 at % 
zirconium are higher than the present results on 
2.5 at % zirconium. Peterson, Russi and 
Nickelson 8) report the as-cast hardness of a 
1 at % aluminium alloy as 84 VPN, which 
accords with 87 obtained here. The room 
temperature hardness of the thorium/3 at % 
uranium alloy showed little variation with heat 
treatment and was in agreement with Bentle’s 9) 
figures for a commercial purity alloy of this 
composition after a variety of heat treatments. 
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The diffusion of hydrogen in the alpha-phase of 
Zircaloy-2 was measured in the temperature range 
260° to 560°C using the gradient technique. The 
diffusion coefficient was found to be 

D = 2.17 x 10-3 exp (— 8380/RT) cm2/sec. 

The terminal solid solubility of hydrogen 

Zircaloy-2 was determined in the temperature range 
260° to 650°C using a modification of the gradient 
technique. The solubility is given by 

Co = 8.50 x 104 exp (— 7600/RT) ppm He 
by weight. 


in 


La diffusion de Vhydrogéne dans la phase « du 
zircaloy-2 a été mesurée entre 260 et 560°C par la 
technique du gradient. Le coefficient de diffusion 
déterminé par cette méthode est 
D = 2,17 X 10-3 exp (— 8380/RT7’') cm2/sec. 
La solubilité terminale dans |’état solide de Vhydro- 


1, 


Introduction 


Zircaloy-2 has the properties of low corrosion 
rate, reasonable strength up to at least 300° C 
and low neutron absorption, all of which make 
it desirable as a structural material in water 
cooled power reactors. One of its chief dis- 
advantages, however, is a drastic reduction in 
impact strength brought by about the presence 
of the hydride phase). A knowledge of the 
hydrogen distribution and the terminal solid 
solubility is therefore necessary to predict the 
useful life of Zircaloy in the reactor. 

As the Zircaloy-2 pressure tubes in the NPD 
(Nuclear Power Demonstration) reactor will 
operate at 260° C it was necessary to know the 
diffusion rate and the solubility limit of 
hydrogen in the a-phase of Zircaloy-2 at this 
temperature. Schwartz and Mallett 2), Mallett 
and Albrecht 3), and Gulbransen and Andrews 4) 
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géne dans le zircaloy-2 a été déterminée entre 260° et 
650° C en utilisant une variante de la technique du 
gradient. La solubilité est donnée par 

Co = 8,50 x 104 exp (— 7600/RT) ppm de He 
en poids. 


Die Diffusion von Wasserstoff in der «-Phase von 
Zirkaloy-2 wurde im Temperaturbereich von 260° bis 
560°C unter Anwendung der Gradientenmethode 
untersucht. Der Diffusionskoeffizient wurde zu 

D = 2,17 x 10-3 exp (— 8380/RT) cm2/sec 
ermittelt. 

Die maximale Léslichkeit von Wasserstoff in 
Zirkaloy-2 wurde zwischen 260° und 650° C bestimmt, 
wobei eine abgewandelte Form der Gradientenmethode 
verwendet wurde. Die Léslichkeit wird in Gew.-ppm 
gegeben durch 

Co = 8,50 x 104 exp (— 7600/RT). 


measured the terminal solid solubility of 
hydrogen in zirconium, but only Mallett and 
Albrecht made measurements at temperatures 
as low as 305° C. Extrapolation of these data 
to 260° C gave terminal solid solubilities of 75, 
45 and 28 parts per million (ppm) by weight 
respectively. Since the last value is similar to 
the amount of hydrogen normally present in 
Zircaloy-2 tubing, it was important to establish 
the terminal solid solubility of hydrogen in 
Zircaloy-2 at 260° C. 

Although the degassing method employed by 
Eborall and Ransley *) for determining the 
diffusion of hydrogen in an Al-Mg alloy was 
tried, it was decided that the technique used by 
Mallett and Albrecht?) for measuring the 
diffusion of hydrogen in zirconium would give 
more reliable results in the temperature range 
of interest. Essentially, the method consists of 
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determining the hydrogen distribution in a 
Zircaloy cylinder annealed in hydrogen for a 
predetermined time. 

The diffusion experiment yielded values of 
terminal solid solubility, but more accurate 
results were obtained by a modification of this 
technique. A thin hydride film was formed on a 
Zircaloy cylinder by reaction with hydrogen 
and annealed for a time sufficient to homogenize 
the alpha phase. Hydrogen analysis then gave 
the terminal solid solubility directly. 


2. Experimental Procedure 


2.1. MATERIAL 


The two sources of Zircaloy-2 with composi- 
tions given in table 1 were in the form of rods 


TABLE 1 


Analysis of the Zircaloy-2 
Alloy analysis in wt percent 


Element Ingot FZ-644 Ingot K-28 
Sn 1.36 1.51 
Fe 0.13 0.15 
Cr 0.10 0.10 
Ni 0.05 0.05 

Impurities in ppm by wt 
N 35 31 
C 90 < 300 
Si ah 30 
Al 38 21 
Hf < 100 72 
Cu < 25 < 20 
Vv < 25 < 20 
Mn 35 < 20 
Mg < 10 < 20 
Pb < 20 < 20 
Mo < 20 < 20 
Co < 10 < 20 
W | < 20 
B < 0.3 = 0.2 
Cd V4 < 0.5 
O < 1200 < 1400 
H 45 32 
swaged to 0.5 in. diameter. Micrographic 


examination showed that the FZ-644 batch was 
in the as-swaged conditions, whereas the K-28 
material had been annealed in the alpha region. 


2.2. APPARATUS 


A schematic diagram of the apparatus is 
shown in fig. 1. Except for several modifications, 
it is the hot-extraction apparatus described by 
McGeary °). 

The sample furnace consisted of nichrome 
ribbon wound on a Vycor tybe, and was capable 
of reaching temperatures up to 1100°C. The 
manual temperature control using a Variac 
produced variations of less than + 1°C. The 
working zone of the furnace had a maximum 
temperature variation of 3° C over 1.5 inches. 
A chromel-alumel thermocouple made from 
National Bureau of Standards calibrated wire 
was placed against the bottom of the sample- 
tube and calibrated to read sample temperature. 
The error in specimen temperature was approxi- 
mately + 5°C. 

Each of the two high-speed mercury diffusion 
pumps | and 2 had a large cold trap sealed to 
the high vacuum outlet. Pump No. 1, used to 
transfer the hydrogen from the sample tube into 
the analysis volume, attained a pressure of less 
than 0.01 microns with a back pressure of 
1500 microns. 


The total analysis volume consisted of the 
1050-ce flask, the 2000-cce flask, the McLeod 
gauge, the tubing and the low vacuum side of 
pump No. 1. This volume was calibrated by 
introducing hydrogen into the analysis volume 
and measuring the pressure with the McLeod 
gauge. The 1050-cc flask, whose volume had 
been accurately determined before assembly, 
was then evacuated and the hydrogen contained 
in the rest of the volume expanded into it. By 
measuring the final pressure and applying the 
perfect gas law the total volume was calculated. 
Suitable combinations of the two flasks A and 
B (fig. 1) together with the associated tubing 
between diffusion pump No. 1 and the McLeod 
gauge enabled one to have four different 
analytical volumes. 


The palladium valve used to introduce pure 
hydrogen into the system was a spiral of 
palladium tubing with a 0.064 cm wall, 0.32 em 
diameter and 30 cm length. One end of this 
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spiral was silver-soldered to the Kovar tip 
joined to a short length of glass tubing, and the 
other end was closed. With the palladium tube 
at 400°C a hydrogen pressure of 100 mm in 
the system could be attained in one hour. 


2.3. THE DIFFUSION EXPERIMENT 


A Zircaloy cylinder 1 cm in diameter and 
3 cm long was degassed in vacuum for six hours 
at 830° C. A residual hydrogen concentration of 
5 to 6 ppm by weight as determined from the 
equilibrium curves, remained in the specimen 
after this treatment. The sample was abraded 
with 600-A grit carbide paper and placed in 
the sample tube of the furnace. It was then 
heated in vacuo for 0.5 h at 830° C to remove 
any oxide film and brought to the diffusion 
temperature. Hydrogen, at a pressure of 200 mm 
was next introduced into the sample tube; this 
pressure was sufficient to maintain a hydride 
film on the specimen surface. The sample was 
kept at the diffusion temperature for a sufficient 
time to give a suitable hydrogen distribution, 
after which the furnace was quickly removed 
and the sample tube quenched in liquid nitrogen. 
Lengths equal to the diameter were removed 
from the end of each sample in order to eliminate 
effects due to the finite specimen length. The 
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remaining portion was machined into seven 
layers of equal weight. The outer layer con- 
taining the hydride was discarded and the 
remaining six analysed for hydrogen. 

A plot of concentration versus radius when 
extrapolated to the sample surface gave the 
terminal solid solubility Co. The diffusion 
coefficient was determined from the plot of 
(C—O,)/(Co—C1) versus r/a where: 

C =Hydrogen concentration at radius r 

C;= Initial uniform hydrogen concentration 

Co=Surface concentration of hydrogen in the 
ax-Zr phase 


a =radius of sample. 


2.4. DETERMINATION OF THE TERMINAL SOLID 
SOLUBILITY 


Since the terminal solid solubilities obtained 
by the above method are susceptible to fairly 
large errors due to extrapolation, a variation of 
this technique leading to more accurate results 
was also used. 

The Zircaloy used in these measurements was 
cold-worked 40 % to eliminate possible effects 
due to differences in the microstructure of the 
two batches of material. The specimens were 
2 cm in length and 0.5 cm in diameter except 


i 
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for those treated at 260° C which were 0.4 cm 
in diameter. 

Two samples, one from each batch, were 
placed together in the furnace and given the 
same degassing and oxide-removal treatment 
as the diffusion specimens. With the samples 
at the desired temperature, hydrogen at a 
pressure of 150 mm was introduced into the 
sample tube. After the samples had absorbed 
enough hydrogen to form a hydride skin 
approximately 0.005 cm thick, the pressure was 
reduced to a value 7) which was just sufficient 
to maintain the hydride film without causing 
further growth. The samples were held at the 
diffusion temperature for a time t=a?/D where 
D is the diffusion coefficient at operating 
temperature. This time was long enough to 
give a uniform hydrogen distribution in the 
alpha phase to within 1 °4. Sample dimensions 
had been chosen to allow homogenization in 
reasonable times. After removal from the 
furnace a layer 1 mm thick was machined from 
the specimen surface to remove the hydride 
film, and the remaining portion analysed for 
hydrogen. The resultant concentration was the 
terminal solid solubility. 


2.5. 


Since the technique used results in effectively 
radial diffusion, the diffusion equation becomes: 


ANALYSIS OF THE DIFFUSION DATA 


E2208) 
where c — concentration 
¢ — time 
ry — radius 
D — diffusion coefficient. 
This equation must be solved using the 
following boundary conditions: 
fy) c=Cprat f=o tort >0 
Ze; lor O--7- a at i=—0. 
As shown in Crank 8) the solution to eq. (1) is 


ghee Tolrn) 
m= 1G D exp (— Don?t) 


2) 


where Jo(x) and J;(2) are first-kind Bessel 


functions of zero and first order respectively, 
and the «»’s are the positive roots of Jo(axn) = 0. 
The solution to eq. (2) can be written in terms 
of the two dimensionless parameters Dt/a? and 
rja. Curves showing (C—C;)/(Co—Ci) as a 
function of r/a for different values of Dt/a? may 
be found in Crank’s book. 

In _ practice the experimental values of 
(C—C})/(Co—Ci) are plotted against r/a and 
compared with the theoretical curves. A value 
for Dt/a? is thus obtained and, since ¢ and a 
are known, D can be calculated. The experi- 
mental conditions are chosen so that a thin 
hydride film is formed on the sample. The 
hydrogen concentration in the Zircaloy that is 
in equilibrium with the hydride, namely the 
terminal solid solubility, is then Co in the 
above equation. 


3. Results and Discussions 
Bad be 


The diffusion of hydrogen in the alpha-phase 
of Zircaloy-2 was investigated in the range 
260° to 560°C. The plot of (C—C\)/(Co—C1) 
against r/a for 358° C, which is typical of the 
curves obtained, is shown in fig. 2. The points 
are the experimental values and the line is the 
best theoretical fit from which D was deter- 
mined. Scatter for small values of (C—C})/ 
(Co —C}) as observed in most of the curves is to 
be expected since the relative error increases 
as C approaches (4. 


DIFFUSION COEFFICIENTS 


TABLE 2 
Diffusion of hydrogen in Zircaloy-2 


Temp. | Time | Batch D x 108 | Co (ppm by wt) 
(°C) | (min) No. (em?/sec~!) | uncorr.| corr. 
261 150 K-28 0.80 75 71 
315 90 K-28 1.70 Ei ss 105 
358 60 K-28 2.94 300 280 
408 30 FZ-644 4.71 180 145 
458 15 | FZ-644 5.49 490 450 
515 1S: K-28 10.3 800 700 
560 15 K-28 15.8 — — 

The diffusion coefficients obtained at the 


various temperatures are tabulated in table 2. 
The values of Co obtained by extrapolating the 


66 A. SAWATZKY 


concentration curves to the sample surfaces 
are also given. 

The assumption is usually made 3) that the 
hydride film is negligibly thin. Using the data 
of Gulbransen and Andrew 9) for the diffusion 
of hydrogen in zirconium hydride, the hydride 
thickness was determined, and a value for Co 
found by extrapolating the concentration curve 
only to the hydride boundary. As shown in 
table 2 the correction was appreciable. These 
corrected values fell reasonably close to the 
Co's obtained by the second method used, 
except for the point at 408°C. The uneven 
hydriding of this specimen may explain the 
very low value for Co. A corrected D was 
obtained using the corrected Cp and assuming 
the average specimen radius to be midway 
between the specimen surface and the hydride 
boundary. The correction here was found to be 
negligible as shown in fig. 2. 

The logarithm of D versus 1/7’ is plotted in 
fig. 3. A least-squares fit of the experimental 
data gives an activation energy Q = 8380+ 400 
cal/mole and a diffusion coefficient 


D=2.7 Xx 10-3 exp (— 8380/R7’) cm2/sec. 
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Fig. 3. 


In the same figure the present data are com- 
pared with those of Schwartz and Mallett?) and 
Mallett and Albrecht®) for the diffusion of 
hydrogen in zirconium. Although the diffusiv- 
ities at any one temperature are quite similar 
in the three cases, the activation energy obtained 
in the present investigation (8380 cal/mole) is 
somewhat higher than those of either Schwartz 
and Mallett (5700 cal/mole) or Mallett and 
Albrecht (7060 cal/mole). It should be noted 
that Schwartz and Mallett’s result is based on 
only three measurements. The two different 
activation energies obtained for hydrogen 
diffusion in zirconium are possibly due to 
different amounts of impurities such as oxygen 
or nitrogen in the zirconium used. Zircaloy-2 
contains the four alloying elements shown in 
table 1, which undoubtedly contribute to the 
difference in activation energy between hydrogen 
diffusion in zirconium and Zircaloy-2. 


3.2. TERMINAL SOLID SOLUBILITIES 


Terminal solid solubilities were measured in 
the range 260° to 650° C by the second method 
described earlier. Table 3 lists the solubilities 
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TABLE 3 
Terminal solid solubilities for H in Zircaloy-2 


Temp. | Time Co (ppm by wt) Co 
(°C) (hours) FZ-644 K-28 (average) 
257 20 79 81 
257 18 76 50 71 
300 1] 93 95 
300 ll 87 
300 12 99 101 95 
356 11 212 218 
356 13 189 191 203 
406 + 273 268 270 
460 3.5 478 468 
457 7 480 465 473 
500 4 635 640 638 
560 4 888 886 887 
600 2 3075 3075 
650 2 3200 3170 3185 


obtained for the two batches of Zircaloy tested. 
In general the reproducibility is quite good, and 
it was also found that the terminal solid 
solubility for the two materials is essentially 
the same. 

The logarithm of Co versus 1/7’ is plotted in 
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fig. 4. A least-squares fit of the experimental 

data yields 

Ay= 7600+ 200 cal/mole and 

Co =8.50 x 104 exp (—7600/RT) ppm He by 
weight. 

Ay represents the difference in partial molar 

heats of solution of hydrogen in the 6(hydride) 

and « phases. The least-squares value of the 

solubility at 260°C is 64 ppm whereas the 

average of four determinations at 257°C is 

71 ppm. 

In fig. 4, the solubility of hydrogen in Zircaloy- 
2 found in the present investigation is compared 
with the published values for the solubility of 
hydrogen in zirconium. The data of Schwartz 
and Mallett 2), Mallett and Albrecht?) and 
Gulbransen and Andrew‘) extrapolated to 
260° C lead to terminal solid solubilities of 75, 
45 and 28 ppm He: by weight respectively. 
Although the present solubility of hydrogen in 
Zircaloy-2 (64 ppm) agrees best with that of 
Schwartz and Mallett, the latter covered only 
the temperature range 400° to 500°C making 
extrapolation to 260° C rather unreliable. How- 
ever, on the basis of probable errors in the 
hydrogen analyses and temperatures, the 
present solubilities are accurate to within 
+ 10%. 

It will be seen that at 600° C and 650° C the 
results fall off the straight line shown in fig. 4. 
Metallographic examination of these two samples 
showed that they had been annealed in the beta 
region. This shows that the «-phase was not in 
equilibrium with the d-phase at these temper- 
atures, in agreement with the equilibrium 
diagram of Ells and McQuillan ”). 

On the basis of this diagram, in a diffusion 
experiment carried out between 560° to 850° C, 
one should get a hydride (é6-phase) layer on 
the specimen surface, followed by a layer of 
B (bec) phase, the rest of the specimen being 
« (hep) phase. If the annealing was carried 
out for a long time, the «-phase region would 
gradually give away to the 6-phase. The speci- 
mens corresponding to the points at 600° C and 
650° C in fig. 4 were annealed for two hours, 
which left the major portion of these specimens 
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in the 6-phase, thus accounting for the observed 
higher hydrogen concentrations. In contrast, 
the two samples at 600° C and 650° C studied 
by Mallett and Albrecht were annealed for 
approximately five minutes. The major portion 
of their samples would be in the «-region, giving 
values of Co nearly equal to the terminal 
solubility at 560°C. As seen in fig. 4, their 
values for Co at temperatures above 560° C are 
slightly lower than expected according to Ells 
and McQuillan. 


4. Conclusions 

The _ diffusion coefficient of hydrogen in 
Zircaloy-2 is 

D=2.17 < 10-3 exp (— 8380/R7') cm2/sec. 

The terminal solid solubility of hydrogen in 
Zircaloy-2 is 

Co = 8.50 x 104 exp (— 7600/RT') ppm He 
by weight. 
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ETUDE AUX NEUTRONS DE LA TEXTURE CRISTALLINE DE BARREAUX D’URANIUM 
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L’étude des textures par diffraction de neutrons a été 
appliquée a un élément de barreau d’uranium filé en 
phase « a 600° C. Dans le cas de détermination de 
texture par diffraction de rayons X, pour les corps 
absorbants, la couche intéressée par le faisceau est 
faible: il est done nécessaire de faire les analyses sur 
un grand nombre de coupes pour obtenir un résultat 
statistique. Les neutrons étant peu absorbés par la 
plupart des substances de poids atomique élevé, il est 
possible d’obtenir en une seule série d’essais une 
image statistique correcte de la texture. L’échantillon 
est un cylindre de 5 cm? environ dont l’axe, perpendi- 
culaire a la direction de filage, coincide avec celui 
d’un faisceau de neutrons monochromatiques (A = 
1,143 A). Le détecteur de neutrons est un compteur 
BF3 qui décrit dans un plan horizontal un cercle 
centré sur léchantillon. Un spectre de diffraction 
effectué sur un élément de barre d’uranium filé en 
phase « et traité en phase fa montré que cet échantillon 
ne présentait qu’une faible orientation préférentielle. 
Par contre, sur l’échantillon brut de filage il apparait 
une texture marquée [110] dans la direction de filage. 
Les résultats obtenus sont en accord avec ceux 
déterminés par diffraction de rayons X par différents 
auteurs sur des barres d’uranium filées en phase « 
ou laminées en phase «. 


The method of studying textures by neutron diffraction 
has been applied to uranium bar extruded in the « 
range at 600° C. When textures are determined by X-ray 
diffraction, the thickness of metal through which the 
X-ray beam can pass is slight in the case of absorbent 
materials; it is therefore necessary to carry out 
analyses on a large number of sections to obtain a 
statistically significant result. Since neutrons are not 
strongly absorbed by most substances of high atomic 
weight, it is possible by means of a single series of 
experiments to obtain a good statistical image of the 
texture. The specimen is a cylinder of about 5 cm3, 
the axis of which, normal to the direction of drawing, 
coincides with the axis of a beam of monochromatic 
neutrons (A = 1.143 A). The neutron detector is a 
BF3 counter which describes a circle in a horizontal 
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plane, centred on the specimen. The diffraction spec- 
trum obtained from a piece of uranium bar, extruded 
in the «-range and treated in the f-range, has shown 
that this specimen has only a weak degree of preferred 
orientation, whereas an as extruded specimen shows a, 
marked [110] texture in the direction of extrusion. 
The results obtained are in agreement with those 
determined by X-ray diffraction by various authors 
with uranium bars extruded or rolled in the a-range. 


An einer Probe aus einem Uranstab, der bei 600° C 
in der x-Phase durch Strangpressen hergestellt wurde, 
ist eine Texturuntersuchung durch Neutronenbeugung 
vorgenommen worden. Fuhrt man Texturbestim- 
mungen mit Réntgenstrahlen durch, so ist fiir absor- 
bierende Koérper die von dem Strahlenbtindel erfasste 
Schicht sehr diinn: daher ist es notwendig, die Unter- 
suchungen an zahlreichen Schnitten vorzunehmen, 
um ein statistisches Ergebnis zu erhalten. Da dem 
gegentiber Neutronen von den meisten Elementen 
mit hédherem Atomgewicht nur wenig absorbiert 
werden, kann man in einer einzigen Versuchsreihe 
ein richtiges statistisches Bild von der Textur erhalten. 
Die verwendete Probe ist ein Zylinder von etwa 
5 em? Volumen, dessen Achse senkrecht zur Strang- 
pressrichtung steht und mit der Achse des Bundels 
monochromatischer Neutronen (A = 1,143 A) zu: 
sammenfallt. Als Neutronenindikator diente ein 
BF3-Zahler, der in horizontaler Ebene einen Kreis 
mit der Probe als Mittelpunkt beschrieb. Aus dem 
Beugungsspektrum einer Probe, die aus dem durch 
Strangpressen in der «-Phase hergestellten und in der 
B-Phase wiarmebehandelten Uranstab stammte, war 
zu entnehmen, dass eine soleche Probe nur eine 
schwache Orientierung aufweist. Dagegen zeigt sich 
bei einer Probe, die direkt nach dem Strangpressen 
entnommen wird, eine ausgepragte [110]-Textur in 
Strangpressrichtung. Die erhaltenen Befunde stimmen 
luberein mit LErgebnissen aus Réntgenbeugungs- 
messungen, welche von verschiedenen Autoren an 
stranggepressten oder gewalzten Uranstaben vorge- 
nommen worden sind. 
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1. Barreaux bruts de filage en phase a 


La déformation de barreaux d’uranium sous 
irradiation est conditionnée, toutes choses 
égales par ailleurs, par deux facteurs essentiels: 
la grosseur des grains et leur orientation. 

De nombreuses expériences ont montré que, 
pour limiter les déformations, il était nécessaire 
(mais non suffisant) d’avoir des barreaux formés 
d’agrégats polycristallins a la fois homogénes et 
exempts d’orientations préférentielles. 

On est conduit pour se rapprocher de ces 
conditions & déformer l’uranium en phase a, 
& le soumettre ensuite a des _ traitements 
thermiques en phase 6, qui éliminent les orien- 
tations préférentielles dues a la déformation, 
puis en phase «, qui provoquent l’affinage par 
recristallisation et éliminent les contraintes 
résiduelles. 

Les résultats ont montré que les déformations 
sous irradiation de barreaux ainsi traités 
étaient effectivement nettement inférieures a 
celles de barreaux d’uranium coulés ou filés en 
phase « ou y sans recuit, sans étre toutefois 
entieérement éliminées. 

Il reste probablement aprés_ traitements 
thermiques un résidu d’orientations préféren- 
tielles complexes, d’ot: lintérét d’étudier |’évo- 
lution de la texture des barreaux d’uranium au 
cours des différents traitements thermiques. 

L’étude des orientations préférentielles dans 
un échantillon polycristallin s’effectue commo- 
dément au moyen des rayons X, dans la mesure 
ou l’échantillon n’est pas trop absorbant. Dans 
le cas contraire, différents procédés fort ingénieux 
ont été mis en ceuvre !-®). Ils présentent tous 
cependant les inconvénients suivants: 

a) le faisceau n’intéresse qu’une faible 
couche superficielle de |’échantillon, 

b) il en résulte que le nombre des micro- 
cristaux diffusants est faible, d’autant plus 
faible que le grain est plus gros, et qu’on doit 
souvent répéter les analyses sur différentes 
coupes de l’échantillon pour obtenir une image 
statistique correcte de la texture. 

A ce point de vue, les neutrons présentent 
Vavantage, par rapport aux rayons X, d’étre 


peu absorbés, sauf cas exceptionnels, par les 
substances constituées d’atomes de poids atomi- 
que élevé. L’étude d’une texture cristalline au 
moyen des neutrons présente alors les avantages 
suivants: 

a) Tous les cristallites d’un échantillon 
massif dont le volume peut atteindre plusieurs 
cm3 participent a la diffusion. 

b) On peut donc obtenir en une seule série 
d’essais une image correcte de la texture, la 
statistique portant sur un nombre de grains 
suffisant. 

La méthode a déja été mise en ceuvre pour 
étudier l’influence du recuit sur la texture 
cristalline de fils métalliques ’). Au moyen du 
spectrométre 4 neutrons de la pile EL 3 de 
Saclay, nous avons étudié la texture cristal- 
line d’une barre d’uranium brute de filage en 
phase «. Cette barre a été obtenue par filage 
a 600° C d’une billette de 90 mm de diamétre, 
le diamétre final étant 27,8 mm. 


2. Dispositif expérimental 


L’échantillon est un cylindre découpé dans la 
barre de maniére que son axe soit perpendicu- 
laire & axe de filage de la barre (fig. 1). Sur une 
face plane du cylindre, la direction de filage a 
été repérée. Ce cylindre a 19 mm de hauteur et 


Axe de Ja barre 


Axe de 
léchantillon 


Fig. 1. Schéma du mode de prélévement de l’échan- 


tillon dans la barre filée. 
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de diamétre (volume 5,4 cm3). Le diamétre 
moyen des grains étant de l’ordre de 100 uy, 
Véchantillon en contient quelques millions. II 
est placé dans un faisceau monochromatique 
de neutrons (A=1,143 A), l’axe du cylindre 
coincidant avec celui du faisceau (fig. 2). Si ’on 
pouvait placer en P normalement a l’axe du 
faisceau (fig. 2) une plaque sensible aux neutrons, 
celle-ci enregistrerait une série d’anneaux de 
Debye-Scherrer de densité non uniforme. En 
fait, le détecteur de neutrons dont on dispose 
est un compteur a BFs3 qui décrit dans le plan 
horizontal un cercle centré sur |’échantillon. 
On réalise néanmoins l’équivalent du dispositif 
de la plaque sensible en laissant fixe le compteur 
au maximum dintensité de la réflexion de 
Bragg étudiée, et en enregistrant les variations 
dintensité quand on fait tourner l’échantillon 
autour de son axe, c’est & dire quand on fait 
défiler devant le compteur les génératrices du 
cOne de Debye-Scherrer correspondant a la 
réflexion étudiée. 


3. Résultats 


Un spectre de diffraction ordinaire a tout 
d’abord été réalisé au moyen d’un échantillon 
ayant subi 4 traitements f suivis d’un traite- 
ment «. Nous nous sommes ensuite assurés que 
Véchantillon ne donnait qu’une faible orien- 
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Fig. 2. Principe du dispositif expérimental. 


tation préférentielle. Ce spectre (fig. 3) montre 
que les raies intéressantes a étudier sont les 
raies nettement séparées (111), (112), (131) 
assez intenses, et éventuellement (020), (022) 
et (130). Le paquet de raies (110), (021), (002) ne 
peut malheureusement donner d’indication bien 
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Fig. 3. Spectre d’un échantillon duranium filé traité ne présentant qu’une faible orientation préférentielle. 
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Variation d’intensité des cinq raies étudiées a: (020), b: 


0 


-180° -90 0 90° 180 


(111), ¢: (130), d: (112), e: (131) en 


fonction de l’angle de la direction de l’axe de fibre avec la verticale. 


précise. Nous nous en sommes néanmoins servis 
en fin d’étude pour confirmer les résultats 
obtenus. 
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Fig. 5. Diagramme qu’aurait enregistré une plaque 
sensible aux neutrons, placée normalement au 
faisceau incident en P fig. 2. 


La figure 4 représente la variation d’intensité 
des cing raies étudiées en fonction de langle fait 
par la direction de l’axe de fibre avec la verti- 
cale. La figure 5 donne une image des résultats 
qu’aurait enregistrés une plaque placée en P 
(fig. 2). La construction d’ Ewald (fig. 6) montre 
que l’angle gnxi entre la direction FF’ de laxe 
de fibre et le rayon qui joint le centre du cliché 
au centre d’une tache T de l’anneau (hkl) est 


relié & Vangle gonx: entre l’axe de fibre et le 
— 
vecteur OMyx; du réseau réciproque ayant 


donné naissance a la tache T par la relation: 


COS Ponki = COS Pnkl + COS Onrl- 
En ne considérant d’abord que les maxima les 
plus intenses des réflexions étudiées, on obtient 


les résultats du tableau 1. Connaissant les 


TABLEAU 1 


hkl Pret Onxr Ponkt 
020 678 iE 1S 67,5° 
WUE DAE 14° 30’ 26° 
1B 52° 18° 39’ D455" 
130 2655. 20° 52’ oon 
131 30a 21° 50’ Sie 
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Anneau hk] 


e2 
Fig. 6. 


Construction d’Ewald. 


a 
angles qonxni des vecteurs OMnxz; du réseau 
réciproque avec la direction de l’axe de fibre, 
cette direction est, dans ce réseau, celle de la 


génératrice commune a tous les cones d’axes 
> 
OMnpx: d’origine O et de demi-angles au sommet 


gonkt. La figure 7 représente la projection 
stéréographique de cette construction, le pdle 
de projection étant le nceud (110). On voit que 
la direction de l’axe de fibre est trés voisine de 
la direction [110] du réseau réciproque. 
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=~ 130 


Fig. 7. Projection stéréographique en prenant le 
noeud (110) comme pédle de projection. 
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On peut alors calculer les angles gonzi que 
fait la direction [110] du réseau réciproque avec 


les différentes directions OM de la forme 
{hkl} et vérifier que les maxima secondaires des 
anneaux de Debye-Scherrer qui n’ont pas 
encore été pris en compte correspondent bien 


X 


& ces nouveaux angles. Le tableau 2 donne le 


TABLEAU 2 


hkl Ponnicale PonkleXP 
ie 

020 34° 
64° 67,5° 

111 pete 262 
ye 52° 

112 46° 54,5° 
64,5° 

130 29,5° 25 
81,5° 76,5° 

131 34° 37° 
So 78° 


résultat obtenu. La proximité des deux angles 
go calculés pour les réflexions de la forme {112} 
explique en particulier la largeur des pics 
trouvés dans le cas de la réflexion (112) (fig. 4d), 
et pourquoi le cercle (112) de la figure 7 ne passe 
pas par le point d intersection des quatre 
autres. Seuls les maxima secondaires de la 
réflexion faible (020) ne s’expliquent pas par 
Vexistence d’un axe de fibre dirigé suivant [110]. 
Par contre, on peut, avec cette hypothése, 
calculer les angles auxquels il faut placer 
Véchantillon pour obtenir les maxima d’intensité 
des raies (110), (021), (002) (0, 68, 90 degrés 
respectivement). La figure 8 donne les formes 
du paquet des trois raies pour ces angles, le 
spectre étant obtenu en couplage 6/26. On voit 
que pour chacun de ces angles, seule la raie 
correspondante apparait, les deux autres ayant 
pratiquement disparu. En outre, comme on 
pouvait le prévoir, le maximum de la réflexion 
(110) est particuliérement intense. 
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Fig. 8. Forme des trois raies (110), (002) et (021) 


pour les valeurs 0°, 68° et 90° des angles auxquels il 
faut placer respectivement |’échantillon pour obtenir 
les maxima d’intensité de raie. 


4. Conclusion 


Jetter et McHargue 8) ont étudié aux rayons 
X la texture de barres d’uranium obtenues par 
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filage & 500° C. Leurs résultats, représentés en 
figures de pdles inverses montrent l’existence de 
deux directions privilégiées, l’une voisine de 
[110] du réseau réciproque, lautre de [010]. 
D’aprés leurs résultats, la proportion de cristaux 
appartenant a la texture [110] croit au profit 
de l’autre quand le rapport de réduction par 
extrusion augmente. En accord avec une étude 
de Calnan et Clews%), Harris!) étudie la 
texture de barres d’uranium laminées et montre 
que le méme phénoméne se produit quand le 
rapport de réduction par laminage et la tempé- 
rature a laquelle s’effectue le laminage aug- 
mentent. La barre ayant été filée a 600°C 
avec un rapport de réduction de 10,5, le fait 
que nous n’ayons pas décelé de texture suivant 
[010] est compatible avec ces résultats. 
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Les pouvoirs thermo-électriques de l’uranium et du 
plutonium ont été mesurés par une méthode différen- 
tielle. En accord avec les mesures précédentes de 
Waldron et Lee, on a trouvé pour le plutonium un 
pouvoir thermo-électrique élevé pour la phase « ainsi 
que les phases f et y, faible pour les phases de haute 
température 6 et «, surtout 6. Ce dernier résultat 
exclut Vexplication du _ coefficient de dilatation 
négatif en phase 6 avancée par Varley. Les variations 
observées du pouvoir thermo-électrique avec la 
température sont un peu différentes de celles rappor- 
tées par Waldron et Lee. 


The thermoelectric power of uranium and plutonium 
have been measured by a differential method. In 
agreement with previous measurements by Waldron 
and Lee, a high value of thermoelectric power has 
been found for the « and f phases of plutonium, and 


1. Introduction 


On sait que dans l’approximation des bandes 
le pouvoir thermo-électrique est directement lié 
a la dérivée en fonction de l’énergie du nombre 
délectrons au voisinage du niveau de Fermi. 

Dans un travail récent !), Varley propose une 
explication du coefficient de dilatation négatif 
observé pour le Plutonium en phase 6, expli- 
cation qui implique une dérivée fortement posi- 
tive du nombre d’électrons en fonction de 
Vénergie, au niveau de Fermi. 

Il était intéressant de confronter cette hypo- 
thése avec l’expérience. Comme nous disposions 
de peu de métal, nous avons utilisé le montage 
décrit plus loin que nous a suggéré Crussard. 
Waldron et Lee ont publié a la Conférence de 
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a low value for the 6 and e phases, particularly the 0. 
This last result rules out Varley’s explanation for the 
negative thermal expansion coefficient of the B phase. 
The observed variation of the thermoelectric power 
with temperature is slightly different from that 
reported by Waldron and Lee. 


Die thermoelektrische Kraft des Urans und des 
Plutoniums wurden mittels einer differentiellen 
Methode gemessen. In Einklang mit den friiheren 
Messungen von Waldron und Lee wurde eine hohe 
Thermokraft fir die « und f Phasen des Plutoniums 
gefunden, und eine niedrige fiir die 6 und ¢ Phasen, 
besonders die erstere. Dieses Ergebnis _ schliesst 
Varleys Erklarung des negativen Ausdehnungs- 
koeffizienten von 6-Plutonium aus. Der Temperatur- 
gang der Thermokraft war etwas anders als der von 
Waldron und Lee berichtete. 


Geneve 2) une série de résultats sur le pouvoir 
thermo-électrique du Plutonium pur que nous 
comparons aux notres. 


2. Principe de la mesure et dispositif expéri- 
mental 


Ktant donné les quantités de métal dont 
nous disposions, nous avons été conduits a 
adopter pour nos expériences le montage 
suivant 3) schématisé sur la figure 1. 

Le fil F dont on veut mesurer le pouvoir 
thermo-électrique, long d’une dizaine de centi- 
métres, est soudé a chacune de ses extrémités a 
2 fils de 2 métaux étalons M, et Me. Les quatre 
extrémités libres de ces fils sont a leur tour 
soudées a des fils de cuivre, ces quatre derniéres 
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Schéma des connexions. 


Fig. 1. 


soudures étant maintenues a la température 
constante de 0° C dans un vase Dewar contenant 
un mélange d’eau et de glace. 

L’ensemble du montage est placé dans un 
four de telle facon que les deux extrémités du 
fil F soient a des températures différentes. 

On enregistre en continu, sur deux potentio- 
métres les différences de potentiel ¢) et dz aux 
bornes des circuits M,;FM, et MokMs. 


— 


Ty 
¢i = J (Sm,—Sr)dT, ¢2 = J (Su,—Sr)dT 
Ty 


| 
8 


(Ta et Ty étant les températures des deux 
extrémités du fil). 

Si lon admet des variations linéaires du 
pouvoir thermo-électrique S dans les inter- 
valles de température considérés, ceci peut 
s’écrire: 


f1=(To—Ta)(Su,—Sr),  $2= (T—Ta) (Sm, — Sr) 


Sp= Su, + (Su, —Sm,) noe 


3. Détail du montage 


L’ensemble du montage décrit plus haut est 
placé dans un bloc de cuivre démontable dont 
une moitié a été représentée sur la figure 2. 
Il a pour réle de régulariser le gradient entre les 
deux extrémités du fil. Les fils sont isolés et 
protégés par des gaines de zircone pure frittée. 

L’ensemble de cet équipage est maintenu 
dans l’axe d’un tube laboratoire horizontal en 
silice, qui constitue l’enceinte & vide. Le four 
est mobile par rapport au tube, ce qui permet 
le cas échéant de modifier la différence (7'q—7'p) 
et de faire en sorte qu’elle soit de l’ordre de 30°. 


Les soudures ont été faites sur une machine 
& souder par points qui a été montée dans la 
méme boite 4 gants que le four. 

Les conditions expérimentales de montage en 
boite & gants nous ont obligés & couper nos fils 
étalons; aux jonctions complémentaires ainsi 
eréées entre les parties fixes des fils et celles qui 
sont soudées sur l’échantillon, les contacts sont 
assurés par serrage. 

Nous avons vérifié que ces jonctions créent, 
dans un circuit constitué par un métal unique, 
entre deux soudures froides a 0° C une différence 
de potentiel thermo-électrique qui reste con- 
stamment inférieure 4 10-2 millivolt et qui ne 
subit aucune variation brutale. 

Nous avons utilisé comme métaux étalons le 
chromel et l’alumel d’une part, le platine et 
Valliage platine-rhodium a 10 % d’autre part. 
C’est pour améliorer la sensibilité de la méthode 
(qui est comme nous allons le voir fonction des 
étalons choisis) et pour éviter les diffusions a 
état solide que nous avons été conduits a 
adopter les premiers pour l’uranium sur lequel 
nous avons fait une série d’expériences pré- 
liminaires, les seconds pour le plutonium. 

Tl nous a fallu procéder a un étalonnage pré- 
alable de l’alumel par rapport au platine pour 
avoir son pouvoir thermo-électrique. 

Les valeurs du pouvoir thermo-électrique 
absolu du platine que nous avons adoptées sont 
celles données par Borelius 4). 


Bouchon 


Fig. 2. 


Détail du montage. 
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4. Sensibilité de la méthode 


L’inconvénient de cette méthode  réside 
évidemment dans le nombre important des 
jonctions, ce qui accumule les sources d’erreur. 
En outre, il est clair qu’une hétérogénéité 
méme faible de ’échantillon peut provoquer une 
erreur importante si le fil présente un point 
chaud dont la température est notablement 
supérieure a celle de son extrémité chaude. 

Cest pour éviter cet inconvénient que nous 
avons placé l’échantillon dans un bloc de cuivre. 

Elle présente par contre l’avantage de donner 
directement Vindication physique intéressante, 
c’est-a-dire le pouvoir thermo-électrique. 

Il est évident que les conditions optima de 
sensibilité sont réunies lorsque le pouvoir 
thermo-électrique de |’échantillon est intermédi- 
aire entre ceux des deux métaux étalons. II est 


alors indépendent du choix des étalons. Par. 


contre si |Sp—S,| > |Sm,—Sm,], la sensibilité 
de la mesure devient, toutes choses égales par 
ailleurs, trés mauvaise. On a en effet, si l’on 
admet sur chaque mesure ¢; etde une erreur 7 
indépendante de la mesure: 


|Su, —Sn,| ae |Su, —Sr| 


et Z MD ees ie 2 
AS = « Pa) 


Ta—T»s 


L’uranium a un pouvoir thermo-électrique 
fort et positif. Platine et platine rhodié ont, par 
contre, des pouvoirs thermo-électriques voisins 
et négatifs: nous leur avons préféré l’alumel et 
le chromel dont les pouvoirs thermo-électriques 
encadrent jusqu’a 600°C celui de Vuranium. 


5. Mesures 
Ol 


Nous avons fait d’une part sur un méme échan- 
tillon une série de mesures en phase « sans 
dépasser la température de transformation 
(665°) et d’autre part plusieurs cycles de 
mesures avec passage en phase y. Dans ces 
derniéres expériences, la diffusion des composés 
du chromel et de l’alumel dans uranium étant 
importante, nous avons, aprés chaque essai 
coupé & chaque extrémité les zones de soudure. 

Les résultats, portés sur la figure 3, sont 
légérement différents au chauffage et au 


URANIUM 


refroidissement et se situent de part et d’autre 
des valeurs données par Dahl et Van Dusen 5). 

Cette différence est sans doute due & une 
hétérogénéité du métal, méme faible, associée 
a un gradient de température variable d’une 
extrémité a l’autre. I] est 4 noter que le fil était 
constitué de trés gros cristaux. Elle disparait en 
phase y; on ne l’observe pas non plus lorsque 
Yon remplace l’uranium par un métal de 
structure cubique, le platine par example. Elle 
est sans doute due au fait que le pouvoir thermo- 
électrique de l’uranium « doit étre anisotrope. 


5.2. PLUTONIUM 


Pour le plutonium, les données de Waldron 
et Lee montraient que l’on pouvait espérer une 
précision raisonnable en adoptant le platine et 
le platine rhodié comme métaux de référence, 
surtout en phase 6 et en phase «. En outre une 
diffusion importante était 4 craindre en-dessous 
de 320° avec l’alumel et le chromel; nous ne 
pouvions d’autant moins l’admettre qu’il nous 
a fallu opérer sur un fil plus court. 

Le fil qui nous a servi d’échantillon provenait 
d’un lingot de métal de haute pureté que nous 
avions filé sous la forme d’un fil de 1 mm de 
diametre. Ce métal présentait globalement 
moins de 250 ppm d’impuretés au cours d’un 
dosage initial. 

Nous avons fait en tout quatre cycles de 
mesures, les deux premiers de 20 a 600°C, le 
troisiéme de 20 & 316° de facon a pouvoir faire 
des mesures de refroidissement en phases y et £, 


Dahl et Van Dusen 
SS 


30 F 


aN 
‘y 
20 - he 


50 
v y 
[ons v/°K id Chauhage 
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ee ent 
7\hefroidissement (cyclage en a) 
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Pouvoir thermo-électrique de Vuranium. 
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Fig. 4. Pouvoir thermo-électrique du plutonium. 


le quatriéme 4 nouveau avec passage en 0. De 
méme que pour l’uranium, nous avons, entre 
le premier et le deuxiéme essai, refait les soudures 
aprés avoir éliminé 4 chaque extrémité du fil 
environ 3 mm de métal. 

Les résultats sont reproduits sur la figure 4. 

Nous avons trouvé pour toutes les phases, 
sauf « ou le pouvoir thermo-électrique est 
assez élevé, des pouvoirs thermo-électriques 
faibles, décroissant faiblement lorsque la tempé- 
rature augmente. 

Nous pouvons en outre penser raisonnable- 
ment que nous commettons sur nos mesures une 
erreur inférieure 4 0,5 ou 1 wV par degré K. 
En effet, si nous revenons a la formule (2), dans 
le cas le plus défavorable (plutonium «) 

(Spt —Spu) + (Sptrn —Spu) 

Spt —Spirn 

est de l’ordre de 4, (Ta—T7) de Vordre de 30 
a 40°, 7 d’environ 10 wV. Cette derniére valeur 
correspond a l’écart entre nos mesures et celles 
de Dahl et Van Dusen (figure 2). Elle est de 
toutes fagons d’un ordre de grandeur raisonnable 
pour un circuit différentiel MFM. 

Cette erreur (0,5 a 1 wV/° K) ne doit d’ailleurs 
varier que faiblement en fonction du temps, 


comme le montre la dispersion faible de nos 
résultats au cours d’un méme essai. 

Il faut noter que nous avons obtenu des 
résultats différents entre le premier essai et les 
suivants en phases «, 6, y qui sont anisotropes. 
Les mesures sont par contre trés reproductibles 
en phases 6 et e. 

Ceci est sans doute a relier a la possibilité 
dune texture qui aurait varié au cours des 
cycles ou peut-étre a la diffusion éventuelle du 
platine qui a donné des effets plus sensibles 
dans les phases de basse température 4 structures 
complexes, propres au plutonium. 

Dans lensemble, les valeurs, trouvées au 
premier essai pour chaque phase dans son 
domaine d’existence sont en bon accord avec 
celles de Waldron et Lee. Nous avons toutefois 
observé une loi de variation du pouvoir thermo- 
électrique en fonction de la température 
légérement différente, surtout en phase «. 

Notons enfin que nous avons observé au 
refroidissement, en dessous de 280°C, une 
légére inflexion de la courbe en phase 6, qui est 
a relier au changement de pente observé dans 
les courbes dilatométriques. On attribue ce phé- 
noméne a un début de transformation 6 — y 6), 
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6. Interprétation des résultats pour le 
plutonium 


A température 7’, supérieure 4 l’ambiante, et 
lorsque les électrons de conductibilité sont 
justiciables de lapproximation des _bandes, 
on sait que le pouvoir thermo-électrique S est 
lié a la conductibilité électrique op par la 
relation *): 


_ «wk Fd log 00(#) 
See ag 


k étant la constante de Boltzmann, e la charge 
électronique; et Hy l’énergie au niveau de Fermi. 

Les fortes valeurs de la susceptibilité de Pauli 
du métal montrent qu’il posséde une bande 
étroite a forte densité d’états (Hr) au niveau 
de Fermi. La conductibilité est limitée alors par 
les transitions des électrons de conductibilité 
vers les états vides de cette bande étroite. Dans 
les métaux de transition usuels, ces transitions 
seraient dues aux diffusions par les phonons. 
On aurait, au dessus de l’ambiante 7) 


= [(#)T m(#) (4) 


ou la fonction f varie lentement avec |’énergie 
EH. D’ot: 


are (5) 


On sait que contrairement a |’équation (4) la 
résistivité du plutonium n’est pas du tout 
proportionnelle 4 la température absolue 7’: au 
dessus de 100° K elle varie peu avec 7' et cela 
dans toutes les phases 8). Ceci suggére que, dans 
toutes les phases, la diffusion par phonons est 
faible et qu’une diffusion par des moments 
atomiques localisés prédomine. Ces moments 
seraient désorientés au-dessus de 100° K en 
phase «, ainsi que dans les domaines de stabilité 
des autres phases. Un comportement analogue 
de la résistivité est observé dans certains autres 
métaux (manganése, métaux des terres rares, 
etc.9)). Aux températures bien supérieures a 
celles ot des moments atomiques cessent d’étre 
couplés, la fréquence de choc 1/t devient 
indépendante de la température. Elle est 


proportionnelle 4 la vitesse des électrons de 
conductibilité et a 11(£). 
D’ot: 
I 1 


= Cx 
00 t Nest 


v(#) m(£) % 
Neu (By) (6) 


Le nombre effectif Neg d’électrons de conducti- 
bilité et leur vitesse v varient peu avec |’énergie; 
on retrouve donc une formule trés voisine de (5). 
Pour des électrons libres d’énergie minima Ko 
par exemple, 


1 
Nett Vv a (H—Ep), 


_ eR dm 1 
aa Sear ~~ Ey a): (7) 


Pour un nombre raisonnable d’électrons de 
conductibilité par atome, leur énergie (Hy — Eo) 
sera une fraction appréciable d’unité atomique. 

Les équations (5) et (7) sont alors pratique- 
ment équivalentes. Elles donnent pour -ai 
les valeurs assez faibles portées au tableau 
suivant: 


GEOR) 325 | 435 | 535 | 635 | 800 
Phase a B y ) € 
S(uV/° K) 12,5 8,2 9242.6 3,8 
dni 
pec ire (Va) 40 19 19,5 | 4,3 5 


Ces valeurs sont positives, comme dans 
Vuranium 1°), mais environ trois fois plus faibles. 
Elles sont done probablement trop faibles pour 
que la théorie de Varley s’applique au Pu. 

D’autre part, dans chaque phase, sauf peut- 
étre «, S n’est pas proportionnel a 7’. Les 


dn , ° 
aie déduites 


précédentes varient done un peu dans chaque 
phase avec la température. Cet effet est peut- 
étre lié & une forte courbure de 1(#) au voisi- 
nage du niveau de Fermi. Une telle courbure 
expliquerait également les faibles variations 
avec la température observées pour la résistivité 
et la susceptibilité. 


valeurs de 


des opérations 
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7. Conclusion 


Méme si l’on tient compte des hypothéses que 
suggére le comportement électronique anormal 
du plutonium, les valeurs que nous avons 
trouvées pour le pouvoir thermo-électrique du 
plutonium restent trop faibles pour pouvoir, 
semble-t-il, s’accorder avec la théorie de Varley. 
Waldron et Lee étaient d’ailleurs parvenus a la 
méme conclusion dans leur communication de 
Genéve. 

On se serait de toutes fagons attendu a 
un pouvoir thermo-électrique particuli¢rement 
élevé en phase 6. Il semble done que les raisons 
véritables du comportement anormal du pluto- 
nium 6 restent a trouver. 


Nous tenons a remercier M. Jacques Friedel 
pour Taide qu'il nous a apportée au cours 
de notre travail, tout particuliérement dans 


Vinterprétation des résultats, et & lui exprimer 
ici notre profonde reconnaissance. 
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The U-Mo equilibrium diagram has been determined 
up to 19 wt % Mo and below 900° C. The eutectoid 
transformation 6 =x«+y has been located at 
639 + 5°C. The limit of the 6 + y field has been 
placed at 4.5 wt % Mo at the eutectoid temperature. 
The y phase undergoes a eutectoid transformation at 
10.5 wt % Mo and 565 + 5°C, to « + 6. The tetra- 
gonal structure of the 6 phase was confirmed, with 
co = 9.854 A, ap = 3.427 A, c/a = 2.876. It is probable 
that 6 phase is formed by a congruent transformation, 
rather than by a peritectoid as reported by earlier 
workers. 


Le diagramme d’équilibre U-Mo a été déterminé 
jusqu’&a une teneur en Mo en poids de 19 % et une 
température inférieure a 900°C. La transformation 
eutectoide B = (« + y) a été située a 639 + 5° C. La 
limite du champ (f+ y) a été situé a 4,5 % en 
poids de Mo a la température du palier eutectoide. 
La phase y subit une transformation eutectoide en 


1. Introduction 


The U-Mo phase diagram has been the subject 
of several earlier investigations. Those prior to 
1958 were summarized by F. A. Rough and 
A. A. Bauer!) and were also reviewed by A. J. 
Carrea et al.2). The diagram of S. T. Kono- 
beevsky et al.?) appeared in 1958 and differs in 
several respects from the earlier work, notably 
in that the latter authors show the 6 (UzMo) 
phase to be formed by a congruent transfor- 
mation from the y phase, rather than the 
peritectoid transformation reported by earlier 
workers. 


2. Experimental Technique 

In the present investigation alloy buttons 
were arc melted, using a tungsten electrode and 
water cooled copper hearth, from high purity 
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(x + 6) pour 10,5 % en poids de Mo et 565° + 5° C. 
La structure tétragonale de la phase 6 a été confirmée 
avec les paramétres suivants: co = 9,854 4A, ao 
3,427 A, c/a = 2.876. Il est probable que la phase 6 
est formée par une transformation congruente plutdt 
que péritectoide comme d’autres chercheurs l’ont 
proposé antérieurement. 


Das U-Mo Gleichgewichtsdiagramm wurde bis zu 19 
Gew % Mo und bis zu 900° C ermittelt. Die eutektoide 
Umwandlung B =« + y liegt demnach bei 639 + 
5° C. Bei der eutektoiden Temperatur reicht das 
B+ y-Feld bis zu 4,5 Gew% Mo. Die y-Phase 
unterliegt bei 565 + 5°C eimer eutektoiden Um- 
wandlung bei 10,5 Gew % Mo in « + 0. Die tetra- 
gonale Struktur der 6-Phase konnte bestatigt werden, 
mit co = 9,854 A, a = 3,427 A und c/a = 2,876. 
Wahrscheinlich wird die 6-Phase kongruent und nicht, 
wie friher berichtet wurde, durch peritektoide 
Umwandlung gebildet. 


electrolytic uranium described by Blumenthal 
and Noland 4) and commercial purity molyb- 
denum. The buttons were in most cases sealed 
in a zirconium alloy sandwich and hot rolled 
at 900° C. The resulting slabs were homogenized 
in the y range and cut up for X-ray, metallo- 
graphic and chemical analysis specimens. All 
alloys were chemically analyzed except a few, 
which were intended exclusively for use in 2 
phase regions. Specimens for X-ray diffraction 
were small needles, roughly 0.4 cm long and 
0.08 cm in diameter. These were wrapped in 
Mo foil, sealed in evacuated Vycor capsules for 
heat treatment, and water quenched. Metallo- 
graphic specimens were processed in the same 
manner. 

All X-ray diffraction work was done at room 
temperature using Debye—Scherrer cameras. X- 
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ray needles were cleaned electrolytically in an 
acetic-perchloric acid bath. Since high precision 
was necessary to determine the lattice para- 
meter of the retained gamma phase, the 
Nelson—Riley extrapolation was applied. 
Metallographic specimens were etched electro- 
lytically in a bath of phosphoric acid, ethylene 
glycol and ethyl alcohol. 


3. The Phase Diagram 
Fig. 1 shows the graph of lattice parameter 


in A vs. composition in at % Mo. A straight line 
having the equation ap = 3.4808 — 0.00314 at % 


U 4 8 l2 16 20 


Big. et 


end of the graph, a distinct inflection is apparent 
corresponding to the y/(y + Mo) phase boundary. 

The phase diagram resulting from this 
investigation is shown in figs. 2 and 3. The 
points shown represent 73 of the 269 compo- 
sitions and temperatures investigated. Lines 
AE and AD were determined by metallographic 
examination of specimens which were quenched 
from transformation temperatures above and 
below the line. This technique, known as the 
method of temperature brackets, is preferred 
in the high uranium end of the diagram (up to 
6 wt % Mo) because equilibrium is readily 


= ee ale 


700 800-7 39005 
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Lattice parameter curve for U-Mo system, measured at room temperature on specimens quenched 


from the y range. 


Mo may be fitted to the data points by the 
method of least squares. The intercept of this 
line on the ordinate, 3.4808 A, is the extra- 
polated room temperature lattice parameter for 
y-U, subject to the assumption that the lattice 
parameter curve remains straight in the high U 
end where y cannot be retained for room 
temperature measurement f. At the high Mo 

+ The lowest Mo content in fig. 1 is 12.7 at % Mo. 
A careful examination of the diffraction pattern 
reveals that the structure of this alloy is retained y, 
which is in agreement with the data reported by 
Ivanov and Badajeva |!) who retained y in an alloy 
of 13 at % Mo, but found a tetragonal phase in a 
12 at % Mo alloy. In work at Harwell, Hills et al.12) 
inferred by metallographic examination that a tetra- 


attained in times ranging up to 4 days, and 
because a small amount of f in y or y in f is 
easily detected. Extrapolation of line AE back 
to pure U gives the value 775° C for the f to y 
transformation in close agreement with the 
value 774.8 -+ 1.6° C reported by Blumenthal 5), 
Figs. 4a and 4b show the microstructure of a 
3 wt % Mo alloy quenched from 665° C in the 
B+y field to consist of 6 particles distributed in 


gonal phase is formed by quenching a 15 at % Mo 
alloy. In the present investigation, as well as in the 
work by Ivanov and Badajeva, it is indicated that 
a b.c. cubic y is retained in the 15 at % Mo alloy. 
This apparent discrepancy may be due to a difference 
in quenching rate. 
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an x’ matrix (y matrix at 665° C). The f particles 
characteristically form along the y grain 
boundaries. The «’ phase, also called martensitic 
o and strained «, appears in at least eight 
uranium-base systems. It. forms on quenching 
from the gamma phase when there is sufficient 
solute to suppress the normal y>fi—>« 
sequence but insufficient solute to permit 
retention of metastable y to room temperature. 
In the present investigation, the lower limit of 
Mo content necessary for martensitic « for- 
mation is placed at 2 wt %. A 2 wt% Mo 
specimen quenched from 696°C exhibited a 
mixture of martensitic « and the Widman- 
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statten « which is obtained at lower Mo 
contents. The maximum Mo content permissible 
for martensitic « formation lies between 4.5 and 
5.5 wt %. The formation of the «’ phase has 
been studied by Harding and Waldron §). 
The eutectoid horizontal CDE was located at 
639+ 5° on the basis of both metallographic 
and X-ray diffraction data. X-ray diffraction 
needle specimens and metallographic specimens 
were capsulated together, given a 1 hour soak 
in the gamma range, then transferred quickly 
to a furnace set at a predetermined temperature. 
As a typical example, specimens of U+1.33 
wt % Mo alloy were held for 3 days at 645° C, 
water quenched, and examined as soon as 
possible after quenching. The diffraction film 
exhibited the pattern characteristic of £ 
uranium, and the microstructure showed y 
particles in a 6 matrix, fig. 5. Successively lower 
temperatures were used until specimens quench- 
ed from 633° C exhibited only an « diffraction 
pattern. To guard against the possibility of 
error due to supercooling of £, a refinement of 
the bracketing method was used. Specimens of a 
0.56 wt % Mo alloy were heat treated at 625° C 
to produce an «+-y structure, shown in fig. 6a. 
Other specimens were heat treated at 645° C 
to produce a f or $+y structure, shown in 
fig. 6b. Pairs of specimens were then capsulated 
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Fig. 4a. Retained 6 network around Martensitic «, 

in U + 3 wt % Mo alloy after quenching from 665° C. 

At 665° C structure was B + y. Polarised light. Etchant: 
phosphoric acid. 200 x. 


Fig. 40. 
as above but bright field illumination. 
phosphoric acid. 200 x. 


Retained 6 network around Martensitic «, 
Etchant: 


Big. 5: 


together for heat treatment at 635 and 645° C. 
It was found that a prior « + y structure 
transformed to f at 645° C and a prior # struc- 
ture transformed to « + y at 635°C. The equilib- 
rium transformation temperature was thereby 
bracketed between these two temperatures. An 
evaluation of all the pertinent data showed 
639° C to be the most probable temperature. 


«’ (formerly y) particles in 6 matrix after quenching a 1.33 wt % Mo + U alloy from 645° C. 
Bright field. Etchant: phosphoric acid. 


400 x. 


The y/(y+«) line, EG in fig. 2, was studied 
intensively. Up to 6 wt % Mo, the line can be 
readily determined by the method of tempera- 
ture brackets. For example, an alloy containing 
5.93 wt % Mo with a prior «+¥y structure 
transformed to y at 625° C, but the same alloy 
previously quenched from the y range trans- 
formed to an «+y structure at 621° C. For Mo 
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Fig. 6a. y particles in « matrix after quenching a 
U + 0.56 wt % Mo alloy from 625°C. Bright field. 
Etchant: phosphoric acid. 400 x. 


contents higher than 6 wt % the transformation 
times necessary to reach equilibrium became 
excessively long, and the method of lattice 
parameters was used instead. In the latter 
method, X-ray needles and metallographic 
controls were transformed in the «+y region. 
The lattice parameters of the y phase were meas- 
ured as accurately as possible, and from fig. 1 
was obtained the corresponding composition of 
the y phase. One advantage of this method over 
that of temperature brackets is that of saving 
time. A 10 wt % Mo alloy for example would 
require at least two months at 570° C to obtain 
metallographic evidence of transformation, but 
a 6 wt % Mo X-ray specimen can be trans- 
formed to «++ at 570° C in less than two weeks. 
Further, the lattice parameter of the y in the 
a+y alloy gives a derived composition directly 
on the y/(«+y) line rather than a point lying 
above or below the line. 

An interesting feature of the «+ field is the 
occurrence within it of a lamellar « + structure 
formed by rejection of « from a supersaturated 
y solid solution when the rejection occurs at a 
temperature that is approximately 15°C or 
more below the solvus line. An example is shown 
in fig. 7, shown by X-ray diffraction to be an 
o«+ystructure. The dark lamellae were identified 


Fig. 6b. Retained 6 structure after quenching a 
U + 0.56 wt % Mo alloy from 645° C. Bright field. 
Etchant: phosphoric acid. 400 x. 


as « by the observation that they darken on 
standing in air; the white matrix is unaffected 
at room temperature. The range of the lamellar 
structure is from 5 to 8 wt % Mo, and 571°C 
to 605° C. Strangely, the lamellar structure is 
not found when the rejection of « occurs 
immediately below the («+y)/y line; here the 
structure consists of spherical « particles in a 
y matrix. It is believed that the lamellar «+y 


IKE le 
U + 5.92 wt % Mo alloy from 581°C. Bright field. 


« lamellae in y matrix after quenching a 


Ktchant: phosphoric acid. 1600 x. 
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structure may have misled earlier workers into 
reporting a eutectoid in the range 5 to 8 wt % 
Mo. Therefore, the fact that proeutectoid « is 
rejected from solution below a sloping line (EG) 
is emphasized. On the other hand, the lamellar 
structure shown in fig. 6a was formed by the 
eutectoid transformation of i >«+y. Of the 
two lamellar structures observed to form in the 
x+y field, the one resulting from eutectoid 
decomposition of 6 consists of y lamellae in a 
matrix of «, whereas the « lamellae precipitated 
from supersaturated y occur in a matrix of 
residual y. 

The lower boundary of the «+y field is the 
line FGH, at the eutectoid temperature. An 
attempt was made to determine the temperature 
of line FGH directly by heat treating X-ray 
needles of a 12 wt % Mo alloy (selected because 
it lay close to the probable eutectoid compo- 
sition) at 555, 560 and 564° C for 1 month. The 
resulting patterns exhibited a mixture of y and 
6 lines, indicating incomplete transformation. 
A more successful technique involved the use 
of a 3 wt % Mo alloy to allow rapid trans- 
formation. A finely dispersed «+6 structure 
was first created and then this structure was 
reheated to successively higher temperatures 
until the 6 changed to y. In this method, X-ray 
needles were transformed to «+6 at 500°C. 
which is near the nose of the TTT curve 7), thus 
permitting rapid transformation. Diffraction 
patterns were taken to confirm the presence of 
6. Due to the relatively low Mo content of 
3 wt %, a strong « and weak 6 pattern was 
obtained. The success of this method depends 
upon the fact that several strong lines of the 
6 diffraction pattern lie in a gap where there 
are no lines of the « pattern, making identifi- 
cation of a small amount of 6 quite certain. 
Reheating the needles to 555°, 560°, 565° and 
570° C for 3 weeks caused the disappearance of 
the 6 lines in the 570° C sample, but not in the 
other three. Therefore, the line FGH is located 
at 565-4 5° C. 

The crystal structure of the 6 phase was 
found to be tetragonal, of the MoSiz type, as 
first reported by Halteman 8). In the present 


investigation, a needle shaped specimen con- 
taining 16 wt °%% Mo was heat treated at 536° C 
and a diffraction pattern taken with Cu 
radiation. The parameters are cp) = 9.854-+ 0.001, 
ay = 3.427-+ 0.0010, c/a=2.876. The micro- 
structure of the 6 phase is shown in fig. 8. No 
attempt was made to locate precisely the limits 
of the 6 phase, but the apparent single phase 
present in the 16 wt % alloy is in agreement 
with the data of Halteman. The data of 
Konobeevsky et al.3) indicate a range of 16.5 
to 17 wt % Mo for the 6 phase. 

The y/(y+ Mo) line was investigated by both 
metallography and X-ray diffraction from 
600°C to 900°C. The metallographic results, 
all of which are shown in fig. 2, bracket the 
line between rather wide limits at 600, 700 and 
800° C. In theory the line can be located by 
means of fig. 1 at 700, 800 and 900° C by the 
intersections of the three iso-parametric lines 
of the y+ Mo alloys with the sloping line of the 
single phase y alloys. However, the small 
variation of ao over this temperature range, the 
relatively small amount of X-ray data in the 
two-phase region and the limited precision of 
the data of fig. 1 all combine to limit the 
accuracy with which the phase boundary can 
be determined in this way. A rather detailed 
analysis of the data indicates that at the three 


Fig. 8. 6 phase formed at 532°C in alloy of U + 
16 wt % Mo. Polarised light. Etchant: acetic-chromic. 
400 x. 
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pertinent temperatures the y/(y+Mo) line lies 
between the following limits: 


700° C: 16.6-18.2 wt % Mo 
800° C: 16.7-18.3 wt % Mo 
900° C: 17.3-18.9 wt % Mo. 


The X-ray data points on fig. 2 are the mean 
values and the boundary has been drawn 
through these. 

The location of the y/(y+Mo) line was also 
confirmed by the method of vanishing lines. 
As an example, a pattern taken of an alloy of 
U+18.19 wt % Mo heat treated at 800°C 
showed a faint (110) line of the Mo pattern, but 
an alloy of U+ 16.32 wt 9% Mo showed no trace 
of a Mo pattern. 

The eutectoid horizontal passing through 
point J of fig. 2 is purely speculative, as is the 
peak of the 6 phase. It is believed that the 
y/(y-+Mo) line is known with sufficient accuracy 
to require a eutectoid at J although no data are 
available. 


4. Discussion 


There is general agreement that the lattice 
parameter of y-U is contracted by Mo additions, 
but there is some variation among the results 
of different investigators. Wilson and Rundle 9) 
presented data for five alloys. Their data have 
been extrapolated to a value of 3.474 A for 
y-U at room temperature. Pfeil 1°) has reported 
parameter data on two alloys, and these data 
lie fairly close to the line determined in the 
present investigation. Konobeevsky et al.°) have 
reported the equation 3.481—0.00333 wt 9% Mo= 
ao, but did not report the number of points on 
which this equation was based, or the purity 
of the uranium. The present investigation is 
believed to give the most reliable lattice para- 
meter data as it is based upon ten alloys made 
from high purity electrolytic uranium. 


The U-Mo diagram in figs. 2 and 3 differs from 
those of earlier workers principally in the lower 
temperature for the 6 — «+ y eutectoid (639° C) 
and in the open maximum formation of 6 rather 
than the peritectoid transformation. The major 
difference between the diagram of Konobeevsky 
et al.3) and the present work is that the former 
authors place the maximum Mo content of the 
B+y field at 6 at % Mo and 635° C compared 
with 11 at °%4 Mo and 639° C in the present work. 
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When f-uranium is transformed into the 
x-phase on cooling, each grain is broken up into 
numerous subgrains, each slightly disoriented 
from its neighbours 1-8), It has been suggested +) 
that these subgrains are formed because of the 
large volume change accompanying the trans- 
formation. It was proposed that the « phase, 
which is softer than the f phase at the trans- 
formation temperature, suffers distortion during 
the transformation; this distortion is a form of 
creep which, at the high temperature concerned 
results in the formation of ‘‘creep cells’? such 
as are formed in aluminium during creep. This 
process may also be regarded as a strain- 
induced polygonization. The faster the rate of 
cooling through the transformation tempera- 
ture, the faster is the effective rate of creep of 
the « grains, and the smaller should be the 
mean size of creep cells. This is in fact observed. 

As a test of the proposed mechanism, the 
orientation relationship between neighbouring 
pairs of subgrains has been determined. If 
subgrains are indeed formed by the thermally 
activated rearrangement of glissile dislocations, 
as in creep, then the sub-boundaries will 
normally consist of arrays of well-defined types 
of dislocations; if the active slip systems are 
known, the nature of the expected misorien- 
tations across the subboundaries can then be 
computed and compared with experiment. 
Creep cells in zinc have been examined ex- 
perimentally in this way 4), and a detailed 
interpretation of the results obtained was 
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given by Ball and Hirsch *) and by Ball ®). 

In the case of x-uranium, the slip systems 
active at high temperature have only recently 
been determined 7). At 600°C, the system 
predominating at room temperature, (010) [100], 
still persists but a new system with a much 
lower critical stress has appeared: (001) [100]. 
(At room temperature this only acts as a cross- 
slip system.) The (110) [110] system also 
operates at 600° C, but with a higher critical 
stress than either of the others. A fourth 
system, (021) [112], has such a high critical 
stress that for present purposes it can safely 
be ignored. 

If a sub-boundary were to consist entirely of 
edge dislocations of the (010) [100] system, it 
would be a simple tilt boundary and the 
rotation axis (i.e. the axis about which the lattice 
of one subgrain must be rotated to bring it into 
congruence with the lattice of its neighbour) 
will be normal to the slip direction and parallel 
to the slip plane, i.e. the axis will be [001]. 
Similarly, a tilt boundary consisting of (001) 
[100] edge dislocations would have a [010] tilt 
axis. A tilt boundary consisting of (110) [110] 
edge dislocations would, again, have [001] as 
rotation axis. 

An actual subgrain boundary formed in the 
way suggested above would normally contain a 
mixture of two different kinds of dislocations. 
The crystallography of this situation is dis- 
cussed by Ball and Hirsch >), The combination: 
(010) [100]+ (001) [100] corresponds to their 
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case 3; here the axis of rotation may take up 
any position in the plane normal to the common 
Burgers vector, i.e any position lying in the 
(100) plane. The combinations (010) [100]+ 
(110) [110], or (001) [100]+(110) [110], both 
correspond to Ball and Hirsch’s case 4, and 
here the rotation axis must be normal to both 
the Burgers vectors; i.e. it must be parallel to 
[001]. 

It follows that, if the present hypothesis as to 
the origin of subgrains in «-uranium is correct, 
all observed rotation axes relating to orientations 
of neighbouring pairs of subgrains should either 
coincide with [001] or else lie in the (100) 
plane (in other words, in the zone defined by 
[001] and [010]). The theoretical treatment 
also predicts the relationship of the rotation 
axis to the plane of the sub-boundary itself; 
however, this aspect was not examined in the 
present experiments. 

Three uranium rods were passed slowly 
through a fairly steep temperature gradient so 
that they transformed progressively from the 
6 into the « phase. Under these conditions 1) 
very large « grains, each decomposed into 
large subgrains, are formed. Several pairs of 
adjacent subgrains were examined in each 
specimen ; some of the sub-boundaries concerned 
were longitudinal, some inclined to the direction 
of the thermal gradient. In each case, a Laue 
photograph was made of the two subgrains 
together, by causing the X-ray beam to straddle 
the boundary. A specially-designed micro-beam 
camera ®) was used which allowed precise 
positioning of the beam. 

12 patterns were clear enough to interpret. 
In each case, the mean orientation of the « 
grain was determined from the Laue pattern, 
and the rotation axis was then determined from 
the relative positions of the various associated 
pairs of Laue spots, by an analytical method 
previously described 4). The results are shown 
in fig. 1, normalised into a quadrant of the 
standard stereographic protection of « uranium. 
Each circle represents the estimated experi- 
mental error in determining the corresponding 
rotation axis (the error depends primarily on 


[001] 


[100] 
Fig. 1. Observed rotation axes for 12 pairs of sub- 
grains in «-uranium samples, after B — oa trans- 
formation. 


the magnitude of the misorientation of the sub- 
grain pairs). Although some of the axes deviate 
slightly from the predicted [001]-[010] zone, 
all are close to it, and within the margin of error 
half of the axes actually lie in the expected zone. 

The results of these experiments, therefore, 
confirm the hypothesis that subgrains in B > « 
transformed uranium are due to a strain-induced 
form of polygonization of glissile dislocations in 
the «-uranium lattice. A corollary of this 
mechanism is that the subgrain configuration 
can be modified by applying an external stress 
during transformation; it is hoped to investigate 
this experimentally. 


We are indebted to Mr. B. R. Butcher for 
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Atomic Energy Authority for financial support. 


References 


1) R. W. Cahn, Acta Metallurgica 1 (1963) 176 

) P. Lehr, Comptes Rendus 241 (1955) 1043 

3) P. Lehr, Saclay (France) rapport CHA (1958) 

) R. W. Cahn, I. J. Bear and R. L. Bell, J. Inst. 

Metals 82 (1954) 481 

5) C. J. Ball and P. B. Hirsch, Phil. Mag. [vii] 46 
(1955) 1343 

8) C. J. Ball, Phil. Mag. [viii] 2 (1957) 977 

‘) HH. Hi. Chiswick, A. EX. Dwight, L.-T. Lloyd, 
M. V. Nevitt and S. T. Zegler, Second Geneva 
Conference (1958) 15/P/713 

8) H. M. Otte and R. W. Cahn, J. 
(1959) 463 


Sci. Instr. 36 


JOURNAL OF NUCLEAR MATERIALS 1, No. 1 (1960) 90-91, NORTH-HOLLAND PUBLISHING CO., AMSTERDAM 


CHANGEMENT DE PHASE ET MISE EN DESORDRE PAR IRRADIATION 
DU COMPOSE U2Mo 


Mme J. BLOCH 
Centre d’Etudes Nucléaires de Saclay, Gif-sur-Yvette (S & O), France 


Recu le 18 février 1960 


Nous avons étudié des plaquettes d’un alliage 
uranium-molybdéne ayant la composition de 
U2Mo. Aprés homogénéisation a 1050° C, 
Valliiage a subi un traitement de recuit de 
sept jours 4 500°C afin d’obtenir la phase y»’ 
quadratique ordonnée. 

Les échantillons ont été irradiés dans un flux 
de neutrons thermiques variant de 21019 a 
2x10! n/em2/sec et a une température infé- 
rieure a 65° C. Des moniteurs d’Al-Co ont servi 
a mesurer les flux intégrés regus par chaque 
échantillon; la gamme étudiée s’étend de 
1016 & 1,210!8 n/em?. Nous avons fait sur 
chaque échantillon, avant et apres irradiation, 
un diagramme de diffraction de rayons X sur 
un diffractometre 4 compteur !). 

Nous avons porté notre attention sur la 
disparition des raies d’ordre et sur le passage 
de la structure quadratique de la phase y’ a la 
structure cubique de la phase y. Cette phase est 
stable 4 haute température. La disparition par 
irradiation de cette phase y’ a déja été constatée 
par Bleiberg d’une part 2) et par Konobeevski 
d’autre part °) sur des alliages U-Mo présentant 
la structure «+y’. 


1. DISPARITION DES RAIES D’ORDRE 


Nous avons pour chaque échantillon évalué 
le degré dordre d’apres le rapport I'/I des 
intensités de la raie de surstructure (101) aprés 
et avant irradiation. 

D’aprés ces résultats 


préliminaires nous 
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avons tracé la courbe donnant J’/I en fonction 
du flux total ¢t, ce qui indique la vitesse 
de mise en désordre du composé considéré. 
(courbe 1) 
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CHANGEMENT DE PHASE 


Nous avons étudié chaque échantillon sous 
des angles de Bragg variant de 10 a 75°. A 
partir de ces diagrammes nous avons calculé la 
meilleure valeure des paramétres pour chaque 
échantillon irradié et tracé les courbes des 
paramétres en fonction du flux total (courbes 2: 
nous avons porté en abscisses a’ et 4c’ car 
dans la maille quadratique de y’ la valeur de 
c est voisine de 3a). 

Nous avons constaté une évolution des para- 
métres de la phase quadratique vers une 
valeur qui est celle du paramétre de la phase 
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cubique, traduisant donc un passage continu 
de la structure quadratique a la structure 
cubique. 


e 
344404} x Shi sede te eS 
34200 | 
34000 
33800 = 
33600 
33400 | 
3.3200 
3.3000 ; 
3.2800 3 

: 10” 108 ee 

Courbes 2 
TABLEAU | 

Flux total en ie = A Le! aheoe 
Shee Sie) (ye | tA (Ask) 
1,26 x 1016 3.4322 | 9,8592 | 3,2864 | 0,1458 
ib x L028) L00 3,4239 | 9,8747 | 3,2916 | 0,1323 
1,71 x 1016 | 42 | 3,4238 | 9,8930| 3,2976 | 0,1262 
dole xe 1026 17,6 | 3,3970 | 9,9621 | 3,3173 | 0,0797 
122° 10" 33 3,4045 |10,0345 | 3,3448 | 0,0597 
Seren Ile 14 3,3930 |10,1259 | 3,3753 | 0,0177 
4.8 x 10}? 3,3890 |10,0770 | 3,3590 | 0,0300 
Mg CANS 0 3,3798 0 


Of 


Wittels a constaté un phénoméne analogue 
dans des composés du type BaTiOs *). 

Tl n’y a, & aucun moment, coexistence des 
raies quadratiques et des raies cubiques comme 
on aurait pu s’y attendre si l’on pensait a un 
mécanisme de transformation qui soit celui 
d’une pointe de fission ou une petite région 
est chauffée au-dessus de la température de 
transition et trempée. C’est par exemple ce qui 
a été observé par Wittels 5) pour de la zircone 
contenant un peu d’uranium. 

La transformation, débutant a des flux aussi 
faibles que 1016 n/em? est achevée pour un flux 
voisin de 1,2 1018 n/cm?2. 
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In their recent paper “The sintering of 
uranium -dioxide in hydrogen at 1350° C” 
M. Bel and his co-authors !) raise a number of 
points which call for comment. 

The French investigators failed to detect any 
beneficial effect of departure from stoichiometry 
upon sintering behaviour in neutral atmos- 
pheres, in contradiction to our findings 2). In 
our studies we used no additions of binder or 
lubricant, and the pressed pellets were sintered 
on alumina trays; M. Bel et al. added 3 % 
camphor and 0.5 % stearine and probably 
sintered on molybdenum trays, although this 
latter point is not clear from the text. 

Our experience using a wide variety of organic 


binders with non-stoichiometric powders is that 
unless the binders are removed carefully, prefer- 
ably in an oxidising atmosphere, the non- 
stoichiometric powders can be reduced to the 
stoichiometric condition before appreciable sin- 
tering has occurred. In the absence of further 
data on the oxygen content of the pellets 
sintered in argon by French investigators it 
cannot be said therefore that they were truly 
studying the sintering behaviour of non-stoich- 
iometric oxides. A second factor which would 
also affect their results would be the use of 
molybdenum trays. We found that non-stoich- 
iometric oxide in contact with molybdenum 
sintered less readily than the same material 


TABLE | 
Powder characteristics > Sintering conditions Miss re 
‘SS ~ So ° 
tet Samed men a ee eto 
| 2 & Sleai wo aS eee come 
ga | 8 | 08 | PS asl] ee | SB | 88 | a3 
Origin € E ms ge | Re g © - | At here | ek =| 52 8 2 
a 8 8 | 8 aa) ¢ mosphere ts (°C i ves 
Be lege, Wendel PO ee cage 
OE Nema eet ee Eee tN FA = : 
oe a ee be ee le 
H 16 -_ 5.5 75 | 5 | 1400] Cracked —78*(2) | 10.4 92.9 | 89.1 
French ammonia 
J 4.5 — 6.3 10.25 OS eOogl 
6 2.18 4.6 —78 to —96 9.65 | 90.6 |110 
British 2 2.08 5.4 300 2 1410 | Hydrogen —96 9.07 | 80.0 | 68 
6 2.24 4.6 —96 8.24 | 76.4 | 78.7 
16 2.67 5.0 | —96 10.4 95.9 |108 


Dewpoint of —96°C = 2 x 10-5 mg/l; —78°C =5 x 10-4 mg/l 


* Assumed from the text to have been dried with magnesium perchlorate. 


eee 
{ ditors’ note: Dr. Deitmas’ reply to this letter will be published in this journal, Vol. 2, no. 2. 
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in contact with alumina. On the basis of the 
data reported, there is no evidence that an 
insensitivity of sintering behaviour towards 
oxygen content has been established. 

The second feature on which I wish to 
comment relates to the sintering of “‘active”’ 
uranium oxide powders in hydrogen. In the 
table 1 are assembled data selected from both 
the French and our own work. Comparison of 
the sintering behaviour of the various powders 
is complicated by their differing compacting 
behaviour and it is necessary to use some 
parameter which takes into account the vari- 
ation in green density. A reasonable, although 
not ideal parameter is the percentage decrease 
in porosity. Ignoring the difference in sintering 
eycles which reflect to the advantage of the 
French powders it will be seen that these 
powders, and the one British powder sintered 
in what we termed “incompletely dried hy- 
drogen”’ showed approximately the same per- 
centage decrease in porosity. Thus there seems 
to be no great difference between the two sets 
of results when rationalized. However when 
comparable British powders were sintered in 
“dry hydrogen’’, i.e. using phosphorus pentoxide 
as a dessicant, their sintering behaviour deter- 
iorated. This apparently large effect of a small 
variation in moisture content of the hydrogen 
upon sintering behaviour was emphasized in 
our original paper and is worthy of further 
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investigation. Judged on the basis of percentage 
decrease in porosity, the powder which sintered 
best in hydrogen was the British U30g which 
had the highest surface area and oxygen content 
of the powders examined. 

The third feature concerns the existence of 
the ‘‘over-sintering’’ phenomenon, i.e. the deter- 
ioration of density with increase of sintering 
temperature above some optimum value. To 
the best of my knowledge, this phenomenon 
was first reported for uranium oxide by Murray 
and Thackray in 1950%) and has been shown 
by ourselves and other investigators to be 
a common occurrence when sintering other 
“active” powders. An explanation of these 
effects could lie in the interaction of grain 
growth and sintering, but explanations in terms 
of grosser mechanisms such as bloating due to 
entrainment of gases in closed pores should not 
be lost sight of; dilatometric studies of the 
sintering process using heating rate as a variable 
would help to eliminate some of the possible 
explanations of the “‘over-sintering”’ of uranium 
oxides. 
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BOOK 


F. A. Rouen et A. A. Baur, Constitution of 
Uranium and Thorium alloys (Addison—Wesley, 
1958. 153 p.) 


Ce livre représente la 2éme édition d’un document 
(BMI 1000) déja bien connu dans les laboratoires de 
métallurgie nucléaire. Il faut lui souhaiter 
diffusion encore plus grande car il rassemble des 
renseignements plus récents, plus complets et la 
présentation a été améliorée. 

Les systémes sont présentés dans lordre alpha- 
bétique (en langue anglaise) des éléments d’addition; 
il n’y a plus de classement a part des systémes dont 
le diagramme d’équilibre n’a pas été établi et la 
consultation de louvrage s’en trouve facilitée. 

Tous les diagrammes comportent les deux échelles 
de concentration: atomes % et poids %, ce qui est 
particuliérement utile pour les alliages a base d’ura- 
nium et de thorium. Un bref exposé résume les données 
existantes sur la constitution des alliages et les 
structures cristallines des composés formés. Comme 
beaucoup de références sont classifiées ou non publiées, 
ce livre constitue une source unique d’informations, 
trés précieuse pour les chercheurs extérieurs a l] AEC 
ou a PUKAEA. 

Un tel ouvrage n’est jamais complet. Le grand 
intérét qu’il présente pour les diagrammes d’équilibre 
et les données cristallographiques fait regretter qu’il 
n’y ait pas un mot ni une référence sur les mécanismes 
et cinétiques de transformation, ni sur les principales 
propriétés des alliages. 

On peut regretter aussi que certaines références 
frangaises ou anglaises aient été omises, et qu’en 
général un plus grand poids soit accordé aux publi- 
cations de source américaine; ¢’est le cas, par exemple, 
du systéme U-Zr pour lequel Rough et Bauer repro- 
duisent le diagramme du Battelle Memorial Institute, 
alors quwune toute récente étude d’Argonne vient de 
confirmer le diagramme différent de Summers—Smith 
(Aldermaston-G.B.). 

La petite échelle des figures empéche souvent de 
distinguer la position des limites de solubilité; mais 
peut-étre un agrandissement eut-il donné une idée 
fausse de la précision des mesures. Et cela se rattache 
a notre derniére observation: 

On trouvera sur ces diagrammes beaucoup de lignes 
en pointillé, mais peut-étre n’y en a-t-il pas encore 
assez. En effet, méme lorsqu’un systéme a été étudié 
par plusieurs auteurs, il arrive fréquemment que 


une 
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chacun de ceux-ci se soit attaché plus spécialement a 
un domaine particulier de température ou de compo- 
sition et qu’ils aient tous reporté une ligne, un liquidus 
par exemple, tracé par un des premiers expérimen- 
tateurs utilisant des métaux de base de pureté 
insuffisante ou des techniques d’élaboration n’excluant 
pas toute contamination des alliages. Ainsi dans un 
diagramme, fruit de nombreuses études, il peut se 
trouver encore des lignes douteuses. 

C’est pourquoi, tout en insistant sur le grand | 
intérét de ce document, nous recommandons de 
Vutiliser seulement comme un catalogue pour une 
comparaison grossiére, ou comme un recueil de cons- 
tantes de provenances connues. Mais pour tout 
travail de détail il est indispensable de se reporter aux 
articles originaux cités par Rough et Bauer, pour 
connaitre les analyses des matériaux utilisés, les modes 
opératoires, la précision des résultats, etc.... Pour 
mieux démontrer cette nécessité, il suffit de signaler 
qu’en reportant sur un méme graphique toutes les 
données publiées concernant un diagramme trés 
étudié comme U-Mo ou U-Nb, on obtient un faisceau 
de lignes qui révéle une incertitude atteignant souvent 
3 at % ou 20°C, si on donne la méme valeur a tous 
les résultats. G. CABANE 


J. A. Lang, H. G. MacPuerson and F. Mastan 
(editors), Fluid Fuel Reactors. (Addison-Wesley 
Publishing Company, Reading, Pennsylvania, 
and London, 1958. 1008 pages. $11.50, 87s.). 


This is the largest (1008 pages) of the Geneva 
Presentation Volumes commissioned by the USAEC 
for the 1958 Atoms for Peace Conference. One may, 
at first sight, be puzzled by the fact that this volume 
is larger than that on Solid Fuel Reactors (864 pages) 
of which there has been much more operational 
experience. The clue lies in the foreword by Dr. Alvin 
Weinberg in which he discusses the two approaches to 
reactor technology—the mechanical engineering ap- 
proach with simplicity of design as its keynote and 
the chemical approach which leads to the concept of 
a reactor as a type of chemical plant utilizing a fissile 
solution as its working fluid. It is unfortunate that this 
comparison is not developed further in this book in 
view of the uncertainty in recent years of the value of 
pressing forward with the development of fluid fuel 
systems. This uncertainty is bound up with the reason 
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for the size of this volume, i.e. the very difficult 
problems arising from the design of plant for the 
generation and utilization of fission heat using a highly 
radioactive and corrosive fluid. The solution of these 
problems has involved a vast amount of research and 
development effort; this book provides a lucid and 
authoritative account of the results of this work. 
The concept of fluid fuel reactors dates back to the 
early days of atomic energy when Halban and 
Kowarski in 1940 performed experiments on sus- 
pensions of U30s in heavy water and also are reported 
to have suggested the use of a uranium-bismuth 


_ solution as a fuel. The development of homogeneous 


aqueous reactors (HAR) has been the most rapid of 
the various fluid fuel types, probably because the 
physics experiments involved comparatively little 
experimental development. A number of reactors have 
been operated successfully, the most recent being the 
5MW HRE-2 at Oak Ridge. In contrast only one fused 
salt reactor, the Aircraft Reactor Experiment (ARE), 
which had a very limited life, has operated and so far 
no liquid metal fuelled reactor (LMFR) has operated. 
This is partly a reflection of the relative degrees of 
experimental difficulty and partly of the effort applied 
reflected in the distribution of space in the book, 
between the three main types of system: HAR 566 
pages, ARE 136 pages and LMFR 244 pages. 

All three systems have the same ultimate objective 
—a breeder reactor, probably working on the Th282/ 
U3 cycle, incorporating integral chemical processing 
and using the fuel as the primary heat transfer fluid. 
Nearly all the systems considered have employed 
highly enriched fuel in thermal cores each using a 
different type of moderator (HAR-D20, ARE-BeFs, 
LMFR-Graphite). These factors lead to compact cores 
with highly rated fuel, from which it follows that 
processing rates will be relatively high and neutron 
fluxes in the core very high. In the case of HAR and 
LMFR the favoured design is for a two-zone system 
with a U3 solution core surrounded by, but physically 
separated from, a thorium slurry breeder blanket. The 
advantages of such systems over conventional solid 
fuel reactors are the absence of a radiation damage 
limitation on fuel burn-up, lack of expense of compli- 
cated fuel fabrication processes, continuous removal 
of fission products and continuous make-up of fuel 
concentration. The HAR system offers the best 
neutron economy but the other two types offer high 
temperature, low pressure operation. In every case 
corrosion difficulties are the outstanding disadvantage. 

In the HAR system the radiation-induced corrosion 
of the zirconium or Zircaloy-2 ““window”’ between the 
core and the blanket has proved to be a major obstacle. 
It is now thought that the radiation-induced phase 
change in the ZrOz protective film causing the film to 
disrupt is probably the main corrosion-accelerating 


process. This had led to detailed studies of corrosion 
under irradiation, and to the development of new 
zirconium alloys and. of titanium cladding as possible 
remedies. The problem has not yet, however, been 
fully solved and is one of the reasons for the consider- 
able diminution of interest in HAR systems over the 
past year or two. 

In the case of the fused salt reactor the corrosion 
problem seems to have been largely solved by the 
development of a new nickel-base alloy designated 
INOR-8 which combines the best characteristics of 
Inconel and Hastelloy B without their disadvantages. 
INOR-8 has excellent corrosion resistance to molten 
fluorides at temperatures well above those expected 
in a reactor. A great deal of work has been carried 
out on the fabrication and on the measurements of 
the properties of this alloy and is fully reported in 
this book. 

Again the main materials problem of the LMFR 
system has been the development of a corrosion 
resistant container material. In this case a somewhat 
different approach is adopted in that a mildly corrod- 
able steel (24 % Cr, 1% Mo) is used in conjunction 
with a corrosion inhibitor, zirconium, which forms a 
protective film of zirconium nitride and carbide on 
the steel surface. Extensive loop testing has developed 
this technique to a stage where it can be used with 
confidence. The study of the mechanism of the 
inhibition process, which also helps to protect the 
graphite moderator, is fascinating and it is a pity 
that it has not been described in this book in the 
detail which it deserves. 

Apart from the work on problems of specific interest 
to these reactors there is much that is of far wider 
interest and application and it is probably this aspect 
which will give the book a fairly wide appeal. One of 
the best examples is the work on the HAR breeder 
blanket. In the first place the data on the physical 
and chemical properties of aqueous suspensions of 
oxides were not very comprehensive nor indeed well 
interpreted when this study was started and much of 
general value has been learned of the preparation, 
stability and flow behaviour of aqueous slurries. As 
part of this work it has proved necessary to investigate 
a variety of methods of producing UOz and ThOs 
particles of controlled size and shape, which is again 
of general interest. 

The development of INOR-8 has been mentioned; 
here is another case where the material developed for 
a specific reactor application can find other uses in the 
chemical industry. 

Many general problems associated with the use and 
handling of molten metals were investigated in the 
course of the LMFR study. In particular the process 
of thermal gradient mass transfer and its prevention 
by soluble inhibitors has been studied in detail. The 
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work has also led to consideration of fundamental 
aspects of liquid metal solubilities. 

One could cite many other instances—development 
of impermeable graphite, the production of pure 
tantalum, remote maintenance of reactor components 
ete.—which illustrate the vast mine of useful infor- 
mation contained in this book. 

Anyone who, like the reviewer, has been associated 
with work on fluid fuel systems will, on reading this 
book, feel a sense of disappointment, if not frustration 
that, on reaching this stage of development, a large- 
scale cut has been made over the past 18 months in 
the effort on such systems on both sides of the 
Atlantic. To the reviewer’s knowledge only two teams 
in the world are currently pursuing integrated 
programmes on fluid fuel reactors; the LAMPRE 
plutonium-solution, tantalum-contained, fast reactor 
at Los Alamos and the Dutch UOs-gas suspension 
(fluidized bed) reactor which, of course, is not 
mentioned in this book. 

The fact that this book has been compiled with a 
separate editor and numerous contributors for each 
section, together with its length, mean that this is 
quite definitely a textbook of the type to which one 
refers for specific information rather than a coherent 
and easily readable account of fluid fuel reactor 
development. Moreover, it is a somewhat uncritical 
survey. Each section has been written by.an enthusiast 
in his special field. A section by an uncommitted 
reactor technologist comparing the relative merits of 
the systems described and their advantages and dis- 
advantages compared with solid fuel systems would 
have enhanced the value of the work. Nevertheless it 
will undoubtedly remain for many years the standard 
reference book on fluid fuel reactors. 

B. T. R. FROST 


SUMMARIES OF PHYSICAL RESEARCH IN METAL- 
LURGY, SOLID STATE PHYSICS AND CERAMICS. 
(USAEC, July 1959. TID-4005 (Pt-1, 5th 
Edition). 204 pages. Obtainable from Office of 
Technical Services, Department of Commerce, 
Washington 25, D.C., at $ 2.75.) 

This is a useful survey of unclassified research work 


in progress at Universities and in industrial labora- 
tories under USAKC sponsorship. About half of the 


150 research projects are concerned with fundamental 
physical metallurgy of very diverse materials; the 
other half deal with the technology of reactor materials 
and with radiation damage. The total range covered 
by the projects is very wide. 

Indexing is rather inadequate in that no page 
numbers are given and it is difficult to find a desired 
abstract quickly. In other respects the presentation 
is of a high standard. 


FUEL ELEMENTS CONFERENCE. (USAEC, August 


1950. TID-7559 (Part 1). 303 pages. Obtainable 
from Office of Technical Services, Department 
of Commerce, Washington 25, D.C., at $ 3.00). 


This is a record of a conference held in May 1958 in 
the U.S.A. 17 of the 43 papers were classified and are 
not included here; in particular the classified group 
includes most of those under the headings “Non- 
Destructive Testing of Reactor Fuels” and ‘Mobile 
Reactors’. Experience with some unconventional 
fuels, including U-ZrH and UO2-ThOs, is discussed, 
and several papers deal with problems of fabrication 
and design of fuel elements, including unusual forms 
such as tubular elements. 

A significant amount of space is devoted to uranium 
oxide. Authors agree that this fuel is very little liable 
to fission gas swelling, and one paper goes in detail 
into the problem of calculating the life of UOe fuel 
elements on the criterion that life is limited by the 
time required to build up a maximum premissible 
released gas pressure. This and several other papers 
are followed by searching discussions which are a 
valuable part of the proceedings. 

A final section consists of a number of brief papers 
in which experts evaluate future trends of fuel 
materials. This includes papers on U and U alloys, 
UOsz, dispersion fuels, fast reactor fuels, canning 
materials, reprocessing and radiation damage. The 
author of the last paper evidently felt something of an 
outsider and concludes plaintively: ‘‘Well, the other 
thing that I should like to anticipate is that some day 
there will be a fuel element conference and radiation 
damage will be the first thing on the program, not 
the last’’. 


R. W. CAHN 


